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The imperfections of many gases may be obtained by 
combining available heats of vaporization and vapor pres- 
sure data with the Clausius-Clapeyron equation. The ac- 
curacy of this method in some cases is comparable to that 
of good direct P.V.T. measurements and therefore may be 
useful in preliminary investigations. If P is the vapor 
pressure, T is the absolute temperature, AH is the heat of 
vaporization, V and V; are the molal volumes of the vapor 
and liquid phases: 

PV —AH B 


RT RQ—Vi/V)[dinP/d/T)] | V 





This equation determines the gas imperfection 8. The value 


of 8 calculated in this manner is compared with experi- 
mental second virial coefficients for HxO, NH3, He, Oz, A, 
CH,, propane, n-butane, n-heptane, ethene, propene, 
CH;Cl, CCl,F, CH;0H, CsHs, and Hg. Values of 6 are 
determined for HI, HBr, HCl, CH;NHo2, (CH;3)2NH, HCN, 
isobutane, tetramethylmethane, 2,2,4, trimethylpentane, 
C2N2, COS, PH;, HS, CH;Br, and C2H;OH. Evidence of 
hydrogen bond formation is shown for CH;0H, C2H;OH, 
CH;NHz2, and (CH3;)2NH. It is shown that HI and HCN 
contain a considerable fraction of instantaneously formed 
double molecules. The gas imperfection of benzene is per- 
fectly normal. The number of instantaneously formed 
double molecules is estimated for all of the above molecular 
species. 





HE imperfections of many gases may be 
obtained by using the Clausius-Clapeyron 
equation together with available experimental 
heats of vaporization and vapor pressure data. 
The accuracy in some cases is comparable to 
that of good direct P.V.T. measurements. In this 
way we can learn a great deal about the degree 
of association and the nature of the intermo- 
lecular forces in gases for which no direct P.V.T. 
measurements have ever been made. However, 
any conclusions based on this indirect method 
should be verified in other ways before they are 
regarded as scientific truth. 
The Clausius-Clapeyron equation is usually 
written : 


dP/dT=AH/T(V—V)). (1) 


Here P is the vapor pressure, T is the absolute 
temperature at which the boiling takes place, 
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AH is the heat of vaporization, and V and V, 
are the molal volumes of the vapor and the 
liquid, respectively. This equation may be ar- 
ranged into the form: 


PV ~AH 
RT R(A—V,/V)(d In P/d(1/T)}] 








(2) 


In elementary treatments the vapor is usually 
assumed to be ideal so that PV/RT =1 and the 
volume of the liquid is neglected in comparison 
with the volume of the vapor. When these 
assumptions are made, it is easy to integrate 
Eq. (2) to express the heat of vaporization 
directly in terms of the vapor pressure at two 
different temperatures. The values for the heats 
of vaporization obtained in this manner with 
good vapor pressure measurements are usually 
in error by around six percent. In this paper we 
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are attributing this deviation to the imperfection 
of the gas and making it serve a useful purpose. 
From the experimental measurements of vapor 
pressure we can calculate dln P/d(1/T). If in 
addition we know AH and V/V, the right-hand 
side of the equation is completely determined and 
we are left with an experimental value of 
PV/RT. 

It is convenient to express the imperfection of 
the gas in terms of a function, 6, defined by the 
equation : 

PV/RT=1+8/V. (3) 


We have tabulated the observed values of 8 for 
a number of substances. Large negative values 
correspond to association, and we can estimate 
the percentage of double molecules present. 

The equation of state of any gas may be 
written in the virial form: 


PV/RT=1+B(T)/V+C(T)/V?+---, (4) 


where B(T), C(T), ---, are called the second, 
third, ---, virial coefficients and are functions of 
the temperature but not of the pressure. In 
ordinary P.V.T. measurements one can deter- 
mine the virials separately by varying the 
pressure while keeping the temperature constant. 
This is not possible in our method since we can 
only consider one pressure for any particular 
temperature. However, most of the experimental 
vapor pressure measurements are made at pres- 
sures sufficiently low that the contribution of 
third and higher virials is negligible for a normal 
gas, and for these we can show that B=B(T). 
For associated gases where @ is exceptionally 
large it is probable that the third and higher 
virials appreciably affect the value of 8. 

The argument that B= B(T) for normal gases 
measured at a vapor pressure of around one 
atmosphere or less depends on the theorem of 
corresponding states. If we divide the second 
virial coefficient by the critical volume V, we 
get a reduced second virial coefficient B,(7/T.), 
which should be a universal function of the 
reduced temperature 7/7,. Similarly, we could 
define reduced third and higher virial coefficients. 
Hirschfelder, McClure, and Weeks! have shown 
that for a number of normal gases the experi- 
mental second virial coefficient agrees reasonably 


1J. O. Hirschfelder, F. T. McClure, I. S. Weeks, J. 
Chem. Phys. 10, 201 (1942). 


well with that which would apply if energy of 
interaction of the molecules satisfied a Lennard- 
Jones potential : 


E(r) =4En[(ro/r)? — (ro/r)°], (5) 


where E,,, the maximum energy of interaction 
between pairs of molecules on collision, is 


En, =0.73T Rk (6) 


and k is the Boltzmann constant. Now we can 
make use of the work of Montroll and Mayer? 
who express the third and fourth virial coeffi- 
cients for the Lennard-Jones gases in the form : 


C(T) = —8/3(B(T))?»2 


D(T) =6(B(T))*vs. (7) 


Here v2 and v3 are known functions of E,,/kT. 
The normal boiling point usually corresponds to 
values of T approximately 0.67, or 


En/kT =0.73T k/0.6T k= 1.2. (8) 


For this argument Montroll and Mayer com- 
puted that v»=0.019 and »v3=1.13. For normal 
gases, in this temperature range, |B/V| is 
always less than 0.05 so that 


|C/V2| <|(8/3)(0.05)2(0.019) | =0.0002 (9) 
and: |D/V*| <|(6)(0.05)9(1.15) | =0.001. 


Therefore, in identifying 8 with B(T) we are 
neglecting terms in PV/RT of the order of 0.001 
or making an error in B(T) of the order of 
22,000 0.001 of 22 cc per mole. This is about 
the limit of the accuracy of 8 which can be 
attained from experimental considerations. 

The ratio of double to single molecules present 
in the gas at any instant, N2/N; is estimated by 
the method developed by Hirschfelder, McClure, 
and Weeks.! In place of B(T) we take 6 and for 
the van der Waals d for the definitional domain 
of the double molecules we take 1.75 V.. Then the 
equilibrium constant for the formation of double 
molecules is: 


N2V 
K= =1.75V.—8. 


7 
Since we are usually dealing with a small degree 


2 According to a private communication from J. Mayer. 
The development of these higher virial coefficients was sub- 
sequently published by E. W. Montroll and J. E. Mayer, 
J. Chem. Phys. 9, 626 (1941). 





GAS 


TABLE I. 








—B 
calc. 


460 
287 


T°K 


373.16 
423.16 
473.16 196 
523.16 139 
573.16 97 
623.16 58 


233.16 
239.68 
243.16 
273.16 
293.16 
313.16 





561 
507 
467 
319 
250 
207 


516 
438 
443 
329 
265 
217 


18.2 164° 
20.4 724 
90.13 206° 
87.29 244! 
99.54 4090 
111.55 2038 


231.04 658" 
272.66 904 
371.51 13003 
169.40 428* 
225.35 575! 
248.94 527™ 
290.40 831" 


. propane 
n-butane 
n-heptane 
ethene 
propene 
CH;Cl 
CFCl; 

. CH;OH 298.16 
323.16 
337.86 
363.16 

, 383.16 


4190° 

2130 

1480 
905 
655 


1390P 
885 
746 

20502 


323.16 
353.16 
383.16 


. Hg 630.0 


Benzene 








* Reference 5. 

» Reference 9. 

© Martinez and Onnes, Comm. Phys. Lab. Univ. Leiden 156b (1922); 
Bichowsky and Rossini, Thermochemistry of Chemical Substances 
(Reinhold, 1936). 

43 Martinez and Onnes, Comm. Phys. Lab. Univ. Leiden 156b (1922); 
Eucken, Ber. d. Deut. Phys. Ges. 18, 4 (1916). 

¢ Cath, Comm. Phys. Lab. Univ. Leiden 152d (1918); Giauque and 
Johnston, J. Am. Chem. Soc. 51, 2300 (1929). 

f Born, Ann. d. Physik 69, 473 (1922); Frank and Clusius, Zeits. 
f. physik. Chemie 42, 395 (1939). 

® Henning and Stock, Zeits. f. Physik 4, 226 (1921); Frank and Clu- 
sius, Zeits. f. physik. Chemie 42, 395 (1939). 

» Kemp and Egan, J. Am. Chem. Soc. 60, 1523 (1938). 

‘ Aston and Messerly, J. Am. Chem. Soc. 62, 1917 (1940). 

i Smith, J. Research Nat. Bur. Stand. 24, 229 (1940); Pitzer, J. Am. 
Chem. Soc. 62, 1224 (1940). 

k Egan and Kemp, J. Am. Chem. Soc. 59, 1265 (1937). 

! Powell and Giauque, J. Am. Chem. Soc. 61, 2367 (1939). 

™ Messerly and Aston, J. Am. Chem. Soc. 62, 886 (1940). 
om” Garner, Doescher, and Yost, J. Am. Chem. Soc. 63, 3496 

© Young, Proc. Roy. Soc. Dublin 12, 374 (1910); Fiock, Ginnings, and 
Holton, J. Research Nat. Bur. Stand. 6, 881 (1931). 

» Dejardin, Ann. de Physique 12, 253 (1919); Fiock, Ginnings, and 
Holton, J. Research Nat. Bur. Stand. 6, 881 (1931). 
. "Menzies, Zeits. f. physik. Chemie 130, 90 (1927); W. Kurbatoff, 
ibid. 43, 104 (1903). 


of dissociation, to the first approximation : 
N,iRT 
P 


IMPERFECTIONS AND VAPOR 


PRESSURE 


So that: 


V. 


N2 PV, B 
—= 1.75 | 
Ni RT 


I. EXPERIMENTAL VERIFICATION OF 
THE METHOD 


Before relying on such an indirect method of 
determining the imperfections of gases it is 
obligatory to make a rather extensive comparison 
between our values of 6 and those of B(T) 
calculated from direct P.V.T. measurements. In 
the cases of NH; and H.O sufficient data are 
available to permit us to make the comparison 
over a considerable temperature range. For He, 
Oz, A, CH4y, CH;Cl, CCl3F, propane, n-butane, 
h-heptane, ethene, and propene we checked our 
values at a temperature in the neighborhood of 
the normal boiling point. In Table I, 8 is com- 
pared with the second virial coefficients obtained 
directly. Most of the values agree within the 
accuracy of the experimental data considered; 
only a few values are poor. 


1. H,O 


Our values of 6 for steam were determined 
from the empirical equation of Keenan and 
Keyes* for the vapor pressure as a function of 
the temperature and from the heats of vaporiza- 
tion of water determined by Osborne, Stimson, 
and Ginnings.* It was necessary to differenti- 
ate the vapor pressure equation to obtain 
d\n p/d(1/T). This is probably the largest single 
source of error in our method. The values of 8 
were compared to those of B(T) calculated from 
the equation of F. G. Keyes.’ The experimental 
values for steam are exceptionally accurate so 
that it is not surprising that the agreement 
between 6 and B(T) is good. 


2. NH; 


Osborne and Van Dusen® determined the heats 
of vaporization of NH; at various temperatures. 


3 J. H. Keenan and F. G. Keyes, Thermodynamical Prop- 
erties of Steam (John Wiley, 1936). 

4 Osborne, Stimson, and Ginnings, J. Research Nat. Bur. 
Stand. 23, 261 (1939). 

5 F. G. Keyes, J. Am. Chem. Soc. 60, 1761 (1938). 

6 N.S. Osborne and M. S. Van Dusen, J. Am. Chem. 
Soc. 40, 1 (1918). 
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Fic. 1. Reduced second virial coefficients as a 
function of the reduced temperature. 


Overstreet and Giauque’ developed an equation 
for the vapor pressure of NH; in terms of the 
temperature, from 199°K to 242°K. They also 
determined the heat of vaporization at 239.68°K. 
This work was done in process of determining 
the entropy calorimetrically. Cragoe, Meyers, 
and Taylor® developed an equation for the vapor 
pressure from —78°C to 70°C. We calculated 8 
from these data. 

Keyes® also gives an equation for the virial 
coefficient of NH3. This equation gives us values 
of B(T) to compare with our values of 8. In 
order to compare our accuracy with that of 


7R. Overstreet and W. F. Giauque, J. Am. Chem. Soc. 
59, 254 (1937). 

8 Cragoe, Meyers, and Taylor, Bur. Stand. Sci. Papers 
16, 1 (1920). 
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various direct P.V.T. determinations, we calcu- 
lated the second virial coefficient from the 
previously accepted data for NHz; as tabulated 
in a publication of the Bureau of Standards.® The 
differences between the two relatively accurate 
P.V.T. determinations of B(T) is of the same 
order of magnitude as that between 8 and B(T). 


3. Simple Non-Polar Gases 


The second virial coefficients for the simple 
non-polar gases were estimated on the assump- 
tion that the molecules obey a Lennard-Jones 
type of potential. For these calculations we used 
previously determined values of the molecular 
constants.’ These values agree reasonably well 
with our values of 8. 


4. Hydrocarbons and Substituted Methanes 


These were compared to values extrapolated 
from experimental data at high temperatures by 
means of equations developed recently.! The 
value of 8 for CH, at 99.54°K is obviously 
erroneous and probably due to small experi- 
mental errors. 


5. Methyl Alcohol and Benzene 


Eucken and Meyer"! have measured directly 
the second virial coefficients for a number of 
organic substances. They included methyl alcohol 
and benzene. Their values are compared to ours. 
The values of B(T) are calculated from their 
empirical equations. This amounts to a short 
extrapolation in both cases. At the higher tem- 
peratures the agreement is satisfactory, i.e., 
within 25 percent. 


6. Hg 


Our value is of doubtful accuracy since the 
heat of vaporization was determined by W. 
Kurbatoff in 1903 by a rather crude method. 
However it is interesting to note that this value 
is in line with the fact that mercury vapor con- 
tains a considerable fraction of double molecules. 
Actually mercury contains such a large fraction 


9 Circular of the Bureau of Standards, No. 142 (1925). 

10J. O. Hirschfelder and W. E. Roseveare, J. Phys. 
Chem. 43, 15 (1939). 

1 A, Eucken and L. Meyer, Zeits. f. physik. Chemie 5b, 
452 (1929). 
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of Hg that its spectrum is well known and 
spectroscopic constants have been calculated for 
it. Hildebrand, Wakeham, and Boyd” deter- 
mined the potential energy curve for Hg2 from 
this spectroscopic data and Robert Bonner" has 
estimated from it the second virial coefficient as 
a function of temperature. The large deviation 
between 6 and the calculated B(T) may or may 
not be significant. Further experimental work on 
the heat of vaporization and vapor pressure of 
mercury would be very desirable. 


II. APPLICATION OF THE METHOD 


Since this method appears to be valid whenever 
the experimental values of the vapor pressure and 
heats of vaporization are reasonably accurate, it 
is interesting to examine data for a number of 
different molecular species for which there are no 
direct P.V.T. data. In Fig. 1 we have plotted 
the values of —8/V. versus T/T, for a num- 
ber of substances. V. and 7, are the critical 
volume and temperature, respectively. If the 
equation of corresponding states were valid, all 
of these points should lie on one curve. Actually 
the points show no such regularity. The reason 
for this is that the particular molecules studied 
here were chosen because of certain unusual 
properties, i.e., tendency to form hydrogen 
bonds, asymmetrical shape, large dipole mo- 
ments, and the like. The reason for the deviations 
from the corresponding states equation is ex- 
plained for the individual molecules. In Table II, 
our results are summarized and references given 
to all of the experimental data. We have in- 
cluded a column for u?2/7T.V. (where yw is the 
dipole moment) since the value of this reduced 
dipole moment should govern the deviation from 
the corresponding states equation due to dipole 
effects. The fact that this is not true for molecules 
forming hydrogen bonds shows that the forma- 
tion of hydrogen bonds must be considered quite 
separately from normal dipole behavior. Un- 
doubtedly many of our values of 8 are unreliable 
because of errors in the experimental heats of 
vaporization and vapor pressure, but it is amus- 
ing to see what implications could be drawn if we 
assume that all of our values are correct. 





® Hildebrand, Wakeham, and Boyd, J. Chem. Phys. 7, 
1094 (1939). 
8 Robert Bonner, unpublished. 


Where possible the critical constants were 
taken from the International Critical Tables. In a 
few cases where no data on the critical volume 
were available, we estimated the value from the 
critical temperature and pressure and a value 
of pv/RT at the critical conditions equal to that 
of a similar substance. The dipole moments are 
those recommended by Smyth."! 











TABLE II. 
N2/N1 
w?/TeVe (ati 
T°K -s T/Te —B/Ve (X10~*) atmos.) 
HI 237.75 13002 = .561 11.50 30 .077 
HBr 206.38 767> = .568 8.14 1.78 055 
HCl 188.07 446¢ 579 5.19 3.80 .039 
CH3NHe2 266.84 5184 .620 4.11 2.79 .034 
(CH3)2NH 280.04 868e .640 4.74 1.15 052 
HCN 298.80 1910! 654 14.16 11.4 .088 
Isobutane 261.44 695¢ 
Tetramethyl- 
methane 278.59 10905 
2,2,4, 
Trimethyl- 
pentane 372.33 1480: 
CoNe 251.95 695i 628 4.65 0 .046 
Cos 214.03 470k .566 3.45 82 .040 
222.87 565! .590 4.15 044 
PH; 185.35 379m ~=.572 3.35 83 .038 
185.72 1432 = .573 1.26 .032 
HS 188.7 410° ~=.505 4.45 2.63 .037 
212.77 1740P = .570 18.9 .109 
CH3Br 276.66 6604 
CH;0H 298.16 4190" ~—.581 35.6 4.67 .180 
323.16 2130 .630 18.1 .088 
337.86 1480 .658 12.6 .061 
363.16 905 -708 7.68 .037 
383.16 655 .747 5.56 .027 
C:H;OH 323.16 2040" .626 12.2 3.35 .087 
343.16 1470 665 8.82 .063 
363.16 1020 .703 6.09 .044 
383.16 715 .742 4.28 .032 
Benzene 323.16 1390s 575 5.43 0 .069 
353.16 885 .639 3.45 .046 
383.16 746 .682 2.90 .038 








The first reference refers to the source of data on the vapor pressure 
of the material. The second refers to the source of data on the heat of 
vaporization. Where only one reference is given, this contains data on 
both. 

a Henglein and Roth, Zeits. f. Physik 18, 64 (1923); Giauque and 
Wiebe, J. Am. Chem. Soc. 51, 1441 (1929). 

b Henglein and Roth, Zeits. f. Physik 18, 64 (1923); Giauque and 
Wiebe, J. Am. Chem. Soc. 50, 2193 (1928). 

¢ Henning and Stock, Zeits. f. Physik 4, 226 (1921); Giauque and 
Wiebe, J. Am. Chem. Soc. 50, 101 (1928). 

4 Aston, Siller, and Messerly, J. Am. Chem. Soc. 59, 1743 (1937). 

e Aston, Eidinoff, Forster, J. Am. Chem. Soc. 61, 1539 (1939). 

f Giauque and Ruehrwein, J. Am. Chem. Soc. 61, 2626 (1939). 

« Aston, Kennedy, and Schumann, J. Am. Chem. Soc. 62, 2059 (1940). 

bh Aston and Messerly, J. Am. Chem. Soc. 58, 2355 (1936). 

i Smith, J. Research Nat. Bur. of Stand. 24 229 (1940); Pitzer, J. 
Am. Chem. Soc. 62, 1224 (1940). 

i Ruehrwein and Giauque, J. Am. Chem. Soc. 61, 2940 (1939). 

k Kemp and Giauque, J. Am. Chem. Soc. 59, 79 (1937); Frank and 
Clusius, Zeits. f. physik. Chemie 42, 395 (1939). 

! Kemp and Giauque, J. Am. Chem. Soc. 59, 79 (1937). 

m Stephenson and Giauque, J. Chem. Phys. 5, 149 (1937). 

» Henning and Stock, Zeits. f. Physik 4, 226 (1921); Frank and 
Clusius, Zeits. f. physik. Chemie 42, 395 (1939). 

© Giauque and Blue, J. Am. Chem. Soc. 58, 831 (1936); Frank and 
Clusius, Zeits. f. physik. Chemie 42, 395 (1939). 

Pp Giauque and Blue, J. Am. Chem. Soc. 58, 831 (1936). 

4 Egan and Kemp, J. Am. Chem. Soc. 60, 2097 (1938). 

r Young, Proc. Roy. Soc. Dublin 12, 374 (1910); Fiock, Ginnings, 
and Holton, J. Research Nat. Bur. Stand. 6, 881 (1931). 

8 Dejardin, Ann. d. Physique 12, 253 (1919); Fiock, Ginnings, and 
Holton, J. Research Nat. Bur. Stand. 6, 881 (1931). 


4C, P. Smyth, Dielectric Constant (Chemical Catalog 
Company, New York, 1931). 











1. HI, HBr, HCl 


For this series HI has by far the largest value 
of —8/V. in spite of the fact that its dipole 
moment is the smallest. Apparently HI is really 
made up of a considerable fraction of Hels. 


2. CH;NH:; and (CH;),.NH 


Dimethylamine has a large value of —8/V. 
and seems to be largely associated. Methylamine 
has a somewhat smaller value of —8 but still 
contains a considerable fraction of dimers. These 
results are not surprising since the amines and 
particularly the diamines have a tendency to 
form hydrogen bonds. A straight dipole effect 
should give methyl amine the larger value of 
—B/V. so that its effect alone should give the 
opposite order. While the above results are 
suggestive, they are not conclusive since they 
might arise from experimental errors. 


3. Methyl and Ethyl Alcohol 


Both methyl and ethyl alcohols have very 
large values for —8/V. which would indicate 
that they contained a considerable fraction of 
dimers in the gas phase. These results are in 
keeping with the tendency for the alcohols to 
form hydrogen bonds. Considerable work has 
been done on the study of hydrogen bonds 
formed in liquid alcohol. 


4. Benzene 


The values of —8/V. for benzene seem per- 
fectly normal and would agree with the corre- 
sponding states equation for the second virial 
coefficients of the hydrocarbons.; Benzene mole- 
cules are shaped like plates and interfere with 
each other’s rotation in the liquid. (This argu- 
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ment is used to explain why the benzene made, 
from deuterium has a greater vapor pressure than 
ordinary benzene.) On this account, we should 
expect an unusually small number of double 
molecules. Direct P.V.T. determination of the 
second virial coefficient of benzene would be 
particularly useful for theoretical interpretations. 


5. HCN 


It appears that HCN gas contains an appreci- 
able fraction of instantaneously formed double 
molecules. This agrees with the conclusion 
reached by Giauque and Ruehrwein" that HCN 
vapor contains straight chain polymers of in- 
definite length. 


CONCLUSION 


1. The Clausius-Clapeyron equation combined 
with accurate experimental values of the heats 
of vaporization and vapor pressures suffices to 
make reasonable estimates of the second virial 
coefficients for a number of substances for which 
no direct P.V.T. data is available. 

2. The most interesting examples are those 
molecules which have a tendency to form hydro- 
gen bonds. Our calculations for the fraction of 
dimers present in HI, dimethylamine, methyl 
and ethyl alcohol, and HCN should be sub- 
stantiated or refuted by infra-red spectroscopic 
observations of the vapors. 

3. This method is particularly suitable for use 
with those compounds whose absolute entropies 
are determined with great precision by calorimet- 
rical methods and accurate heats of vaporization 
and vapor pressures are determined as by- 
products. 


16 W. F. Giauque and R. A. Ruehrwein, J. Am. Chem. 
Soc. 61, 2626 (1939). 
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The Photolysis of Dimethylhydrazine' 


WILLIAM L. Kay AND H. AusTIN TAYLOR 
Department of Chemistry, New York University, New York, New York 


(Received June 11, 1942) 


The ultraviolet absorption spectrum of dimethylhydrazine has been shown to exhibit a 
continuum extending from 2800A to beyond 2000A. Despite the absence of structure the 
continuum is probably experimental and indicates predissociation. The photolysis using 
various light sources, at several pressures and wave-lengths, has been studied. The over-all 
quantum yield is 0.3. The principal initial rupture of the molecule involves the scission of a 
hydrogen atom. Analytically three molecules of the hydrazine give one molecule each of He, 
Ne, CHs, NH3s, CH3sNH2, and (CH;N : CH2)2. The low quantum yield is accounted for by 
addition of hydrogen to the radical produced in the initial split. 


HE reaction of atomic hydrogen with azo- 

methane was shown by Henkin and Taylor? 
to yield s-dimethylhydrazine alone at room tem- 
perature. At 110°C there was, in addition, a 
rupture producing methane and methylamine. 
Whether this rupture was solely of the azo- 
methane, the dimethylhydrazine, or of both 
could not be clearly deduced. The photochemical 
decomposition of hydrazine has been studied but 
no work has been reported on s-dimethylhydra- 
zine. The absorption spectrum and over-all 
quantum yield has therefore been measured at 
various wave-lengths. 


EXPERIMENTAL 
Materials 


Dimethylhydrazine prepared from the dihy- 
drochloride by treatment with pellets of potas- 
sium hydroxide was drawn through a drying 
column filled with fused barium oxide and con- 
densed in a trap cooled with dry ice. The liquid 
was freed from dissolved gases by repeated 
evaporation and condensation under vacuum. 
The boiling point, determined by a micro method, 
was 80.0°C in agreement with that reported by 
West,? 80.5°C. The refractive index at 20°C 
using a sodium D lamp was 1.4204; West reports 
1.4201. The density was 0.826 as against West’s 
value 0.827. 


Abstract from a thesis presented in partial fulfillment 
of the requirements for the degree of Doctor of Philosophy 
at New York University, March, 1942. 

(1940) Henkin and H. A. Taylor, J. Chem. Phys. 8, 1 
: W. West and R. B. Killingsworth, J. Chem. Phys. 6, 1 
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Propylene used in the tests was taken from a 
commercial tank and repeatedly fractionated ‘n 
vacuum from dry ice to liquid nitrogen tempera- 
tures. Analysis showed the presence of 5 percent 
propane and 94 percent unsaturate. 

Uranium oxalate was prepared by cooling a 
hot solution of reagent quality uranium nitrate 
and oxalic acid. The resulting crystals were 
thoroughly washed, dried overnight in vacuum 
and then in air for three hours at 110°C. Oxalic 
acid and sodium oxalate were recrystallized 
several times, dried overnight in vacuum and 
finally in air at 105°C for three hours before use. 


Apparatus 


The reaction cell of fused quartz was a cylinder, 
20-cm inside length and 2.5-cm inside diameter. 
It was joined through a ground glass joint to a 
Pyrex manifold evacuated by a mercury diffusion 
pump and oil pump. The cell was removed for 
cleaning prior to each experiment. The cell was 
connected to a manometer and a series of traps 
through which the products of the reaction could 
be withdrawn by a mercury piston and delivered 
to a gas buret for measurement of volume or 
for collection over mercury for subsequent 
analysis. A small removable trap between the cell 
and mercury piston was filled with yellow mer- 
curic oxide prior to each experiment. 


Radiation 
Various light sources were used ir the work: 


Arc 1: A “cold” mercury vapor lamp operated 
from a transformer drawing 110 volts in the 
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primary and delivering 15 volts at the secondary. 
Arc II: A “hot” mercury arc by the Hanovia 
Company of the Alpine sunlamp type. 

Arc III: A General Electric Type H4 lamp with 
the outer Pyrex envelope removed, operating 
from a transformer drawing 2 amp. in the 
primary at 120 volts. 

Arc IV: A low pressure vapor lamp used for the 
quantum yield determinations. It was wound asa 
helical spiral of 10 turns with an internal di- 
ameter of 30 mm and operated from a 0.265-kva 
transformer. 

Spark: Several runs were made using a condensed 
spark between electrodes of Zn, Cd, or Al. The 
power was furnished by a 10-kva transformer 
drawing 30 amp. at 250 volts and delivering 
10,500 volts on the secondary. A Pyranoil 
capacitor of 0.0375 uf was placed parallel with 
the transformer. The electrodes were metal 
collars 10-cm diameter and 5 mm thick with 
V-shaped edges mounted on 6-cm disks carried 
by water-cooled shafts rotating at about 80 r.p.m. 
The electrodes ran at right angles to each other 
and the adjustable gap between them was about 
4 mm. The set-up is similar to that used by Wiig 
and Kistiakowsky.‘ 


Absorption Measurements 


The absorption spectrum of dimethylhydrazine 
was determined using a Hilger E-2 spectrograph. 
The light source was a hydrogen discharge tube 
supplied by the Hanovia Company operating at 
about 2500 volts and 0.3 amp. Eastman 50 and 
ultraviolet sensitized plates were used. The in- 
tensities of the photographed spectra were meas- 
ured on a Kipp recording microphotometer. 
Three absorption cells were used, of 10, 20, and 
120-cm length. 


Procedure 


The system being thoroughly evacuated (10-° 
mm or less on a McLeod gauge), the dimethyl- 
hydrazine, freshly pumped out at dry ice tem- 
perature, was allowed to warm up and fill the 
reaction cell to the desired pressure which was 
always less than the saturated vapor pressure at 
room temperature. Exposure to the arc or spark 


‘E. O. Wiig and G. B. Kistiakowsky, J. Am. Chem. Soc. 
54, 1806 (1932). 


KAY AND H. A. 
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which had been permitted to reach a steady state 
while shielded was then made. Pressure changes 
during the run were followed on the manometer ; 
the final pressure was read after the light had 
been turned off for several minutes. The products 
of reaction were drawn by the mercury piston 
through a trap in liquid nitrogen and their vol- 
ume measured in the buret at atmospheric 
pressure. A second fractionation was made using 
melting sec-butyl chloride at —131°C and a 
third from dry ice at —78°C. The gases were 
collected over mercury and analyzed in a Blacet 
and Leighton micro gas analyzer.® 

The gas drawn from dry ice was condensed in 
the yellow mercuric oxide trap to remove am- 
monia after its volume had been measured. The 
trap and its contents were allowed to warm to 
room temperature and again cooled several times 
until no further decrease in volume occurred. 
The residual gas was withdrawn and tested with 
an acid bead for complete absorption and taken 
to be methylamine. 

The gas condensed in the dry ice trap was 
warmed to room temperature and returned to 
the reaction cell where a measurement of the 
pressure permitted the calculation of its volume. 
This was taken as the undecomposed dimethyl- 
hydrazine. 

Quantum yield determinations were made 
using a uranyl oxalate actinometer, 0.01. uranyl 
oxalate and 0.05M oxalic acid. Measurements 
were made using Arc IV, a fused quartz lens 
serving to focus the light through an iris dia- 
phragm on the cell, and, since the lens only 
transmitted down to 2400A, it filtered out the 
1850A radiation. The actinometer vessel was of 
fused quartz 4 cm in length with planar end 
plates 2.5 cm in diameter. The solution was 
vigorously stirred during photolysis by a plati- 
num ribbon. A blank of the same solution was 
protected from light and titrated with permanga- 
nate at the same time-as the photolyzed solution. 
Two runs were made with the actinometer placed 
behind the empty reaction cell and then behind 
the cell filled with dimethylhydrazine. The 
normal procedure, however, was to make a run 

5F. E. Blacet and P. A. Leighton, Ind. Eng. Chem., 
Anal. Ed. 3, 266 (1931); 5, 272 (1933); 6, 344 (1934); 9, 
44 (1937). 


6G. S. Forbes and L. J. Heidt, J. Am. Chem. Soc. 56, 
2363 (1934). 
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Fic. 1. Absorption coefficients of dimethylhydrazine. 


with dimethylhydrazne and then replace the 
reaction cell by the actinometer vessel. 

In focal isolation the full light from the zinc 
spark was focused by two crystal quartz lenses 
and two small diaphragms each at optimum 
position. The lenses were plano-convex and 
carried stops of black paper mounted centrally 
on their plane faces. A small quartz spectrograph 
was used to determine the correct focus for the 
desired lines. Spectrograms showed the zinc 
lines 2030, 2060, 2100, and 2140A could be iso- 
lated as a group. 

In determining the absorption spectrum, ex- 
posures were first made for various times through 
the empty cell and then through the filled cell 
for the same times. This was repeated for various 
pressures of dimethylhydrazine. A copper spec- 
trum was used for comparison. 


RESULTS 
Absorption Spectrum 


The absorption spectrum of dimethylhydrazine 
was determined at pressures ranging from 0.1 to 
50 mm in cells 10, 20, and 120 cm in length. 
The spectrum showed continuous absorption 
commencing at 2800A and continuing through 
the limit of the plate. Spectrograms taken using 
the higher resolution of the Hilger E-1 instru- 
ment showed no structure either on the plate or 
in the photometer tracings. 

The absorption coefficient k, given by I/Io 
= 10-*«' where c is the concentration in moles per 
liter and J is the length in cm, was calculated 
using the characteristic curve for the plate de- 
termined from the photometer tracings of the 
exposures through the empty cell. In Fig. 1 data 
are presented for 4.2-mm pressure in the 120-cm 
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cell. It is clear that at 2537A absorption will be 
effectively complete under the experimental con- 
ditions of the photolysis. 


Decomposition Products 


The products of photolysis were shown to con- 
sist of hydrogen, methane, nitrogen, ammonia, 
methylamine, and a liquid. The first three, not 
condensed by liquid nitrogen, were analyzed by 
combustion with oxygen. Hydrogen was deter- 
mined in a separate sample using a cuprate bead. 
Nitrogen was obtained from the difference be- 
tween the total sample and the hydrogen and 
methane found. 

The volume of the sec-butyl chloride fraction 
was never more than 1/15 the volume of the 
liquid nitrogen fraction. Analyzed by combustion 
it consisted principally of methane although in 
runs using the hot arc placed close to the reac- 
tion cell, when the amount of dimethylhydrazine 
decomposed was large, small amounts of ethane 
were found. Introduction of a fuming sulfuric 
acid bead into this fraction before combustion 
caused a contraction of 3—4 percent. This might 
be caused by small amounts of ammonia being 
drawn over. 

The presence of ammonia and methylamine in 
the fraction drawn from dry ice was established 
in several ways. In the first place Tanret’ says 
that the ammonia reaction of Nessler’s reagent is 
not given by methyl or ethylamine. Charitschkoff® 
claims that the amine hydrochlorides which are 
decomposed by aqueous potassium iodide give 
the reaction. It was found in synthetic methyl- 
amine-ammonia mixtures that with appreciable - 
amounts of methylamine a precipitate resulted. 
If the solution is greatly diluted, however, it is 
possible to determine the ammonia colorimet- 
rically. 

An aqueous solution of the gases drawn from 
dry ice gave a yellow precipitate with Nessler’s 
reagent. A solution of pure methylamine gave a 
much lighter colored precipitate which darkened 
on addition of ammonia and could be matched 
with the former. Dimethylamine gives a white 
precipitate which even on addition of large 
amounts of ammonia could not be made to 


7C. J. Tanret, J. pharm. chim. 5, 28, 498 (1893). 


8K. W. Charitschkoff, J. Russ. Chem. Soc. 39, 230 
(1907). 
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TABLE I. Products of photolysis using mercury arcs. Gases are in cc at N.T.P. 




















Expt. number 28 29 34 36 39 40 41 49 50 52 
Light, arc and dist. in cm II, 3 II, 3 Cdarc I,6 IlI,*8 I,*8 I,8 IV, 8 IV, 8 IV, 8 
Initial (CH3;NH)- 4.85 5.23 4.00 4.40 4.55 1.10 3.92 0.95 1.70 0.35 
Time (min.) 140 22 70 30 53 325 25 60 61 60 
Press. increase (mm) 20.6 8.5 15.0 14.0 6.2 24 3.2 2.0 3.1 0.7 

2 1.74 0.72 1.32 1.05 0.45 0.22 0.27 0.17 0.32 0.10 
CH, 1.59 0.50 1.02 1.00 0.32 0.15 0.16 0.12 0.22 0.06 
C2He 0.20 0.05 0.10 0.08 0 0 0 0 0 0 
Ne 1.92 0.49 1.32 1.02 0.35 0.22 0.16 0.17 0.32 0.11 
NH; 0.55 0.46 0.38 0.50 0.30 0.16 0.30 0.05 0.10 0.05 
CH;NH, 0.45 0.40 0.35 0.37 0.25 0.01 0.27 0 0 0 
(CH3;NH)2 used 2.80 1.46 2.50 2.00 1.10 0.49 0.75 0.22 0.50 0.20 
% decomposed 58.0 28.0 62.0 45.5 23.2 44.5 19.1 23.2 29.4 57.0 

* Crystal quartz lens used. 

TABLE II. Dimethylhydrazine quantum yields using Arc IV. Gases are in cc at N.T.P. 

Expt. number 46 47 48 48A 51 53 53A 54 55 56 
Initial (CHs;NH).2 4.17 4.03 3.90 3.50 3.20 4.03 2.50 4.13 4.02 4.03 
Time (min.) 60 60 60 60 60 10 10 90 5 60 
Press. increase (mm) 5.0 8.4 tA 6.4 6.6 1.2 1.5 14.0 0.5 6.3 
Hz 0.47 0.92 0.74 0.60 0.60 0.17 0.14 0.95 0.10 0.52 
CH, 0.25 0.58 0.46 0.39 0.38 0.07 0.05 0.74 0.04 0.32 
Ne 0.24 0.57 0.46 0.40 0.42 0.12 0.11 0.78 0.08 0.31 
NH; 0.30 0.49 0.41 0.40 0.35 0.10 0.12 0.37 0.08 0.32 
CH;NH: 0.25 0.37 0.35 0.35 0.25 0.08 0.10 0.20 0.04 0.26 
(CH;NH)>2 used 1.15 1.70 1.40 1.20 1.22 0.30 0.25 1.87 0.14 1.10 
¢ on (CH;NH)>2 0.28 0.32 0.33 0.29 0.29 0.42 0.35 0.30 0.40 0.26 
¢@ on H 0.23 0.34 0.35 0.29 0.29 0.47 0.39 0.30 0.60 0.26 
% decomposed 27.6 41.2 36.0 34.3 38.1 7.4 10.0 45.2 3.5 27.3 








TABLE III. Products of photolysis using spark sources. 








Expt. number 33 33A 42 43 44 
Light Full Full Full Zn Full 
Zn Zn Zn F.Q. Cd 
Init. (CH;NH)> 440 4.35 4.35 3. 3.86 
Time (min.) 105 105 161 450 143 
Press. increase (mm) 8.3 8.4 1.5 0 8.7 
2 0.54 0.56 0.15 0.04 0.58 
CH, 0.35 0.35 0.05 0.01 0.34 
Ne 0.33 0.36 0.09 0.02 0.38 
NHs; 0.42 0.47 0.15 0.04 0.45 
CH;NH,2 0.40 0.44 0.110 0O 0.40 
(CH;NH)2 used 1.20 1.20 0.35 0.10 1.25 


Full 
Cd 
4.25 
140 
8.8 
0.60 
0.35 
0.36 
0.43 
0.40 
1.12 


67 


444A 45 57 58 58A 59 65 66 

Full Zn Full Full Zn Zn Zn Full 
Zn F.Q. Al Al FI. F.Q. FI. Al 
440 3.49 2.96 3.50 3.68 3.70 3.80 3.90 
150 310 194 200 300 340 £780 60 
7.5 5.9 7.3 6.8 3.3 3.0 3.5 2.0 
0.51 0.22 0.54 0.52 0.11 0.24 0.20 0.22 
0.30 0.12 0.33 0.32 0.05 0.12 0.10 0.10 
0.36 0.20 0.36 0.33 0.07 0.20 0.15 0.15 
0.41 0.27 0.27 0.30 0.18 0.25 0.10 0.17 
0.41 O.25 0.21 0.23 0.06 0.18 0.10 0.10 
1.09 0.50 1.00 1.02 0.24 0.60 0.48 0.40 








F.I., focal isolation; F.Q., fused quartz lens (>2400A); gases are in cc at N.T.P. 


match the precipitate from a run. Equal amounts 
of dimethyl and methylamines and ammonia 
gave a lighter color and only with traces of 
dimethylamine could a color match be obtained. 
Dimethylhydrazine gives a brown precipitate 
with Nessler’s reagent which is easily distin- 
guishable. 

The presence of a primary amine was proved 
by the pink color of the Rimini test using acetone 
and sodium nitroprusside. The absence of a 
secondary amine was proved by the failure of 





the Simon test using acetaldehyde and sodium 
nitroprusside. 

The separation finally adopted was based on 
the Francois® observation that ammonia forms a 
compound with yellow mercuric oxide while 
methylamine does not. Tests on mixtures con- 
firmed the observation. The mercuric oxide 
should be kept evacuated since exposure to 
mercury vapor was found to reduce its effective- 
ness. This method of analysis was checked 


9M. Francois, Comptes rendus 144, 857 (1907). 
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against that used by Emeleus and Jolley” of 
distilling the mixture of methylamine and am- 
monia in vacuum from — 118°C to liquid nitrogen 
temperature, the ammonia being volatile. The 
two methods agreed reasonably closely. Com- 
bustion of the material not absorbed by mercuric 
oxide showed it to be methylamine. All tests for 
the presence of cyanide, possibly as hydrogen 
cyanide, failed. 

During photolysis small amounts of a brown 
liquid formed on the sides of the reaction cell. 
After removal of the gases after an exposure, 
this liquid was found to be soluble in water, the 
solution being alkaline to methyl red. The solu- 
tion reduced Fehling’s solution in the cold while 
dimethylhydrazine only reduced it on warming. 
Methyl methyleneimine would hydrolyze to form- 
aldehyde and thus reduce Fehling’s solution in 
the cold. 

Table I presents a summary of the experiments 
using various mercury arcs under varying con- 
ditions. 


Quantum Yield 


Since with Arc IV ninety-eight percent of the 
ultraviolet radiation is concentrated in the 2537A 
resonance line and the fused quartz lens elimi- 
nated the 1850A line, the amount of radiation at 
2967, 3126, and 3132A which could be absorbed 
by uranyl oxalate must be extremely small. For 
complete absorption, as at 2537A, by both uranyl 
oxalate and dimethylhydrazine, the quantum 
yield of dimethylhydrazine is given by the ratio 
of the moles of the hydrazine decomposed to Io, 
the average number of quanta per sec. entering 
the actinometer solution; this is given by the 
relation Ib =NM/¢t where N is the Avogadro 
number, M is the number of moles of uranyl 
oxalate decomposed in ¢ sec., and @ is the 
quantum yield, given by Brackett and Forbes" 
as 0.63+0.03. 

Quantum yields were also calculated on the 
basis of atomic hydrogen by considering each 
hydrogen molecule found by analysis as repre- 
senting two hydrogen atoms from two dimethyl- 
hydrazine molecules. Table II presents a sum- 
mary of the data obtained. 


(1938) J. Emeleus and L. J. Jolley, J. Chem. Soc. p. 1612 
mn F. P. Brackett and G. S. Forbes, J. Am. Chem. Soc. 
55, 4459 (1933). 





It is apparent that the over-all quantum yield 
at 2537A is 0.3. 

From the data in Tables I and II it may be 
seen that the amount of dimethylhydrazine de- 
composed being a linear function of the time at 
a given initial pressure and intensity of illumina- 
tion is proportional to the total light absorbed at 
least up to 40 percent decomposition. In Fig. 2 
data from Tables I and II are presented graphi- 
cally as the volume of product formed per unit ~ 
amount dimethylhydrazine decomposed against 
the percentage decomposition. The small volumes 
of products and of dimethylhydrazine decom- 
posed for runs at low pressures justify their 
exclusion from the graph. For the same reasons 
runs with 10 percent decomposition or less are 
not reliable. The hydrogen, methane, and nitro- 
gen values are all high under these conditions 
suggesting that the error lies in the value of the 
hydrazine decomposed rather than in the gas 
analysis, since then the nitrogen, determined by 
difference, should be low for high methane and 
hydrogen values. 

A material balance calculated for the more 
reliable runs in Tables I and II, of the differences 
between the hydrazine decomposed and the C, N, 
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and H found in the products, expressed as a C» 
compound since carbon appears in a minimum 
number of products, gave CeHs,.1+0.4N1.1+0.1- 
The lowest value for H in any run was 3.5 and 
the highest 5.7, while N varied between 0.6 
and 1.4. 


Effect of Wave-Length 


The products of photolysis using the spark 
between electrodes of various metals are pre- 
sented in Table III. Runs were made using the 
zinc spark with and without the fused quartz 
lens between cell and source to study the effect 
of the zinc spectrum above 2400A as compared 
with the full spectrum. No significant variations 
are observable. This was also true when the zinc 
spectrum between 2030 and 2140A, obtained by 
focal isolation, was used. The aluminum spec- 
trum is rich in lines below 1900A, a region in 
which both ammonia and methylamine are 
known to decompose rapidly. This would account 
for the low values of these products found in 
runs with aluminum electrodes. 


Effect of Propylene 


Since propylene does not readily form active 
intermediaries or enter into a chain reaction it 
has been used to detect the presence of hydrogen 
atoms in the photolysis of hydrazine and amines. 
Bamford” found a large reduction in the amount 
of hydrogen and the production of propanes and 
isomeric hexanes. Several runs were made in 
presence of added propylene. The data are given 
in Table IV in which the values of propylene 
added, residual propylene, and propane have all 
been corrected for the propane originally present 
in the added propylene. It was not possible to 
determine the hexanes directly ; the amount was 
estimated as half the difference between the 
propylene used and the propane formed. The 
fractionation of the product gases from these runs 

‘was made as before. The sec-buty! chloride frac- 
tion was treated with fuming sulfuric acid to 
remove propylene. Combustion of the residue 
gave a carbon value slightly less than that for 
propane, probably due to dissolved methane. 
The dry ice fraction was treated as before to 
remove ammonia and methylamine, and then 


2 C, H. Bamford, J. Chem. Soc. p. 17 (1939). 
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TABLE IV. Products of photolysis in presence of propylene. 
Gases are in cc at N.T.P. 








Expt. number 37 38 60 61 63 
Light HgArcIl HgArcI FullZn ZnF.Q. Full Zn 
Initias (CH3NH):2 4.15 4.77 3.24 3.70 3.89 
Time (min.) 125 90 140 189 175 
Press. increase (mm) 10.5 20.5 4.8 1.4 0.3 
Cs3H¢ added 10.50 5.30 7.85 10.90 17.70 
CsHe left 9.25 4.00 7.10 10.55 17.25 
2 0.29 0.92 0.11 0.03 0.01 
CH, 0.78 1.51 0.32 0.12 0.14 
C3Hs 0.25 0.33 0.15 0.05 0.05 
CeHis (estim.) 0.5 0.5 0.3 0.15 0.2 
Ne 0.85 1.81 0.25 0.20 0.12 
NH: 0.50 0.56 0.23 0.20 0.16 
CH3NH2 0.60 0.70 0.30 0.2 0.19 
(CH3NH): used 1.9 2.9 0.8 0.5 0.5 








subjected to combustion. The carbon value 
found was 3. 

With large amounts of propylene, a higher 
stationary concentration of propyl radicals is to 
be expected and a lower hydrogen atom concen- 
tration. The table shows larger amounts of 
hexane relative to propane at higher propylene 
concentrations. Comparison with the results ob- 
tained in the absence of propylene shows that 
hydrogen production is very considerably re- 
duced and that the over-all decomposition is also 
retarded. 


DISCUSSION 


The absorption spectrum of ammonia! shows 
a region .of double-headed diffuse bands _be- 
ginning at 2400A and extending to 1665A. Below 
1665A the bands are sharp and show rotational 
fine structure. True continuous absorption begins 
around 1200A and below 1150A becomes so 
intense as to obscure the bands. The spectrum 
is attributed to an excitation of a non-bonding 
electron associated with the nitrogen atom. The 
spectrum indicates the existence of a predissocia- 
tion phenomenon over the range in which photo- 
chemical studies have been made, with the 
primary reaction postulated as 


NH;+hv—-NH;*—NH2+H. 


In the alkylamines absorption begins at 2300- 
2500A and extends in a series of diffuse bands to 
shorter wave-lengths. The intensity and diffuse- 
ness increases as the wave-lengths decrease ; with 
dimethyl and trimethylamine the bands even- 
tually give rise to a continuum. It appears from 
the similarity of the spectra of the amines and 


13 A. B. F. Duncan, Phys. Rev. 47, 822, 866 (1935) ; 49, 
211 (1936). 














THE 


ammonia that the cause of both is the same and 
Bamford” has postulated the primary rupture 


to be 
hv 


CH;NH.—-CH;NH + H. 


In hydrazine also, faint diffuse absorption 
bands, probably due to predissociation, begin at 
2370A. Bamford," from a study of hydrazine 
photolysis in presence of propylene and of nitric 
oxide, has furnished strong evidence that the 
primary rupture involves the scission of a hydro- 
gen atom from the hydrazine molecule. 

By analogy then, the continuum in dimethyl- 
hydrazine is probably ‘experimental’ and in- 
volves a predissociation. Light absorption results 
in the formation of an excited molecule which 
might conceivably undergo one or more of the 
following dissociations 


Ju +R (a) 
CH;NHNHCH;—CH;+R’ (b) (1) 
\2CH,NH. (c) 


Reaction (1a) would seem, by analogy, to be the 
most probable reaction and is supported by the 
reduction in hydrogen production when propylene 
is present, as also by the agreement between the 
quantum yields calculated directly and on the 
basis of H. Inspection of Fig. 2 shows that 
hydrogen production exceeds that of any other 
product. The presence of relatively large amounts 
of methane, as also of small amounts of ethane in 
the runs made with the hot arc for longer periods, 
indicates the probability of the presence of 
methyl radicals. The photolysis of methylamine 
does not give appreciable quantities of methane 
whence it is necessary to’ consider (1b) as a 
source of these radicals. Reaction (1c) is un- 
likely since it would be expected, initially at any 
rate, to yield larger amounts of methylamine 
than of the other products. A dissociation of an 
activated dimethylhydrazine molecule directly 
into stable molecules likewise seems improbable 
since larger amounts of methane relative to 
hydrogen than are observed would be expected. 

Little can be said definitely concerning the 
fate of the atoms or radicals produced in reac- 
tions (1a) and (1b). Hydrogen and methane 


“C. H. Bamford, Trans. Faraday Soc. 35, 568 (1939). 
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might result from such reactions as 


H+CH;NHNHCH;—H2+R, 


(2a) 


CH;+CH;NHNHCH;—CH,+R. = (2b) 


An alternative, though less likely reaction, for 
hydrogen production would be by recombination 
of hydrogen atoms. The relative weakness of the 
N-—N bond as compared with the N—C and 
N —H bonds suggests that the reaction of hydro- 
gen atoms and dimethylhydrazine might cause 
a rupture such as 


H+CH;NHNHCH;—CH3;NH2+CH;NH. 


Henkin and Taylor? present evidence that the 
energy of activation of this reaction is about 8 
kcal. The fate of the radical CH;NH might be 
pictured as 


2CH;NH—-CH;N : CH.+NHs. 


It is seen from Fig. 2 that while He, No, and 
CH, increase with increase in the percentage 
decomposition, NH; and CH3NHe decrease. This 
decrease is to be expected since each undergoes 
photolysis at the wave-lengths used and gives 
hydrogen and nitrogen as the chief products. 
The increase in methane with increased percent- 
age decomposition points to the accumulation in 
the system of some product which is a better 
potential source of methane than dimethylhydra- 
zine itself. In the absence of these secondary 
decompositions the curves should be horizontal 
lines. In any case they must start out, at small 
decompositions with a horizontal tangent. It is 
not unreasonable from the figure that all the 
curves might originate at a value around 0.3; 
that is, stoichiometrically three molecules of 
dimethylhydrazine are involved in the produc- 
tion of one molecule each of the measured 
products. A material balance between three 
molecules of dimethylhydrazine and one mole- 
cule each of Ho, Na, CH,, NHs, and CH;NH,. 
leaves a discrepancy of 4C, 10H, and 2N atoms. 
Methyl methylene imine CH;N :CH2 is known 
to polymerize to a dimer, its presence has been 
indicated among the products, and-the actual 
material balance on runs in Tables I and II 
supports this possibility. If reactions of types (1) 
and (2) are correct, the significance of the three 
molecules of dimethylhydrazine involved must 
mean that the radicals R and R’ produced in 
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reactions (1) and (2) require a third molecule of 
the hydrazine for ultimate decomposition giving 
over-all 


2R+CH;NHNHCH;=CH,+N2+NHs 


+CH;NH2+ (CH3N : CHa), 
or 


2R’+CH;NHNHCH;=H.2+ No+NHs; 
+CH;NHo+ (CH3N ° CHa)s. 


The observed quantum yield of 0.3 as opposed 
to the analytical yield of 3 is probably due to a 
reversal of (1a) or its equivalent. Henkin and 
Taylor? in the study of the addition of hydrogen 
atoms to azomethane found complete conversion 
of the azomethane to dimethylhydrazine even 
when the ratio of atomic hydrogen to azomethane 
was much less than 0.5 whereas, had the addition 
occurred entirely by steps involving atomic 
hydrogen, the ratio would have to have been at 
least two. It was obvious that some of the hydro- 
genation was occurring through hydrogen mole- 
cules and the reaction of the radical formed by 
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addition of one hydrogen atom to azomethane, 
with hydrogen : 


H.+CH;NHNCH;—~CH;NHNHCH;+H, 


appeared very plausible as a Walden inversion 
type reaction. This reaction is equivalent to a 
reversal of (1a) and, since the maximum energy 
of activation of addition of hydrogen to azo- 
methane was found to be 3-4 kcal., will occur 
quite readily reducing the quantum yield as 
found. 

The recent observation of Birse and Melville” 
that the efficiency of the primary process in 
ammonia photolysis is 0.58 rather than one indi- 
cates, since no fluorescence has been detected in 
ammonia and the pressure dependence of the 
quantum yield is less than would be expected for 
any collisional process, that there is some mecha- 
nism for internal degradation of the energy. This 
might also be present in the hydrazine molecule. 


16E, A. B. Birse and H. W. Melville, Proc. Roy. Soc. 
A175, 164 (1940). 
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The Diffraction of X-Rays by Liquid Oxygen* 


P. C. SHARRAH AND N. S. GINGRICH 
University of Missouri, Columbia, Missouri 
(Received May 17, 1942) 


The x-ray diffraction pattern of liquid oxygen at 89°K shows an intense peak at sin 6/A =0.157, 
and very weak peaks at 0.35 and at 0.5; that at 62°K shows an intense peak at 0.159, and 
weak peaks at 0.35 and at 0.5. The atomic distribution curves have peaks at 1.3A, 2.2A, 3.4A, 
and 4.2A for 89°K, and at 1.25A, 2.15A, 3.2A, and 4.1A for 62°K. The first peaks are due to 
the interatomic distance in the diatomic molecule, the second may be due to triatomic oxygen, 
and the others may be due to higher molecular aggregates. 


INTRODUCTION 


STUDY of the diffraction of x-rays by liquid 
oxygen is of interest not only for its contri- 
bution to the field of x-ray diffraction, but also 
for the information it may supply to our knowl- 
edge of molecular aggregates in liquid oxygen. 
The present work was undertaken to obtain as 
complete a diffraction pattern of liquid oxygen 
as possible and to analyze this pattern for the 


* A preliminary report of this work was presented at the 
December, 1940 meeting of the American Physical Society. 





atomic distribution curve. This has been done 
for liquid oxygen at 89°K and at 62°K. 

One of the earliest studies in the field of x-ray 
diffraction by liquids was done with liquid oxygen 
by Keesom and de Smedt. In their first work,! 
copper radiation filtered through nickel was used 
in an evacuated camera of 2.75-cm radius. Two 
peaks in the intensity curve were reported, the 
main peak at sin 6/A\=0.151 and a weak second- 


1W. H. Keesom and J. de Smedt, Proc. Amst. Akad. Sci. 
25, 118 (1922). 
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ary peak at sin 6/A=0.254. In their second work,? 
molybdenum radiation filtered through zircon- 
ium was used, enabling them to include larger 
values of sin 6/A. They reported, in this case, the 
main peak at sin @/A=0.154, a weak second peak 
at sin 0/A=0.239 on six out of nine films, and a 
somewhat questionable peak at sin @/A=0.35. 
Using Ehrenfest’s relation,? a=7.72\/47 sin 6, 
and the angle @ corresponding to the main peak, 
they calculated ‘‘a”’ and correlated this with the 
distance between molecules in a close-packed 
system. The primary result of this work was that 
it indicated that the intense peak was probably 
due to interference of x-rays diffracted by neigh- 
boring molecules. 

In view of the recent advances in both experi- 
mental technique and theoretical interpretation 
in this field, a repetition and extension of the 
early work of Keesom and de Smedt is clearly 
desirable. This attempt at an application of these 
newer techniques and interpretation has been 
made in the present work. 


EXPERIMENTAL 


Two types of system were used to obtain a 
sample of liquid oxygen, one which will be re- 
ferred to as an open system, and one which will 
be referred to as a closed system. In using the 
open system, a large amount of liquid oxygen was 
prepared from tank oxygen with no further puri- 
fication, and from pressure and temperature 
measurements, assuming it to have nitrogen as 
an impurity, its purity was found to be better 
than 96 percent. This liquid was transferred to a 
Dewar flask fitted at the bottom with a very 
thin-walled Pyrex glass tube 7.7-mm diameter in 
an evacuated camera of 9.53-cm radius, as shown 
in Fig. 1. As the oxygen evaporated, more was 
added, and for the low temperature work, a pump 
reduced the pressure over the oxygen. In using 
the closed system, which was nearly the same as 
that used for argon,‘ tank oxygen was given some 
further purification and introduced into the 
liquid air-cooled capillary at approximately one 
atmosphere. In this case, the reduced tempera- 


* W. H. Keesom and J. de Smedt, Proc. Amst. Akad. Sci. 
26, 112 (1923). 

°P. Ehrenfest, Proc. Amst. Akad. Sci. 17, 1184 (1915). 
(1940) Eisenstein and N. S. Gingrich, Phys. Rev. 58, 307 
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Fic. 1. A diagram of the evacuated camera using 
the ‘‘open system.” 


ture was attained by pumping on the liquid air 
in the Dewar flask. : 

Mo Ka radiation reflected from a _ rocksalt 
crystal was collimated in part by the crystal, and 
in part by two slits and the width of the beam at 
the sample was about 1.5 mm. A photographic 
film recorded the diffracted x-rays and exposures 
of 20 to 30 hours gave adequate densities on the 
films. A Leeds and Northrup recording micro- 
photometer was used to convert film densities to 
relative intensities, in the usual manner. 

The relative intensity curves were corrected 
for polarization’ and for absorption. The absorp- 
tion correction used was that for a cylindrical 
sample completely bathed in a parallel beam of 
x-rays, as discussed by Blake.* In the present 
case, the sample was only partially irradiated, 
and hence this correction is approximate. How- 
ever, it has been shown’ that for sodium, with 
cylindrical samples having diameters of from one 
to five times the width of the beam, there is little 
difference in the ratio of the intensity of scatter- 


5S. Katzoff, J. Chem. Phys. 2, 841 (1934). 
6 F. C. Blake, Rev. Mod. Phys. 5, 180 (1933), 
7G. A. Mitchell, unpublished work. 
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Fic. 2(a). (Left.) The x-ray diffraction patterns of liquid oxygen at 89°K and 
at 62°K. (b). (Right.) The atomic distribution curves for liquid oxygen at 89°K 


and at 62°K. 


ing at the main peak to that at angles as large as 
were used here. It is believed, therefore, that this 
approximation has introduced no serious error. 
The curve corrected in this way and for inco- 
herent radiation was fitted to the independent 
coherent scattering curve for making an analysis 
of the pattern. 

Seven exposures were made with the oxygen 
at 89°K, and six were made at 62°K, this lower 
temperature being attained by pumping on the 
liquid oxygen or on the liquid air. Some of the 
films were heavily over-exposed in order to dis- 
close features of the pattern at large angle where 
the intensity is normally very low. Typical curves 
are reproduced in Fig. 2(a). At each temperature, 
the pattern consists of one intense peak and two 
secondary prominences more in the nature of 
plateaus than of peaks. At 89°K, these occur at 
values of sin 6/X of 0.157, 0.35, and 0.5, while in 
the 62°K case, they occur at 0.159, 0.35, and 0.5. 
The position of the main peak at 89°K, 0.157, is 
slightly higher than the two values given by 
Keesom and de Smedt, 0.151! and 0.154,? but it 
is believed that with the larger camera here used, 
the present value is likely to be more nearly 
correct. Their second peak at 0.254 or 0.239 could 
not be found and it is believed that there is no 
peak near this position. Their third peak? at 0.35, 
however, does check well with a prominence 
found in the present work. 





ANALYSIS OF PATTERNS 


Fourier analyses of the patterns were made 
according to the method of Debye and Menke® 
giving the atomic distribution of atoms about 
any atom in the liquid. The expression used is 
derived elsewhere® and is 


4R*p(R) =47Répy +2R/x f Si(S) sin Rsds, 
0 


where i(s) is a function obtained from the experi- 
mental data, S=4z7 sin 6/X, 0 is half the scattering 
angle, p) is the average density of atoms, p(R) is 
the actual density of atoms, and R is the distance 
in angstroms from the origin atom. Four analyses 
were made for 89°K oxygen and five for 62°K 
oxygen. These various analyses were made for 
different patterns and for different fittings on a 
given pattern. In about half the cases, the inte- 
gral was evaluated by plotting Sz(S) sin Rs 
against S and integrating graphically, and for 
the remaining ones, the integral expression was 
evaluated by the method of trigonometric inter- 
polation.'® As is to be expected, the two methods 
gave very closely the same result for a given 
curve and fitting. The large number of analyses 
was made because of the unusual appearance of 


8 P. Debye and H. Menke, Physik. Zeits. 31, 797 (1930). 

9B. E. Warren, J. App. Phys. 8, 645 (1937). 

10C, Lanczos and G. C. Danielson, Phys. Rev. 55, 242 
(1939). 
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the distribution curve with many prominent 
peaks. An averaged curve for each temperature 
is shown in Fig. 2(b). In all the curves obtained 
from these analyses, most of the essential features 
shown in the figure were reproducible up to a 
little over 5A, but beyond that distance the vari- 
ous analyses differed considerably, and hence, 
no confidence is placed in the curves beyond 5A. 
At less than this distance, the first four peaks 
appeared in all cases at approximately the same 
positions but with somewhat different heights. 
This feature appeared to be unusually dependent 
upon the fitting of the corrected experimental 
curve to the independent coherent curve for the 
peaks between 3A and 5A. The first peak was 
reproducible to within about 5 percent of the 
reported value for position and to within about 
10 percent of the reported value for the area and 
the weak second peak persisted under all fittings 
with variations in position and area slightly 
greater than for the first peak. 


DISCUSSION OF RESULTS 


The positions of the peaks for the 89°K case 
are 1.3A, 2.2A, 3.4A, and 4.2A, and for the 62°K 
case they are 1.25A, 2.15A, 3.2A, and 4.1A The 
area under the first peak, for 89°K, is 1.08, and 
that for 62°K is 1.18, indicating that there is one 
nearest neighbor, or that liquid oxygen has a 
fundamental diatomic aggregate closely similar 
to the diatomic molecule in the gas, for which 
each atom has one nearest neighbor at about 
1.22A."." From this analysis, one cannot be sure 
that there is a significant difference between 89°K 
and 62°K, in the positions of this first peak, i.e., 
between 1.3A and 1.25A. The area under these 
peaks is slightly more than one, but it is believed 
that the precision and reproducibility of this 
peak is hardly great enough to give this difference 
much significance. Liquid oxygen is the only 
element thus far to show a second discrete peak, 
which occurs here at 2.2A and at 2.15A in the 
cases of 89°K and 62°K, respectively, and the 
interpretation of this second peak, if it is real, 





" A. H. Compton and S. K. Allison, X-Rays in Theory and 
Experiment (D. Van Nostrand Company, 1935), p. 161. 

® G. Herzberg, Molecular Spectra and Molecular Structure 
(Prentice Hall, 1939), p. 491. 


is somewhat difficult. It is true that from the 
assigned structure of ozone,'* O3, the first inter- 
atomic distance is between 1.25A and 1.30A and 
that the second distance is about 2.2A, but the 
amount of QO; present in liquid oxygen is normally 
very small. On the other hand, irradiation of 
liquid oxygen is known" to enhance the amount 
of O; present, which would contribute both to 
the presence of this peak and to an area under 
the first peak greater than one, since one atom 
of O; would have two neighbors at about 1.3A. 
Thus the distribution curves here obtained are 
consistent with the assumption that there was 
O; in the samples used, though the data are not 
considered sufficiently well known to warrant 
making quantitative deductions from this result. 

The third and fourth peaks appear to be re- 
lated to distances between oxygen atoms in 
larger aggregates of oxygen molecules. From 
studies of the magnetic susceptibility and of the 
absorption spectra of oxygen, it appears prob- 
able that O, may exist in the liquid in a con- 
siderable amount. Vegard!® has determined the 
crystal structure of the y-form of crystalline 
oxygen, and he states that the x-ray diffraction 
pattern was typical of that for a cubic lattice with 
rotating molecules. The O2 molecules are grouped 
into pairs along the trigonal axis and the distance 
between the two molecules of a pair is given as 
3.48A. In the present work, two distances, about 
3.4A and 4.2A, appear fairly predominantly, and 
it is not unlikely that the former may have some 
relation to Vegard’s value, 3.48A, though an 
unambiguous interpretation of both these peaks 
appears somewhat difficult. 


ACKNOWLEDGMENTS 


It is a pleasure to acknowledge financial sup- 
port of this work by the American Academy of 
Arts and Science, the American Association for 
the Advancement of Science, and the Elizabeth 
Thompson Science Fund. 


13 W. G. Penney and G. B. B. M. Sutherland, Proc. Roy. 
Soc. A156, 672 (1936). 

14 J. W. Mellor, A Comprehensive Treatise on Inorganic 
and Theoretical Chemistry (Longmans, Green and Com- 
pany, 1922), p. 880. 

15 M. and B. Ruheman, Low Temperature Physics (Cam- 
bridge University Press, 1937), p. 228. 

16. Vegard, Nature 136, 720 (1935). 














AUGUST, 1942 


JOURNAL OF CHEMICAL PHYSICS 








On the Absorption Spectrum of Some Polymethine Dyes 
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(Received April 29, 1942) 


The energy levels of certain polymethine dyes are calculated approximately. The dyes consist 
of a chain of conjugated double bonds, with end groups which might or might not differ from 
each other. Ions and neutral molecules are treated. The calculations are performed both with 
the valence bond method (resonance between different structures) and the LCAO method, 
but without antisymmetrization and without taking electronic interaction into account. The 
effect of chain length and of the nature of the end groups is discussed. 





I. INTRODUCTION 


N a series of papers, Brooker and Sklar! discuss 
the color, or in general, the first electronic 
transition in a series of dyes, mainly of cyanines. 
The present paper has as its purpose an ap- 
proximate quantum theoretical discussion, which, 
while insufficient to permit an actual calculation 
of the wave-length of the absorption, would at 
least permit us to understand the dependence on 
the chain length and on the nature of the end 
members of the molecules. 

As is well known, there are essentially two 
methods to treat the electron states, the Heit- 
ler-London-Pauling-Slater method (abbreviated 
HLPS) and the molecular orbital method 
(Hiickel-Hund-Mulliken, abbreviated LCAO). 
Beforehand, we cannot know which will fit 
better. It turns out, however, that the first part 
of the formal calculation is the same in both; only 
the numerical value of the constants involved 
and the final application of the formula are 
different. 

Let us, as an example of the structures dis- 
cussed in the HLPS method, consider the ion 


& 
H:.N=C—(C=C),—NH. © (1) 
H HH 


It contains N=2n+1 CH groups. According to 
the HLPS viewpoint, there are two ground struc- 
tures, namely besides (1) 


H,N—(C=C),—C—=NHp, (1’) 
HH H 


the (negative) energy of which we call e. Besides 
the two ground structures there are N+1 higher 


1 For the literature, see the following paper by A. L. Sklar. 
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structures 


H,N—C—(C=C)—NHp, (2) 
H HH 


H.N—C—C—C€—(C=C),_.—NH». (3) 
H HH HH 


Formulas (2) and (3) do not have exactly the 
same energy, but in first approximation we ne- 
glect this difference and ascribe to the excited 
structures the energy f, which is smaller than e 
(i.e., the excited structures lie above the ground 
structures). Any structures lying still higher are 
neglected. This would not be justified for a nu- 
merical calculation, but, we hope, does not matter 
for the qualitative discussion attempted here. 

Furthermore, it is assumed that only neigh- 
boring structures interact, (1) with (2), (2) with 
(3)---(N+1) with N+2. The interaction energy 
is called a and is assumed to be the same between 
any two states. 

When the energy levels—of which there are 
N+2—are calculated, the HLPS method takes 
the lowest one as the ground state and the next 
higher one as the first excited state. The first ab- 
sorption band is therefore given by the energy 
difference of the two lowest states. On the other 
hand, the molecular orbital method starts from 
the skeleton 


+ + + + + 
H.N—C—C ... C—NH, (4) 
H H H 


into which N+1 electrons are put. There will be 
(N+2) states, of which (N+1)/2 are filled. The 
first transition takes an electron from state 
(N+1)/2 and puts it into (V+3)/2. However, 
the constants have entirely different meaning and 
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value in the HLPS and the LCAO method, even 
if they appear formally in the same manner in 
the equations. 


The LCAO method consists in forming a 
molecular orbital of the whole molecule out of a 
sum of atomic orbitals for the different ions, 


+ + 
C or NH. Mathematically, this task is identical 
H 

with that of the HLPS method. 


II. SYMMETRICAL DYNES, N>1 


For dyes like the ion (1), which have equal 
ends, the energy is determined by the determi- 
nant of order N+2 


e—E a 0 0 0 0 
a f-E a 0 0 0 
0 a f-E a 0 0 
0 0 =0. (5) 
a 0 
a f-E a 
0 a e—E 








The problem of the dienes e=f has been ex- 
tensively treated.” There, however, the difference 
between the end members and the chain members 
is neglected, while the influence which this differ- 
ence exerts is one of our principal interests. On 
the other hand, Lennard-Jones, Turkevich, and 
Coulson have taken into account the fact that 
the chain double bonds might alternatively have 
different lengths, a fact we neglect. 

To solve (5) we write the classical differ- 
ential equations of a system of coupled oscillators 
for which (5) is the secular equation with w? in- 
stead of E. Calling the vibration coordinate of 
the end members X and Y and that of the sth 
chain members x,, we have * 


X=-eX—on, (6) 
¥)= —fx1—aX —axo, (7) 
€,= —fx,—a(Xs1+%5-1), S=2,---N-1 (8) 
&n= —fxn—aV—axn-1, (9) 
V=—eY—axp. (10) 


*E. Hiickel, Zeits. f. Physik 70, 204 (1931); 76, 628 
(1932); Zeits. f. Electrochemie 43, 752 (1937). R. S. 
Mulliken, J. Chem. Phys. 7, 364, 570 (1939). J. E. Lennard- 
Jones, Proc. Roy. Soc. 158, 280 (1937). J. E. Lennard- 
Jones and J. Turkevich, Proc. Roy. Soc. 158, 297 (1937). 


—_ Proc. Roy. Soc. 164, 383 (1938); 169, 413 


This can be interpreted as the equation of motion 
of a vibrating chain in which the end members 
are tied to positions of equilibrium in space 
by quasi-elastic forces —(e—a)X or (e—a)Y, 
the chain members with quasi-elastic forces 
— (f—2a)x,. In addition there act between neigh- 
boring particles destabilizing quasi-elastic forces 
+a(x,—Xs1), a(X —x), ete. 

It turns out that the system of differential 
Eqs. (6)—(10) is completely equivalent to one 
treated by Bartholomé and Teller.’ These authors 
treat the vibrations of a carbon chain (quasi- 
elastic force constant k,, mass m,) with valence 
angle 8, to which end members (quasi-elastic 
force constant ke, mass mtg) are attached. (6)- 
(10) go over into their equations if one identifies 
as follows: 





1 1 2k. k. 
E:w?; e:ke —+—),; f:—; a=—cos8. 

Mp mM, mM. mM. 
Furthermore, nowhere do they make use of the 
fact that cos B=1, so that no limitation is put on 
the ratio of a/f. Bartholomé and Teller have not 
actually written the equations for the energy 
levels or discussed them. That shall be done here. 
One finds as a general result, which can easily 
be confirmed from the following, that there are 
N+2 solutions, that (V+3)/2 are symmetric 
and (N+1)/2 antisymmetric, that the sym- 
metric and antisymmetric levels alternate, with 
a symmetric one at the bottom and at the top. A 

symmetric solution is one for which 


X=Y, xqw4n/2}-2=Xw41)/2] 48 (11) 


N-1 
et ee 
2 
an antisymmetric one has a minus sign on the 
left of both Eqs. (11). It follows that for all cases 
in the HLPS method and for an odd number of 
double bonds in the LCAO method, the first 
transition takes an electron from a symmetric 
state and puts it into an antisymmetric one.‘ 
We introduce the quantity 


R=a/(e-f), | (12) 


3 E. Bartholomé and E. Teller, Zeits. f. physik. Chemie 
B19, 3 (1931). 

‘For an even number of double bonds, the highest filled 
one-electron state is odd, and the first excited state even. 
For the state of the ion as a whole, however, the statement 
in the text is true whatever the number of double bonds. 
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which gives the ratio of the interaction energy to 
the energy difference between the ground struc- 
tures and the excited structures. Depending on 
the question whether R>1, =1, <1, the discus- 
sion will be different. 

To solve (6)—(10), one writes (A ~e*‘) for 
symmetric states 


N+1 
x,=A cos (——-s)e, (13) 
2 
for antisymmetric ones 
N+1 
x,=A sin (——-s)e. (14) 


k is the ‘“wave vector” of the standing wave, and 
corresponds to 7/\ in Bartholomé and Teller’s 
paper. (In addition, they count s from the middle 
instead of from one end.) 
Both (13) and (14) inserted in (8) give 
E-—f=2a cos k. (15) 
X is then found, from (6) as 


X=xy— (16) 
=X _ 6 
'E —e 
One gets, by inserting (15), (16), and (13) or 
(14) in (7), for the symmetric states 
N+1 
k 


2 
——————— a R, 17 
was (17) 


2 








cos 





cos 


for the antisymmetric ones 
N+1 
k 
2 
N+3 
k 


sin 
2 





sin 


=R. (18) 





If R21, these equations give the appropriate 
number of solutions, but if R<1, the number is 
too small by two. The additional solutions are 
found by putting® 


Xe=A {ylN+D/2—-s4y-(N+H/2+8} — (13’) 


5 This can be expressed, by writing v=e~*, as 
N+1 


xs=A cosh (AF -s)e, xs=A sinh (“F-s)e. 


(15’) is then E—f=2a cosh k. 
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where the upper sign is for the symmetric, the 
lower for the antisymmetric solution. (15) is then 
replaced by 
1 
B-f=a(v+-). (15’) 
(15’) has, for a given E, two solutions v which are 
reciprocal. We decide to call v the one which 


is <1. 
Then, (17) and (18) are replaced by 


(17’) 





One can easily prove again that the symmetric 
solution is the lower one and that both solutions 
of (17’) are below® those of (17) and (18). 

However, for R near, but smaller than 1, (17’) 
does not give an antisymmetric solution. In fact 
one finds, using the negative sign and putting 
v=1 (which is the largest value v can have), 

2 


=i —-—_. 


For R smaller than this value, (17’) provides both 
a symmetric and an antisymmetric solution. For 


2 
jo LE <i, 
N+3 


(17’) provides the symmetric solution, while an 
additional antisymmetric solution comes from 
(18), with 
N+1 
——k<.-. 
2 2 


Th. Foerster,’ in a very interesting paper, has 
given the formula (13), (14) (15), (17), (18). He 
evaluates them for R=0, 1, «, which can be 
done exactly, and discusses the trend of the ab- 
sorption frequencies from his results. He does 
not attempt to give approximate expressions. 

The rest of the discussion will be made for 
R<1i, R=1, R*1. 

1. R<1 or R<1, N>1 or both conditions. 

In the HLPS method, we are interested only 
in the lowest two solutions, i.e., the two given by 
6 y+1/y is larger than 2, except for v=1; 2cosk is 


smaller than 2, except for k=0. 
7 Th. Foerster, Zeits. f. physik. Chemie B47, 245 (1940). 











’ ee om — _ 


rr -— 70 <7 * 





Pe 


) 
ct 
1g 


th 


or 


an 
m 


ily 
by 


+0). 











POLYMETHINE DYES 511 


(17’). Then, in zero approximation, the right side 
is zero and 
v 
1——=0. 
R 


For the next approximation, this is inserted in 
the right and v‘*! in the denominator neglected. 
Then® 


9 


1-—=+(1-—R?)R®*', 
R 
v= RF (1—R?)R**?, (19) 
1 a 
E-f=a(—-+R) +—(1—R?*)?R‘t! 
R R 
and finally 
a? 
eo (20) 
eé=— 


so that we see: the center of the two lowest levels 
is depressed below e by the amount a?/(e—f), 
which is independent of N, while the splitting 
depends on N. The energy difference between 
these two levels, which in the HLPS view gives 
the first transition, is 


AE=2a(1— R*)?R%. (21) 


(21) says that the wave-length of the absorption 
band goes to infinity exponentially with N. 

Our results have a close connection with the 
assumption of Pauling® that the large number of 
intermediate structures can be replaced by an 
increased overlap. That is so according to our 
calculations, but the quantitative side needs 
some discussion. In the present case, the de- 
pendence on N is of the same form which it would 
have in the case of an increased overlap integral. 
The latter usually has the form C exp (—7r/p) 
where C and p are constants and ¢ the distance 
between the overlapping atomic orbitals. If we 
now write for the chain r=r,+Nre, the overlap 
integral would be (C exp —r:/p)(exp —7r2/p)% 

* This is a good approximation provided R¥+! is small 
compared to one. E.g., if R=} and N=3 (a single C=C 
group), RY+1= 4 and the zero approximation v= R is only 
wrong by 5 percent. The split, of course, is wrong by a 
greater error (about 20 percent). For N=5, the error in 
the split is less than 7 percent. It would be easy to improve 
the calculations for small N but not worthwhile, because of 


the simplified assumptions. 
*L. Pauling, Proc. Nat. Acad. Sci. 25, 577 (1939). 


which is of the form (21) with R instead of 
(exp —?2/p). However, Pauling is not interested 
in the dependence on chain lengths, but in the 
influence of different chain endings, and therefore 
he expects the variation with e to be proportional 
to (e—f), i.e., expects a transition frequency 
~R for fixed N and a. None of the formulas dis- 
cussed in this chapter behaves quantitatively in 
this manner; only those found with the assump- 
tions of Chapter V do. 

In the HHM view, the lowest transition is the 
one which takes an electron from the topmost 
full band, the (N+1)/2th one, and puts it into 
the lowest empty one, i.e., the (V+3)/2th one. 

We have therefore to investigate the solutions 
of (17) and (18). If Ris very small, the numera- 
tors of (17) and (18) will be very small; therefore 


N+1 7 
en (r=1---N) (22) 


r odd for symmetric solutions (i.e., (17)), even 
for antisymmetric ones. 
Then (17) and (18) can be written in the same 


form 
N+3 rr 
sin 6=R sin ( 6+ ). (23) 
N+1 N+1 





For the highest filled state, r=(N—1)/2 and 
E=f+2acos { —— -—-+——— 
N+12 N+1 
w—26 
=f+2asin ———, (24) 
N+1 


N-1i-r ~~) 





N-1i2 N+3 
sin i= sin ( ) 


> 61 
N+12 N+1 


T N+3 
=+R cos ( = i). (24’) 
N+1 N+1 





For the first empty state, r=(N+1)/2 and 








E=f-2 262 
=f—2a sin = (25 
N+1 
N+3 
sin 6.=R cos "hi (25’) 


Up to now, the formulae are exact. 
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If R is very small and WN arbitrary, (24’) and 
(25’) give approximately 








6,:= R cos ——, 5.=R, 
T 2R 
AE= 2a} sin ~ 
N+1 N+1 
T 2R | 
Xcos* ——-+———_ 
N+1 N41) 
T 2R ra 
= 2a sin + sin? ——-._ (26) 
N+1 N+1 N+1 


If N is very large and R has any value <1, 
(24’) and (25’) give approximately 


TT TT 
tg 6:=R cos (_)(14R sin —_) 
N+1 N+1 


tg bo=R. 


The energy difference is, in first approximation, 


T 62— 461 
(142 ) 
+1 T 
or, if R is small 


rg R 
AE= 2a—"_( 1 +2), 
N+1 N+1 


In both (26) and (27), the wave-length goes to 
infinity”’ with increasing N. The principal term, 
which is independent of the nature of the ends, 
goes with the first negative power of N+1. 

2. For R=1, (17) and (18) can be solved 
exactly, by 








AE=2a 
N 


(27) 





N+1 N+3 
—_—h=rzy———k 
or (28)" 
as 
k= , r=0,1---N+1 
N+2 





10The expression for AE is really a power series in 
(N+1)~". If the wave-length is calculated, the first term 
is proportional to N+1, the next term is independent of N. 
The first two terms of the expressions for the wave-length, 
taken together, are therefore proportional to N+C, where 
C is a constant independent of N. The following term is 
~(N+1)7. 

1 y= N+2 is odd and k=z is not a solution of (18). 
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where r even is valid for the symmetric, r odd 
valid for the antisymmetric states. 
The two lowest states are 


E=f+2a=e+a (symmetric) (29) 


and 





T 
E=f+2a cos 
N+2 


T 
+0(2 cos —-1) (antisymmetric). (29’) 
N+2 


The energy difference for the first transition ac- 
cording to HLSP is 


T T 
) = 4a sin®?————-._ (30) 
N+2 2(N+2) 


For very large N, the wave-length of the absorp- 
tion band goes to infinity like N? [(more pre- 
cisely, it is hc/4asin~*2/2(N+2)]. For the 
LCAO method the lowest AE is given by the 
difference between the (V+1)/2th level and the 
(N+3)/2th one, i.e., between the levels with 
r=(N—1)/2 and (N+1)/2; therefore 


(N+1)x —— 
a ae a 
2(N+2) 2(N+2) 





AE= dal 1 — cos 


AE=-—- 2a cos 


TT TT 
= 4a cos ———— si 


(N+2) 


The wave-length of the absorption band goes to 
infinity with increasing N, but with the first 
power. 

3. For R>1, (17) and (18) provide all the 
energy levels, i.e., (13) and (14) all the solutions. 
It is now better to take the reciprocals of (17) 


and (18). 














N+3 
cos k 
1 
—_—_——_—_— =— (symmetric), (17”’) 
N+1 R 
cos k 
2 
N+3 
sin k 
1 
—_——_—=— (antisymmetric).  (18’) 
N+1 
sin k 











Id 


9) 


31) 


, to 
irst 


the 


ns. 


17) 


18’) 








We expect now 


N+3 T 
——k=r-+6 (32) 
2 2 


with odd ¢ for the symmetric, even ¢ for the anti- 
symmetric states. 

Then (17’’) and (18’) are again expressed in the 
single formula 


1 rr N+1 
sin j= —— sin ( _ s). (33) 
R N+3 N+3 


r=0 and r=N-+3 give no solution 6<2/2 which 
(32) presupposes. 7 goes therefore from 1 to N+2, 
the two lowest states being r=1 and r=2 
(HLPS) while for HHM the highest occupied 
state is r= (N+1)/2, the lowest unoccupied state 
r=(N+3)/2. We approximate again for R>1. 
(33) is then 

















1 rr N+1 1 rr 
5= —— sin + — sin (34) 
R N+3 N+3R N+3 
and the energy 
Tr 26 
E=f+2acos | + 
N+3 N+3 
rT 2 1T 
f-+2a cos". 1 sin’ 
N+3 R?(N+3)? N+3 
4a Tr 





+ sin 
R(N+3) N+3 








Tr N+1 1 YT 
i | + — sin ; 
N+3 N+3R N+3 
Or 











‘Tr e—f ™r 
E=f+2a cos +4 (sin ) 
N+3 N+3 N+3 


1 N Tr 
x{14+— cos | (35) 
R N+3 N+1 








The first two terms would appear in the pure 
polyene chain of N+2 members. We have then, 
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for HLPS 


3 of T 
AE =4a sin — sin 
2N+3 2(N+3) 











¢= T TT 
—4 sin? —~=\4 cos? - 
N+3 N+3 N+3 


1 WN 7 


-14+— ——_ coos ——_ 


R N+3 N+3 


T 2 
x ( 400s- cos —.~1)}. (36) 
N+3 N+3 





Now the wave-lengths of the absorption band go 
to infinity with (V+3)? with a correction propor- 
tional to N+3 (just as for R=1). 











For LCAO 
fa 
XE=2e( sin. ) 
N+3 
N 2 T 
{14+— - — cos? —— 
(N+3)? R? N+3 
4(e—f) cs 
- f sin? . (37) 
N+3 N+3 


There we have again, for large N, the same be- 
havior as for R=1. 

4. For values of R and N, for which our ap- 
proximations are not good, numerical calcula- 
tions have been performed. E, and therefore k or 
v, were assumed and the corresponding R calcu- 
lated from (17), (18), (17’). The results were 
plotted and selected values read, which are given 
in Table I. 

5. Change in interaction energy of the ends. 
We can also treat the case in which the ends of 
the ion have a different interaction energy, a’, 
with the C atoms than the C atoms among them- 
selves. Then (8) remains unchanged, while (6) 
and (7) are replaced by 


(E—e)X =+a'x, (6’) 
(E—f)x=a'X+ax2 (7’) 


with corresponding changes in (9) and (10). For 
the lowest states, (14’) and (15’) remain un- 
changed, while the energy equation (replacing 











(17’)) takes the form 

















v 1-—v? 
1 —— = yv? tvNt1— (17’”) 
({ty)** 
with 
a? 
y=1-—. 
a’? 
Similarly (17) and (18) are replaced by 
N-1 
cos k 
1 a”? 2 
2 cos k-—-— = — ———"—_- (18’’) 
R & N+1 
cos ; k 


and a corresponding equation with sinus, (13), 
(14), (15) remaining unchanged. The calculation 
shall be limited to the two lowest states, given 
by (17”) for R<1. It is better to divide (17’’) by 
v. One then finds 


1 1 
~= [1+ (1-49?) c01(1 2) 





v (1+4yR2)) 
Call!” 
(1+4yR?)?—1 
=é€. 
27R 
Then 
. a 
gg eneng hI] gl) 
ve (1+4yR?)! 
and 
= ¢-+————¢"#1(1 —e?). 
(1+4yR?)! 


Therefore 


1 Ra 
E=j+a(—+e) +——_—_—_———¢"-1(1—*)? (20’) 
€ (1+4yR?)! 
and the energy difference of the first transition 
according to HLPS 


2Ra 


E=—————(1 —e”) *e4!" (21’) 
(1+4yR?)? 


(21’) has the same form as (21), only the 
constants are different, ¢« replacing R and 
2Ra/e(1+47R?)! replacing 2a. The coefficients of 
the different structures in each state, according 
to HLPS, are given by (13), (14), (13’), and (16). 
In the LCAO interpretation these are the co- 


12 It is assumed that | yR| <1, but y may be negative. 
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TABLE I. AE/a, accordingly to HLSP. 











R N=5 N=7 N=9 
0.70 0.105 0.021 
0.75 0.120 0.058 0.032 
0.80 0.135 0.072 0.044 
0.85 0.17 0.085 0.052; 
0.90 0.18 0.100 0.061 
0.95 0.19 0.105 0.074 
1.00 0.20 0.120 0.081 
1.5 0.28 0.185 0.135 
2.0 0.325 0.220 0.160 
2.9 0.345 0.235 0.170 

10 0.410 0.270 0.195 








efficients of the atomic orbitals in the molecular 
orbital. In both cases, these coefficients are not 
yet normalized. 

The transition moments have been calculated 
as ~N+2 by Mulliken" (for polyenes). 


Ill THE ENDS ARE UNEQUAL 


In this case we again assume the same inter- 
action energy between X and x1, as between Y 
and x,, but take the energies of X and Y to be 
e; and é2 instead of e. (10) has to be changed ac- 
cordingly. The states are now not anymore 
purely symmetric or antisymmetric. 

Therefore we try to integrate (8) by 


X= avl N+ /21—s 4 hy-(N+D/2) +8 (38) 


which takes the place of (13) and (14) if v=e* 

and the place of (13’) if v is real. The difference 

between a and b expresses the fact that the two 

progressive waves which make up the standing 

wave (13) and (14) are not any more equal here. 
Insertion of (38) into (8) leads again to (15’). 
Then, (6) and (7) give 

















a? Xe a? 
E-e, X41 E-e, 
apl\N+1)/21-2.4 hy UN +1)/2)42 
a . . (39) 
aviN+1)/2)—-1 4 hy lN+1)/2} +1 
Similarly, (9) and (10) result in 
a? xn-1 89a? 
E-f= +a = 
E-eés XN E-e. 
| apl\N+0/21-(N—D) 4 by- UN +0214 (N-1) 
+a . (40) 


avl(N+D/2)-N 4 hy UN+D/2)4N 


~ 8R, §, Mulliken, J. Chem. Phys. 7, 20, 364 (1939). 








v 








ar 
ot 


ed 


i 


be 
Ac- 
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38) 


1ce 
wo 
ing 


cr 


5’). 


(40) 








Eliminate b/a between these two equations, and 
the following equation is found from a rather 
lengthy calculation: 

(E-f a 1 
vy-%) oineascale | 


a E-e, v 








E-f «a | 





a E-e. v 





“Pr. 


E-f a E-f a 
+ -»}| + 


a E-e, 








a E-e, 


Using again (15’) one finds 


U 


R 
y?V+2(1—vR)(1—vR’) = | 1 — | 1 -—}, (41) 


v v 
where 
a a 
R= , R=——. 
ei—f éo—f 
For R= R’, this reduces to (19), as one can easily 
show. We call the lower structure e)(e;> és, 
R<R’). 
1. If R, R’<1, we find, as in II, two solutions 
with real v. In this case zero approximation gives 
(by setting the left side zero) 


v=R, v= R’. (42) 


This inserted in the left side of (41), leads to 





v= 3(R+R’)—3{(R’—R)? . 
+4R2%+4(1—R?)(1—RR’)}3, (43) 
v' =3(R+R’)+3{(R’—R)? 
+4R’2N+4(1—R”)(1—RR’)}4. (43’) 


One can with the same approximation, write 1/v 
and 1/v’. Abbreviate 


o?=4R*'(1— R?)(1— RR’) (44) 


and call o’ the result of exchanging R’ and Rin o. 
One finds then for the lowest levels 


ei +e2 


2 





Qa 
By=——+ (RR) +41 (01 e2)* 0%"! 


R?)! 
—+3RR’ alta P (45) 


_ ‘The numerical limits within which the approximation 
is good are the same as given in footnote 8. 
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eites a 
2= FUR+R)— Fi (eres)? 
tate!) 43RR| (€;—€2)?’ + aa"? : | . (45’) 
Therefore, the split is 
AE=}3{(e1—e2)?+a7o"}} 
+21 (€1—e2)? +a%o’?}} 
R?}3 
—4RR'| (e.-e1)*+ate | 
R” 
12,3 
—4RR’| (e.~es)*+a%0" | (46) 


In (43), (43’), (45), (45’), (46) the absolute value 
of the square root is meant. If N is sufficiently 
big so that (e;—é2)?>a’o?(R?/R”)... (which is 
more stringent than the condition for the ap- 
plicability of (43), if e:—e2 is small compared to 
e—f), the formula can be simplified to 


AE= (e;—e2)(1—RR’) 


a? 1 R? R® 
+— |@"(1-—) (1-—)}, (46’) 
4 €\;—€2 R’ R 


which for N=, approaches the finite value 
(€:—e2)(1— RR’). If R, R’ are larger than 1, and 
for the higher levels of the case where R, R’ are 
smaller than 1, one has to substitute for v in (38) 
e* and gets then 
RR’ sin (N+3)k—(R+R’) sin (N+2)k 
+sin (N+1)k=0. (47) 
k=0 and k=z7z, which formally satisfy (47), 
are not solutions of the problem.'® This also is 
true for the corresponding value v= 1 in Eq. (38). 


(47) takes the place of the two Eqs. (17), (18). 
We now solve again separately for the two cases: 


Rei, <i and R>1, R’>1. 





4 From these values of k (or v=1) follows that x, would 
be independent of s. Then Eqs. (6), (7) give - 


ae 


(E-f)= ste 
while (9), (10) lead to 





2 
E-f=;—,,ta, 


which cannot both be true for e; e.. 








In the former, we would have, for R=R’=0, 
sin (V+1)k=0, or k=rx/(N+1). We therefore 


set 








rT 
k= +5, r=1---N (48) 
N+1 
and 
2rr §=609wNN+3 
RR’ sin | a i) 
N+1 N+1 


rT N+2 
—(R+R’) sin | oe 
N+1 N+1 





5} sin 5=0. (49) 


If N is very large one can write, to the first 
approximation 





tT 
tg b= | R+R'-2RR’ cos 


N+1 
rT 
x (sin 
N+1 





rT 
)} 1—(R+R’) cos —— 
N+1 


2rr 


+ RR’ cos 
N+1 





q . (50) 


The highest filled state, in the LCAO method, 


is again 





or 





T 
tg i= | R+R’—2RR’ sin 
N+1 








Tv TT 
x (cos )) 1—(R+R’) sin 
N+1/ | N+1 


— RR’ cos 





2r 7 
or 


T 
tg is~ | R+R’—2RR'—_| 
N+1 


x {1-RR'(R+R)"_} (51) 
N+1 


or even more roughly 


R+R’ 





tg 6,= (51’) 


1—RR’ 





516 KARL F. HERZFELD 


For the first empty state r=(N+1)/2 
tg 62=(R+R’)(1— RR’), (S1”’) 


which, to this approximation, agrees with tg 4;. 
The energy of the transition is then 





N-1i7 
AE=2a cos| “+04 
N+12 





N+17-7 
~2a cos “+61 
N+12 


= 2 (1 RR')\(14 R914 RI. (52) 
N+1 


This formula is correct for sufficiently large N, 
without demanding that R and R’ have to be 
very small compared to 1. 

Therefore, for LCAO we find here, just as for 
equal ends, that the wave-length of the first 
transition goes to infinity as N+1. 

2. Instead of the special case of R=1 of II, we 
have here the special case RR’=1 which can be 
solved exactly. Then (47) takes the form 


1 
(sin (v4+2)8)|2 cos k— (r+) |-0 (53) 


with the solutions: 
k=rnr/(N+2) 


and, in addition, the lowest lying term, given by 
v=R (Eq. 42); v= R’=R™ describes the same 
state. Therefore the two lowest states are 


r=1---N+1 (53’) 


1 
Ey=f+a( R+—) =f+(erter—2f) (54) 


and 





E,=f+2a cos i . (54’) 
N+2 


The first transition in the HLSP approximation is 


TT 
(55) 
N+2 


and approaches, for N= ~, the finite value 





1 
AE=a| R+—— 2 cos 
R 


1 
SE. =ol R+——2| ete, 2f—20 


=([(ei—f)'—(e2—-f)'}. 55’) 








S 


Si 








2) 


be 


or 
rst 


we 


be 


53) 


3’) 


me 








For a LCAO description, the highest filled 
state is given by r=(N—1)/2, the lowest empty 
one by r=(N+1)/2. These have the energy 


N-17 


T 
E,=f+2a cos —=f+2asin , (56) 
N+12 N+1 








N+i-r 


E.=f+2a cos —— —=f (56’) 
N+12 





and the energy of the first transition is 


E=2a sin 





(57) 


with a wave-length going to infinity proportional 
to N+1, as for R=1 in a symmetrical ion. 
3. For R>1, R’>1 we rewrite (47) as 


1 1 
sin w+3k-(+—) sin (V+2)k 
gs 


1 
+— sin (N+1)k=0 (47’) 
RR’ 


and solve by putting 











rr 
k= +6, r=1---N+2 (58) 
N+3 
or 
1 1 (V+2 Yr | 
sin o-|-+—] gia § -——-§ ---—— 
R RJ \|n+3 N43!) 
1 N+1 2rr 
+ sin | ——§— | =o. (59) 
RR’ N+3- N+3 


For large N, one finds now 

: 1 1 2 YT 

ind=-| (—4+— - “COS 
R R‘/ RR’ N+3 


TK 1 1 Tr 
x (sin \P-(G+5 cos - 
N+3 = x N+3 


1 
+—— cos 
RR’ 

















2rmr ]7! 
. (59’) 
+3 


Therefore, for ground and excited state in 
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DYES 
HLPS (r=1, 2), 
1 1 2 
wolfe 
gk KF BF 
.- 3 '; T 
x | i-—_-_+—_| -—_, @ 
R R’ RR’) N43 
62= — 26, (60’) 
TT 
AE=2a cos] +0] 
+3 
2r 
~2a cos| +251] 
N+3 








7 2 : & 2 
en iecuae 
N+3 R R RR 


14 #177) 
x|1-—-— + z (61) 
R R' RR’ 





Therefore the wave-length goes to infinity like 
(N+3)? as for equal ends. 

If the LCAO gives the better description, then 
we are interested in the transition between the 
terms with r=(N+1)/2 and r=(N+3)/2. 

The corresponding values of 6 are, in the zero 
approximation, for both levels 


i 1 
ancndhinanss 
R R 

1 
Scena 
RR’ 





tg 6=— (62) 


The levels themselves are 


E,=f+2a sin | a] (63) 
(N+3) 


E,=f—2a sin 6 (63’) 


and the transition energy 


2ar 2ar 1 
AE= cos 6= (1-- ) 
N+3 N+3 RR’ 


1\- +“ 
x(1+-,) (1+—) . (64) 
R? R” 


This is identical with the expression (52), 
which gives AE according to LCAO for R<1, 
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R’ <1, except that the absolute value is to be 
taken and that we have N+3 instead of N+1, 
which difference has no significance for the pres- 
ent accuracy. 

4. Finally there exists the possibility 


R<il R’>1 RR’ +1. 


Then the lowest state is given by (41), with v 
near R, or in the first approximation 


, 


v=R+R24+4 (1—R?), 
R 





, 


| 


_— R2N+2 





1 
\ R 


“ 





l 1 _ RR’-1 | 
E,=f+a)R+—— R**+? (1—R*)*. (65) 
ie R'-R 


The higher states are contained in (47), which, 
divided by R’+R, can be written 
R’ 


sin (N+2)k—R 





sin (V+3)k 
R 


, 





1 
sin (V+1)k=0. (66) 
R 


/ 


The approximate solutions are 
(N+2)k=rr+6 r=1---N+1, (67) 


of which there are the proper number, V+1. 
Therefore 


























R’ r N+3 
sin 6=R sin m+ 6 
R'+R N+2 N+2 
1 ‘Tr N+1 
= sin ; - 6;. (68) 
R+R |N+2 N42 J 
For large N or small R, 1/R’ 
RR’-1 r 
sin 6= sin TT, (68’) 
R’'+R N+2 
1K 
E=f+2a) cos 
1 RR’-1 1T 
— sin? . (69) 
_N+2 R’+R N+2 
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Therefore, according to HLPS, the first absorp- 
tion corresponds to the energy 
ra 
AE=E,-—EF2= — 2a cos ——— 
N+2 
1—R?*)? 


1 ( 
+a(R+—) +a| ee 
R’ R’—R 


1 1 
+ ( ) sin vd | RR’), (70) 
N+2 \R’+R N+2 


while the LCAO method gives (transition from 
r=(N—1)/2 tor=(N+1)/2) 








ra ra 
sin 
2(N+2) N+2 
1 RR’— 1 3m 


+ sin? - 
N+2 R+R’L 2(N+2) 


AE = 4a cos 











° 9 
—sin- 





T 
—| (71) 
2(N+2) 


Therefore AE approaches, for increasing N, a 
finite limiting value, if one or both of the R’s 
are smaller than 1 and if the HLPS method is 
used, in all other cases it goes to zero. 


IV. THE NEUTRAL MOLECULE" 


The neutral base of the dye has two ground 
structures, one of energy é1. 


H.N—(CH+=CH),—_-CH=NH ( 


and one of energy é2 


~~ 
bo 


H.N=CH—(CH=CH),—NH, (73 
| 


and N excited structures of energy f 


NH,—CH—(CH=CH),—NH. (74 


Formally, therefore, the situation is exactly the 
same as the one treated in III. For R, R’<1, the 
transition energy is given by (46), (46’) in HLPS 
and by (52) in the HHM method. 

For RR’ =1, it is (55) for HLPS and (57) for 
LCAO. 

For R, R’>1, it is expressed by (61) in the 
HLPS method and by (64) according to LCAO. 


16 See also U. Fano, Phys. Rev. 61, 736(A) (1942). 
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Finally, if R<1, R’>1 and RR’+1, AE is 
found from (70) for HLPS and from (71) for 
LCAO. 

Just as for the unsymmetrical ion, AE goes to 
a finite limit with increasing N only if at least one 
of the R is smaller than 1 and the HLPS method 


is used, in all other cases it goes to zero. 


V. A DIFFERENT ASSUMPTION ABOUT THE ION" 


In the preceding paragraphs, we have given to 
the structures of type (2) and type (3) the same 
energy f. It is possible that it would be a better 
numerical approximation to take the energy of 
(1) and (3) equal, namely e. In this case, the 
structures have alternately energies e and f, 
(N+3)/2 having e, (V+1)/2 the energy f. The 
equivalent differential equations are 


#1 = —ex,;— Bxo, (75) 
i2= —fr2— B(xs+x1), (76) 
a ae oY (77) 
#oe41= —€X2s41—B(X2e+X2542), (78) 
dys2= —eXn42—BXn41.- (79) 

For integration, one writes 
Xe =a yl(N+3)/ 2)-— si. abla dd: (80) 


=b{1 {lv4+8)/2)— 2s— sure iabeae (81) 


Ww sien the Eqs. (77), (78) in the form 


Xos-1+X2841 








E-f=s———_—__, (77’) 
X25 
' XoestXeos42 
E-e=6 : (78’) 
X2s4+t 


the insertion of (80), (81) and the requirement 
that the right side be independent of s, leads to 


Then (77’), (78’) take the form 


1 
B-s=9-(0+-), (82) 


1 
E-e= s-(0+-) (82’) 


" By mistake, the letter 8 has been used in this chapter 
for the quantity previously called a. 
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or 


1 2 
(E-f)(E-e) =8°( 0+ ) , (83) 


which replaces (15), (15’). 

Furthermore, writing (75) and (76) in the 
analogous forms to (77’), (78’) and multiplying 
leads to 

Xst+x1 
(E—f)(E-—e) =B*——. (84) 


x1 


The use of (80) and (83) then gives 


yl(N+3)/2)+1 4 g—((N +3) /2)+1) 





=—1. (85) 


yl(N+3)/2]—-1 4 y—((N +3) /2]—-1) 


This has no solution for a real v<1. Putting 
therefore again v=exp.(ik) one finds, for the 
symmetric solutions, 


N+5 N+1 
cos ( , bk) = cos - k (86) 








and for the antisymmetric ones 











N+5 N+1 
sin = —sin k, (87) 
2 
which leads to 
rr 
k=- , (88) 
N+3 


r odd for symmetric, even for antisymmetric 
solutions. 
A short transformation shows that 


a rr 
x;= | sin ( ), (80’) 
b N+3 


where the coefficient is a or 6, depending on the 
question whether j is odd or even. 
One sees that r=1---N+2, which at first 
glance seems to give the correct number of states. 
The energy follows from (83) as 





rr )} 


(e—f)?+ 168" cos* 
N+3 





(89) 


r 


We have now to investigate the sign in (89). 
Physically it seems clear that the positive sign 
should go with (N+3)/2 states, which lie near 
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e for small 8, the negative sign with the (V+1)/ 
states lying near f. 

But from the mathematical standpoint, the 
double sign seems to give, at first glance, two 
energy values for each state. 

We see clearer, if we remember that xj, i.e., 
the coefficients in the wave function, is not de- 
fined by the sine alone, but also by the ratio a/d. 

Compare the state 7 with the state N+3 
—r=r'. 


We find 















a 7 
x= sin [+ —w+3-n] 
b N+3 





a 
= (—1)*'sin (90) 









That is to say: If a/b has the same sign in 7 and 
r’, the phase of every alternate “‘particle”’ differs, 
in the state 7’, by 180° from the phase of the same 
“particle” in the state r (analogous to the “‘op- 
tical branch”’ in crystal vibrations). If, however, 
a/b has the opposite sign for 7 and 7’, the alter- 
nate “‘particles’’ have the same phase in 7 and r’ 
(acoustical branch). In other words, what matters 
is the ratio 










Tr} ™r 
: b sin 
N+3 
In the energy formula (89), cos 2r/(N+3) 


has the same value for 7 and 7’. 
It follows from (75), (89) that 








(j+1), j=2s41. 


a sin 






x f-e 
E-—e=—— 





x1 2 






(91) 





1 as 
+-—} (e—f)*+16? cos? 
2 N 





Therefore, with the positive sign, x2/x1 is posi- 
tive; with the negative sign, negative. To sum up: 


r need only be taken from 1 to (V+3)/2. 
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For r=(N+3)/2, we have only one solution, 
namely (see 82’) 


E=e a=0 (positive sign in (89)). 


The negative sign in (89), giving E=f, to 
which would correspond }=0, does not satisfy 
Eq. (82). To E=e belongs the wave function 


Xe41=(—1) 
Xo,= 0. 


For each r which obeys 


N+1 


’ 


2 


1 


IIA 


r= 





there are two solutions. 

The positive sign in (89) gives an energy larger 
than e. The corresponding wave functions have 
the character of the ‘‘acoustical branch.” 

The negative sign in (89) gives an energy 
smaller than f, the corresponding wave functions 
resembling the ‘‘optical branch.” 

The energy difference between the two lowest 
energy values is 


T } 
N+3 


2r 3 
— | (e=$)*+168 cos? —_| 
N+3 





1 
az=| (e—f )*+ 168? cos? 


1 ae 
=| | (e=s)*+168°— 168° sin ented -| 
7 N+3 


| Qn )? 
~} (e—f)?+166?—16 sin? — . (92) 
| N+3 


which, for large N, takes the approximate form 
3x? 


(N43)? 





AE=46?| (e—f)?+1687}—* (92’) 


This investigation was started at the sugges- 
tion of A. L. Sklar and has profited by many 
discussions with him. 





AUGUST, 1942 JOURNAL OF CHEMICAL PHYSICS VOLUME 190 


Energy Levels and Color of Polymethine Dyes 


A. L. SKLAR 
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(Received April 29, 1942) 


The secular determinant for the energy levels of the unsaturation electrons of a polymethine 
dye is discussed in both HLSP and LCAO approximations. The roots of the secular determi- 
nant, which were obtained in the preceding paper by Herzfeld, are applied to a discussion of the 
longest wave-length electronic band of symmetrical and unsymmetrical polymethines. The 

*LCAO approximation gives good numerical values for the dependence of both the transition 
energy and oscillator strength on the length of the polymethine chain in symmetrical dyes, 
but cannot handle the questions which depend sensitively on a small difference between the 
groups attached to the two nitrogen atoms at the extreme ends of the chain. Although the 
HLSP method yields correct qualitative results for symmetrical dyes, it does not give good 
numerical values. The HLSP method, however, is sensitive to a difference between the two 
ends of the dye molecule and affords a simple explanation for a number of properties in which 


unsymmetrical dyes differ from symmetrical ones. 





HE spectra of the polymethine dyes have 

been discussed briefly from the atomic 
orbital point of view by Pauling! and from the 
molecular orbital viewpoint (LCAO) by Mulli- 
ken? and Foérster.* The connection between the 
molecular structure and the low lying energy 
levels of these dyes will be discussed qualitatively 
in a forthcoming paper;* here it will be treated 
more quantitatively in order to arrive at a firmer 
basis for the purely qualitative considerations. 
The lower energy levels of the plane ion I will be 
approximated by both the atomic orbital (HLSP) 
and molecular orbital (LCAO) methods. Ion I 
represents a large class of polymethine dyes, if 
the T groups are 


T; ° T> 


\ / 
N+=CH—(CH=CH),—N’ 


% 
Ty’ T! 


understood to mean either individual groups or a 
nucleus which bends around and joins to the 
conjugated chain, as in the ion II. 


'L. Pauling-Gilman, Organic Chemistry, Vol. 2, p. 888. 

*R. S. Mulliken, J. Chem. Phys. 7, 570 (1939). 

* T. Forster, Zeits. f. physik. Chemie B47, 245 (1940); 
48, 12 (1941), 

*A. L. Sklar and L. G. S. Brooker (to be published). 
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Since the excitation of the in plane o electrons 
yields the molecular Rydberg series below 2000A, 
we may restrict our discussion concerning the 
visible and near ultraviolet absorption to the 
energy values of the unsaturation electrons (z or 
P.) which do not interact with the o electrons. 

In the LCAO method, as is well known, we 
write a linear combination of atomic orbitals for 
the molecular orbital of each unsaturation elec- 
tron in the field of the multiply-charged residue 
formed by stripping ion I of its unsaturation elec- 
trons. Ion I has (2n+4) unsaturation electrons 
but we will consider that two of these are perma- 
nently localized, one on each nitrogen atom. Our 
system, then, consists of (27+2) electrons in a 
potential field made up of a chain of (2n+3) 
singly charged nuclei. If we write the one-electron 
orbitals as a linear combination of the (2n+3) 
atomic orbitals and neglect all overlaps and also 
interactions between electrons on non-neighbor- 
ing atoms, the energies of a single electron are 
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given by the roots of the (2n+3) rowed secular determinant III, whose solution has been dis- 


cussed in the preceding paper.’ 





la-W; am; O ; 0; 
a ; f-W; a; 0; 
0 ; a ; f-W; a; 
™*, 
0 » 3 a oe 
‘ ‘\ 
; 
0 ; 0 ; S 3 
® . % 0 ; 
| Oo; O; 0; 


In this determinant W is the one-electron en- 
ergy; a is the interaction energy between two 
electrons on neighboring carbon atoms; a; and 
a are, respectively, the interactions between two 
electrons on a carbon atom and on either the left- 
or right-hand nitrogen atom; and f, e:, and eg are, 
respectively, the ionization potentials of an un- 
saturation electron localized on either a carbon 
atom, the left- or right-hand nitrogen atom. In 
this approximation the f’s will be somewhat 
different for the various carbon atoms in the 
chain. It will, however, be sufficiently accurate 
for our purposes to take all the f’s as equal,® 
especially since the variation in f in this approxi- 
mation is due primarily to the fact that the inter- 
action between the unsaturation electrons is 
neglected. In the ground state of the ion I the 
(2n+2) electrons will fill the lowest (n+1) orbi- 
tals and the longest wave-length electronic band 
will be associated with the transition of an elec- 
tron from the (n+1)st to the (n+2)nd orbital. 

In going over to the Heitler-London-Slater- 
Pauling method it is to be noticed that the ground 
structure I is doubly degenerate in that the posi- 
tive charge may be shifted to the other nitrogen 
atom and the double bonds readjusted. Ion I is 
thus analogous to benzene which has the two 
Kekulé structures as doubly degenerate ground 
structures. A difficulty, however, arises in carry- 
ing over from benzene to ion I the view that the 


5K. F. Herzfeld, J. Chem. Phys. 10, 508 (1942). 

6 We also neglect energy differences due to the replace- 
ment by other groups of the hydrogen atoms attached to 
the conjugated carbon chain. 
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longest wave-length electronic absorption is asso- 
ciated with a transition between the two molecu- 
lar states which arise from a resonance splitting 
of the degenerate ground structures.’ The elec- 
tron distributions in the two I structures, which 
differ in that an electron has moved from one end 
of the molecule to the other, do not overlap ap- 
preciably. The split of these two structures would 
thus be expected to be negligible instead of the 
observed value of the order of a few volts. 

Pauling,’ however, has pointed out that one 
must consider, in addition to the two structures 
I, the following set of intermediate structures 
X 2p41(p =0 to n). Since the overlap of electrons on 
atoms separated by even two interatomic dis- 
tances is very small, we will also have to include 
the set of structures Yo,(p=1 to n) 


T; 
\ + 
N—(CH=CH),—CH— T» 
A ff 
Ty’ _(CH=CH),_»—N 
Xops1 T;! 
T) 
\ 
N—(CH=CH),-1—CH— J 
ry : / 
—CH—CH—(CH+=CH),_,—N 
\ 
Yo T;! 








7A. L. Sklar, J. Chem. Phys. 5, 669 (1937); H. Sponer, 
G. Nordheim, A. L. Sklar, and E. Teller, J. Chem. Phys. 
7, 207 (1939). 

8L. Pauling, Proc. Nat. Acad. Sci. 25, 577 (1939). 
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For every structure in the set consisting® of I, 
X, and Y there is a second structure which differs 
from the former in that an electron has been 
moved to the neighboring atom on the left. Thus, 
although the two I structures do not overlap 
appreciably directly, one I structure interacts 
with Xj, X; with Yo, Y2 with X3, and so on until 
Xen41 interacts with the other I structure. The 
two ground structures I can split through this 
chain of interactions with the intermediate struc- 
tures X and Y. 

Although it is not easy to estimate numerically 
the relative energy of I and X, it is clear from 
chemical considerations” that the I structures are 
considerably lower in energy than the X and Y 
structures. The Y structures are energetically 
higher than the X structures since the former 
have one less double band. Acting against this 


TABLE I, 








HLSP LCAO 





R&I. 2a(1 —R2)R (21) 2esin aT ; (26) 


R=1 4a@ sin? [r/2(N+2)] (30) 4a cos (31) 


—_ T 
N +2 °" (N42) 


3 of , 7 ” 
i a 2 —— 
R> 1 4a sin 2Na3 (36) 2a sin Ni3 (37) 








* The set of structures X and Y are chosen in preference 
to the following set, Z.(k=1 to 2n+1), on energetic 
grounds. 


T; 
N=CH—(CH=CH),-,—CH— T: 
oe 
Tr,’ —(CH=CH),:—CH=N 


* 
Zap. T2’ 


It is not obvious, though, that the interaction energy 
among the Z structures is not sufficiently larger than that 
among the X, Y structures to compensate for the higher 
energy of the former. This question is, however, immaterial 
for the purpose of this paper since we will not try to 
evaluate the interaction energy but only its variation as n 
is increased or as the character of the T groups is changed. 

10 The essence of the chemical arguments is the fact that 
whereas ammonium ions N* are frequently encountered, 
carbonium ions C* are never met in any quantity, although 
they probably do occur in minute amounts in the course 
of reactions. As an example, we may compare the basicity 
of ammonia (NH;+H:O—NH,*++OH~-) with that of 
ethylene (CH2—CH2+H:0—CH;—CH.2*+0H7). Except 
for the very small difference between the strength of a CH 
and NH bond, the second reaction requires more energy 
than the first by about the same amount that is required 
to change structure I into an X structure. (The hydration 
energies balance out to a first approximation.) The fact 
that ethylene is enormously weaker as a base than ammonia 
would, then, suggest that X is considerably higher in 
energy than I, since there is no reason to suspect the 
entropy to be vastly different in the two cases. 


effect is the additional possibility of resonance 
possessed by the Y structures due to the fact that 
the C+ breaks the carbon chain into two odd 
segments which may be written several ways, as 
for example: 


—(CH=CH),—CH—=—CH—(CH+=CH),—. 


To each Y, then, corresponds a group of struc- 
tures and we may consider that, in writing one 
structure for this group, we have written the 
stabilized result of the resonance among the 
whole group. Since two compensating factors 
enter in the relative energy of the X and Y struc- 
tures, and since a small difference in their energy 
will not affect our discussions, we will assume, for 
mathematical: simplicity, that all intermediate 
structures, X and Y, are of equal energy. 

In all we have 2n+3 structures, two I, (n+1)X, 
and mY structures. If we neglect all Coulomb and 
exchange integrals which are smaller than the 
product of the first power of the overlap integral 
of electrons on neighboring atoms times the or- 
dinary Coulomb and exchange integrals, respec- 
tively, and also neglect terms containing overlaps 
of electrons on non-neighboring atoms, the mo- 
lecular energies resulting from resonance among 
the 2n+3 structures are given by the roots of a 
secular determinant which is formally the same 
as that reached in the LCAO method, namely 
determinant III. The meaning of the quantities 
in the determinant is, of course, now quite differ- 
ent. a, a1, @ are, respectively, the interaction 
energies between an intermediate structure, and 
either a second intermediate or one of the two 
ground structures; e; and é2 are the energies of 
the two ground structures I; and f is the energy 
of an intermediate structure, X and Y. Since W 
is the energy of all the unsaturation electrons in 
the molecule, the transition energy in question is 
the difference between the two smallest values 
of W. 

The roots of determinant III have been dis- 
cussed as a function of the number of rows in the 
preceding paper by Herzfeld’ for various ranges of 
the parameters R=a/(e,—f) and R’=a/(e.—f). 
The results will now be applied to a discussion of 
the qualitative conclusions which were reached 
in reference 4, since purely qualitative considera- 
tions of complex situations, even when applied 
with caution, are occasionally misleading. 
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Fic. 1. The solid points in the figure give the wave- 
length of the absorption peak for ion II as a function of 
the number of double bonds between the nitrogen atoms 
(NH). The empty circles give the same data for an ion 
which differs from ion II only in having an acetoxy group 


VA 
(CH;—C—O—-) in place of the hydrogen atom on the 
sixth carbon atom of the polymethine chain. 


SYMMETRICAL IONS 


We will first treat those dyes in which the same 
groups are attached to both nitrogen atoms as, 
for example, in ion II. For such symmetrical ions 
the secular determinant is simplified since a;= az 
and ¢;=é2. Since a and a; are of the same order 
of magnitude and since a small difference be- 
tween them would not be expected to affect the 
general questions treated in this paper, except 
for very small values of 2, we will further sim- 
plify the mathematical treatment by setting a; 
equal to a. 
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Table I gives the transition energies to the 
first order found in reference 5 from the roots of 
the secular determinant, for various ranges of 
the parameter R=a/(e—f). In Table I the num- 
bers in parentheses are the numbers of the equa- 
tions in reference 5 and N=(2n+1) is the num- 
ber of carbon atoms in the chain between the two 
bounding nitrogen atoms. 

It is clear from the formulas in Table I that 
the magnitude of R plays a determining role in 
the HLSP but not in LCAO approximation, ex- 
cept for very small values of NV. This becomes 
clear, physically, if we regard a as constant and 
vary R by varying e—f. In the HLSP method the 
energy of the degenerate ground structures, e is 
the zero approximation to the energies of the 
two lowest states of the molecules. Therefore, 
the transition energy depends on the magnitude 
of the split of the degenerate energies of the 
ground structures. Since these two structures 
interact only through a chain of second-hand 
interactions with the intermediate structures X 
and Y, which are higher in energy than the 
former by an amount (f—e), it is clear that!! the 
greater is (f—e), the smaller is the transition 
energy. Table I of reference 5, obtained by a 
numerical solution of Eqs. 17 and 18, shows this 
dependence if a is considered to be constant; the 
transition energy increases as (f—e) decreases 
(R increasing) and approaches for small values 
of (f—e) the limiting value, independent of R, 
required by Eq. 36. 

In the LCAO method, on the other hand, e is 
the zero approximation, not to the energy of the 
molecule, but to the lowest energy state of one 
electron of the molecule. Since the states associ- 
ated with e, as well as a number of states whose 
energies in zero approximation are f, are full, 
both in the ground and excited states, there is no 
longer any reason to expect (f—e) to play a large 
role in the transition energy. 

It is clear from Table I that in both approxi- 
mations and for all ranges of R, the transition 
energy should decrease monotonously with in- 
creasing chain length Vas is found to be the case.” 

11 f—e may be changed experimentally by changing the 
character of the groups 7; which are attached to the 
nitrogen atom. For example, the ionization energy of a 
electron from methyl amine, (CH3)-NHz is smaller than 
that of aniline, (CsH;)2NH. 


2L. G. S. Brooker et al., J. Am. Chem. Soc. 62, 1116 
(1940). 
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The algebraic relationship between the peak 
wave-length and WN differs in the two methods 
and, in the HLSP approximation, depends on R. 
The usual linear dependence of peak wave-length 
on chain length shown in Fig. 1 agrees with the 
results found for large values of N by the LCAO 
approximation for all ranges of R. In the HLSP 
approximation the dependence on N is quadratic 
for RS 1, but is exponential for very small values 
of R; only in the region of R slightly smaller than 
unity is there a possibility of a linear dependence. 

The value of the absorption frequency for fixed 
N is directly proportional to the interaction 
energy a for constant R, but the ratio of the fre- 
quencies of a dye with a given value of N and 
that with an NW larger by one, is not explicitly 
dependent on a. The value of this ratio affords a 
test for the two methods. In the LCAO method, 
the ratio does not depend sensitively on the value 
of R and is therefore completely determined by 
N; in HLSP approximation, the ratio does de- 
pend on R. We can compare the theoretical 
values with the experimental" value of 1.17 (cf. 
Fig. 1) for the ratio of the peak absorption of a 
dye with a chain of seven carbon atoms (n= 3) to 
that of one with a chain of nine carbon atoms 
(n=4). The LCAO approximation, in excellent 
agreement with experiment, yields the value 1.19 
for this ratio when R is very large; 1.20 when R 
is unity, and 1.24 for very small values of R. In 
the HLSP approximation, however, the above 
ratio varies between 1.2 and 1.4 as R varies 
between one-half and ten.% 

In addition to yielding too high a value for the 
above ratio, the HLSP method demands an ab- 
normally high value for the exchange energy a 
to explain the absolute values. The transition 
energy for a five-membered carbon chain (n= 2), 
which is actually about three volts, is only two- 
tenths a in the region R=1, which seems to be 
the most likely region in this method. In the 
LCAO approximation the transition energy, 


The data in Fig. 1 were obtained at the Eastman 
Kodak Laboratories and communicated to the author by 
Dr. L. G. S. Brooker. 

“It has been pointed out (reference 12) that the wave- 
length of the peak of most symmetrical polymethine ions 
increases by about 1000 angstrom units when m increases 

y one. The constancy of this value is to be expected 
according to Table I except in the HLSP approximation 
when RK1. 

® Reference 5, Table I. 
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roughly independent of R, for a five-membered 
carbon chain, is about nine-tenths a, which is 
much more reasonable. 

It is seen that although both methods tell us 
that the absorption frequency should decrease as 
the chain lengthens, the quantitative features of 
this decrease are not given very well by the HLSP 
method, but are given very well indeed by the 
LCAO approximation. 

We may also inquire into the dependence of 
the intensity on chain length. Although the chain 
is a bent one, we may assume that each increase 
of m in ion I by unity lengthens the chain by a 
constant amount. The coefficients of the orbitals 
in the LCAO approximation are given, except for 
a normalization factor, in Eqs. 13 and 14 of 
reference 5. With them we can calculate the 
transition moment, which, as Mulliken? has 
pointed out, is polarized along the chain. In first 
order, the transition moment turns out to be 
proportional to V when N is large. Since the in- 
tensity is proportional to the product of the 
frequency and the square of the transition mo- 
ment, and since the frequency in the LCAO ap- 
proximation is inversely proportional to N, the 
intensity should vary linearly with N. The in- 
tensity of the shorter members of the series 
represented by ion II (small values of 2) is shown 
in Table II.1® Although the first member is out of 
line, the last three fp values vary linearly with the 
chain length as indicated by the last two succes- 
sive differences. Since the theoretical discussion 
on the variation of intensity with chain length 
was limited to large values of N, the excellent 
agreement of the last three values is probably in 
part fortuitous. It would, accordingly, be useful 
to obtain f values for the higher members of this 
series but these, unfortunately, are not very 
stable and rapidly change with time. 

In addition to yielding the correct dependence 


TABLE II. 

















fo is the oscillator of the longest wave-length electronic band of ion 
II and N =2n-+1 is the number of carbon atoms in the chain. 


16 The data in Table II were calculated from the absorp- 
tion curves of Fig. 2. 
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of intensity on chain length, the LCAO method 
also has yielded interesting results concerning the 
relative intensities of the set of absorption bands 
of a given dye. On the basis of the LCAO method, 
Mulliken has reported? that in a series of N-V 
transitions, as are our transitions, in a maximally 
elongated chain, the longest wave-length transi- 
tion is by far the strongest. Since the ions II have 
bulky nuclei on the ends, their chains may be 
expected to be in an elongated configuration. 
Figure 2 shows,’ in accord with Mulliken’s pre- 
diction, that the higher energy transitions are 
very weak compared with the lowest energy 
transition. For example, the oscillator strength 
of the band having a peak at 3600A (N=9) is 
0.23, the long wave band having an oscillator 
strength of 1.85. 


UNSYMMETRICAL IONS 


If the nuclei 7; which are attached to the right- 
hand nitrogen atom in ion I, differ from those 
attached to the left-hand one, the ionization 
energies of the two nitrogen atoms are no longer 
equal and ¢é; is, in general, different from e2 in 
both methods. In the HLSP approximation; we 
would expect a sensitive dependence of the color 
on the difference between the two e’s since the 
resonance splitting of the energies of the two I 
structures is both affected by, and superimposed 

17 The absorption curves of Fig. 2 were taken at the 


Eastman Kodak Laboratory and communicated to the 
author by Dr. L. G. S. Brooker. 
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Fic. 2. The absorption curves of ion II for n=1 to 4; N=3 to 9. 
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on, the energetic difference (e:—é2). In the LCAO 
approximation, on the other hand, this difference 
clearly plays a minor role in the transition energy. 
The special sensitiveness of the HLSP method to 
(€,—€2) is due to the fact that in this method e:, 
é2 are the zero-order approximations to the ener- 
gies of the two levels involved in the transition; 
whereas, in the LCAO approximation, the two 
one-electron levels approximated in zero order 
by é1, é2 are filled in both the ground and excited 
state of the molecule. 

That differences in kind between symmetrical 
and unsymmetrical ions do exist may be seen 
experimentally by comparing the absorption 
peak of an unsymmetrical ion with the mean of 
the absorption peaks of the two parent sym- 
metrical ions of which the unsymmetrical ion 
may be considered a hybrid. The existence of 
important differences between the absorption 
peak of unsymmetrical ions and the mean of 
those of its two parents has been emphasized by 
Brooker'® and called the “deviation.” Isolated ex- 
amples! of the deviations are given in Table III. 

As discussed above we can only hope to get an 
explanation of the significant differences between 
unsymmetrical and symmetrical ions through the 
HLSP approximation. However, before any con- 
fidence can be placed in results which are yielded 
by the HLSP but not by the LCAO method, one 


18 L. G. S. Brooker and Sprague, J. Am. Chem. Soc. 63, 
3202 (1941). 
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TABLE III. % is the wave number of the absorption peak of the ion in column I, %,(%g) is that for the symmetrical ion 
having the nucleus on the left (right) in common with the unsymmetrical ion, and the last column is the deviation. All 
ions occurring in this table have N=5. Me is a CH; group, Et a C2H; group, and Ph a phenyl group. It is to be noted 
that in the ions given in the last two lines of the table the end groups containing conjugated double bonds are joined to 
an inner member of the chain. It is not clear how far the present formulae apply to this case. However, the table shows 


that the rules discussed do apply and that the deviation is particularly large. 








Ion 


Dp 





(Me): 

Cc S 

-_* oN, ) 
C—CH=CH—CH=C 


‘N—Ph 


C—CH=C vy =C 


on, 


18,150 


17,140 


23,950 


22,730 17,510 








what closer to the LCAO as judged by the 


. of must ask whether the results are artefacts of the 

m- HLSP approximation, especially since the LCAO numerical results concerning the dependence of 
ion method gave better numerical results on sym- the transition energy on chain-length. Since the 
: of metrical dyes. actual situation is intermediate between that de- 
ion That one can trust, in a qualitative manner, — scribed by the two methods, the approximation 
_ of the conclusions which the HLSP method gives which overemphasizes a given effect is certainly 
| by concerning the effect of the difference (e:—e2) can more reliable for finding, in a qualitative way, 
ex: be seen by comparing the direction of the errors the consequences of that effect than is the ap- 
III. made by the two methods. As discussed before, proximation which practically neglects the effect 
t an €1, €2 affect the transition energy sensitively in in question.!® 

reen HLSP because they are zero-order approxima- Our discussion of unsymmetrical ions, e:+ é2, 


the 
con- 


ided 


one 


Se 63, 


tions to the energies of the two molecular states 
involved in the transition; whereas, in the LCAO 
approximation, the difference (e:—e2) has no ap- 
preciable effect on the molecular states. The true 
molecular state, however, is somewhere between 
those described by the two methods, being some- 


19 Since, in the HLSP approximation, ¢, e2 are important 
primarily because the two structures I have more chemical 
bonds, one might hope that the dependence of the molecu- 
lar energies on (€:—é2) will reappear in the molecular 
orbitals method if antisymmetrical molecular orbitals, 
including spin, are used and the interaction of electrons 
explicitly included in the Hamiltonian. See M. Goeppert 
Mayer and A. L. Sklar, J. Chem. Phys. 6, 645 (1938). 
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will accordingly be limited to the HLSP approxi- 
mation and the results will depend on the magni- 
tude of R. 

The same factors which cause e; to differ from 
é2 also introduce a difference between the inter- 
action energies a; and a2. However, a distinction 
between the roles played by an inequality of e; 
and é2 and that played by an inequality of a 
and a2 appears if one compares an unsymmetrical 
dye with the two parent symmetrical dyes of 
which the unsymmetrical dye may be considered 
a hybrid. Both ground structures are degenerate 
of energy é: in one parent symmetrical dye and 
of energy é2 in the other parent; whereas, in the 
unsymmetrical dye the ground structures are no 
longer degenerate, but have the energies e; and 
és in first order. Since the e’s are zero approxima- 
tions to the energies of the molecular states in- 
volved in the transition, we might anticipate, as 
will be shown later, that the transition energy of 
the unsymmetrical dye will be larger than that 
of either parent symmetrical dye. That is, on the 
basis of a difference in e; and é2, we should expect 
the unsymmetrical dye to differ from any mean 
of its parent symmetrical dyes. On the other 
hand, if e; were equal to é., even if a, were not 
equal to a2, we would expect the properties of 
the unsymmetrical dye to be intermediate be- 
tween those of its symmetrical parents. Further 
theoretical investigation will be required to de- 
cide which type of mean one should expect if 
a, differed from a2 but if e: were essentially 
equal to é2. 

The properties, then, which are peculiar to un- 
symmetrical dyes, and which distinguish them in 
kind from symmetrical dyes, may be found ex- 
perimentally by comparing the absorption peak” 
of the unsymmetrical dye with those of the sym- 
metrical parents, and may be found theoretically 
by setting a1= a, but keeping e; different from es. 

We will, then, set a; equal to ae,”4 but keep e1 
different from é2, in the secular determinant III, 
and thereby obtain a transition energy for the 
unsymmetrical dye which is to be contrasted 
with a mean of the transition energies of its two 


20 Since the absorption curves are obtained in solution 
the significance of the wave-length of peak absorption is 
open to question. This can perhaps be settled by a study, 
now under way of the structure of the absorption bands at 
low temperatures. 

a in the case of symmetrical dyes, we also equate a 
and a. 
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symmetrical parents. The results depend upon 
the order of magnitude of R=a/(ei:—f) and 
R’=a/(e.—f). 


Case I. R<1 and R’<1 


The easiest case to discuss qualitatively is that 
of small values of R and R’ which is the case when 
the intermediate structures, X and Y, are con- 
siderably higher in energy than the ground struc- 
tures I; the qualitative discussions of reference 4 
were limited to this case. Equation 46 of reference 
5 gives the transition energies when R and R’ are 
small as: 


2AEn = { (€1—e2)?+(aon)?}? 
+ { (€:—e2)?+ (aoy’)}?, 


oy? =4R?% (1— R?)(1—RR’). 


(46) 
where 


oy’ has R and R’ interchanged. acy, acy’ differ 
from the transition energies of the parent sym- 
metrical dyes only in that the factor (1— RR’) 
replaces a factor (1— R®) which is a difference of 
higher order for the case of small R under 
discussion. 

A number of interesting results follow from 
Eq. 46. First, it is clear that the effect of a small 
non-degeneracy (€:—é2) causes the transition 
energy AEy to be larger than (oy+oy’)a/2, 
which is essentially the mean of the transition 
energies of the parent symmetrical dyes. Thus, 
it is seen that unsymmetrical dyes should absorb 
at shorter wave-lengths, that is, be ‘‘lighter”’ in 
color, than related symmetrical dyes. This has 
been found by Brooker” who, in a series of papers, 
gives values of the ‘‘deviation,’”’ which are always 
found to be toward short waves except in one or 
two cases when they are extremely small. Ex- 
amples of the deviation are shown in Table III. 

Although the transition energy is larger, the 
resonance stabilization in the ground state is for 
unsymmetrical dyes smaller than the mean of 
their related parent dyes. This may easily be 
seen from Eqs. 20 and 45 of reference 5, by com- 
paring the difference between the energy of the 
lower ground structure and that of the ground 
state for the various cases. One should accord- 
ingly expect that unsymmetrical dyes, by virtue 
of their larger resonance stabilization, should be 

2L. G. S. Brooker et al., J. Am. Chem. Soc. 62, 1116 


(1940) ; cbid., 63, 3192 (1941); ibid., 63, 3203 (1941) ; idid., 
63, 3214 (1941) ; zbid., 64, 199 (1942). 
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Fic. 3. ‘Converging’ and 
“non-converging”’ series. The 
dye in the center of the figure 
differs from the one at the top 
only in that the acetate radicle 


Et 


ots 


O T(-NHPh) 
CH cf $ 
~ = 
Omen 

is replaced by a hydrogen atom. 


more stable than unsymmetrical dyes. They 
should, for example, be less readily decomposed. 

An interesting confirmation has been called to 
the author’s attention by Brooker. The cyanine 
dyes act as indicators both in an acid and alka- 
line range. Now, either the addition or subtrac- 
tion of a proton will be resisted by the resonance 
stabilization of the cyanine chain since the reso- 
nance is present only before either of these 
processes. On this basis of a larger resonance 
stabilization of symmetrical dyes, one should 
expect the applicability to the polymethine dyes 
of a rule which has been found by Schwarzen- 
bach* and which states that the difference be- 
tween the pK values” (pH at half-transformation) 
is larger for symmetrical than for unsymmetrical 
dyes. 

It may also be seen from Eq. 46 that as N in- 
creases, ov, and accordingly AEy, decreases. Now 
in a symmetrical dye AEy would go to zero as N 
is increased indefinitely, but in an unsymmetrical 
dye, AEy decreases to the magnitude of the 
energy difference (e:—e2). It has, indeed, been 
found by Brooker” that for many series of un- 
symmetrical dyes the members of which differ 
only in the value of NV, the absorption peaks ap- 
pear to converge to a finite wave-length limit as 
N increases, whereas the peaks for a correspond- 
ing series of symmetrical dyes do not converge 
to a finite line when plotted on a wave-length 

3G. Schwarzenbach, Zeits. f. Electrochemie 47, 40 (1941). 
Schwarzenbach’s examples are not strictly analogous to 
our case but will be shown in a later paper to have the 
same basic explanation. 

** Since the difference in the pK’s is a measure of the 
free energy change involved in the reaction {2I-—HI*++B} 
where B is the neutral molecule obtained by taking a 
proton away from ion I, one would not expect the entropy 
change in the case of unsymmetrical dyes to differ markedly 
from that of symmetrical dyes, except perhaps for an effect 
due to the symmetry number. The symmetry number 
comes into question because both HI** and B may be 
unsymmetrical whether I is symmetrical or unsymmetrical. 


This effect, however, is at most small (of the order of a few 
tenths of a pH unit), and even in the wrong direction. 
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scale. Many unsymmetrical dyes form series 
which behave like series of symmetrical dyes; the 
former presumably have small values of the 
difference (e:—e2). Figure 3 shows an example of 
a ‘‘non-converging”’ series (middle case) and two 
examples of ‘‘converging”’ series. The dye at the 
bottom of Fig. 3 is not an ion but its analogy to 
the present case is discussed in the last section. 

Equation 46 also affords an explanation of the 
“sensitivity rule’ found by Brooker,” which 
states that the deviation produced by introduc- 
ing a given alteration in the molecular structure 
of a symmetrical molecule is much smaller than 
the change in deviation produced by making the 
same structural alteration in a highly unsym- 
metrical molecule. This may be seen to be a con- 
sequence of the fact that the resonance inter- 
action and the energetic difference (e:—e2) enter 
into the transition energy as a sum of squares. 
When the difference (e:—e2) is small compared 
to acy, a given increase in (€,—é2) affects AEy 
much less than would the same increase if (e; — é2) 
were of the order of magnitude of acy or larger. 
One set of examples from Brooker’s papers is 
given in Table IV to illustrate the “sensitivity 
rule.” 


Case II. RR’ =1 


In order to see whether the results discussed in 
Case I are valid even if the intermediate struc- 
tures are not very much higher in energy than 
the ground structures, we may consider an un- 
symmetrical dye in which RR’ =1. This case can 
be solved exactly; Eq. 55 of reference 5 gives the 
following transition energy when RR’=1: 


(55) 


T 
AEy =a} R+R’—2 cos 
N+2 


One should expect, even in Case II, a finite 
convergence limit for the absorption peaks of a 
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TABLE IV. 7 is the wave number of the absorption peak in methyl alcohol of the dye in column I, %z, 5g are those of 
the parent symmetrical dye, which have both ends like the left (right) side of the ion in column I; Dev. is the deviation 
and AD is the change in deviation produced by introducing a nitro group. All ions in Table IV have five CH groups be- 
tween the nitrogen atoms. Note that the ion in the second row differs from that in the first row only in that a hydrogen 
atom is replaced by a nitro group; the fourth ion differs from the third in exactly the same way, but because the third 
ion is unsymmetrical (deviation 55), the effect of the introduction of the NO» is much greater than for the symmetrical 





ion in the first line. 
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homologous series of dyes as N increases in- 
definitely, though a smaller one than in Case I. 
This may be seen from equation 55’ which applies 
to Case II. 


Lim AEy={(ei—f)!—(ex—f)4]*%. (55’) 


For values of (e:—e2) which are small compared 
to (e,—f) this limit reduces to: 


(€:—€2)"/4(e1—f), 


which is smaller than (e;—e2), (Case I). 

The finite convergence limit implies that one 
should expect that unsymmetrical ions, even in 
Case II, should exhibit a ‘“‘deviation,”’ at least for 
large values of N. 

That the “sensitivity rule’’ should hold even 
in Case II may also be seen by considering the 
convergence limit rather than the deviation. The 
energy difference (e;—é2) is a function of e; since 





RR’ =1. If we now differentiate the convergence 
limit (R+R’—2) with respect to e:, we find that 
the derivative is proportional to the energy differ- 
ence (€:—é2). That is, we should expect that 
there exists a ‘‘sensitivity rule’ for convergence 
limit which would seem to imply the same for the 
deviation, at least for large values of N. 

The discussions concerning resonance stabiliza- 
tion cannot be treated in the same way as in the 
preceding section because only the symmetrical 
case in which R=1 has been solved exactly. Since 
the cases when R is slightly above or below unity 
have not been treated exactly, it is not simple to 
compare the behavior of the unsymmetrical ion, 
RR’ =1, with the mean behavior of the two sym- 
metrical parent ions. 

That the resonance stabilization is less for un- 
symmetrical than for symmetrical ions, even in 
Case II, can be seen, however, in the following 
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way, at least when R does not differ too much 
from R’. In this case, although the behavior of 
an unsymmetrical ion cannot be compared with 
the mean behavior of its symmetrical parents, it 
can be compared to the behavior of a symmetrical 
ion having an R equal to the mean of those of its 
symmetrical parents since if RR’=1 and the R’s 
are almost equal, the mean of R and R’ is unity. 
The resonance stabilization, E;— 1, is equal to a 
in the symmetrical case, R=1 (cf. Eq. 54, refer- 
ence 5), and equal to R’a in the unsymmetrical 
case, RR’=1 (cf. Eq. 54, reference 5) where R’ is 
the smaller of the two R’s and hence less than 
unity. Since the resonance stabilization of un- 
symmetrical ions is less than that of symmetrical 
ions when R and R’ are not very different, one 
would expect this same to be true, a fortiori, 
when the R’s do differ considerably. 


Case III. R>1; R’>1 


In case both R and R’ are large, the intermedi- 
ate structures are close to the ground structures 
and one should expect this case to differ in kind 
from Case I. From chemical considerations” one 
would not expect this case to obtain in the poly- 
methine dyes. 

Equation 61 of reference 5 gives the transition 
energy in Case III as: 


nr \? 1 1 2 
et) fie) 
N+3 R R RR 


1 1 ; \* 
x(1-=-—+ ) f (61) 
R R RR 


In this case, Eq. 61 shows that there is no 
finite convergence limit to the wave-length as NV 
increases indefinitely. This is not surprising since 
the ground structures, whose energetic difference 
caused the convergence limit in Cases I and II, 
are now very bad approximations to the two 
lowest states of the molecule. 


Case IV. R<1, R’>1 


This case is useful in considering the neutral 
dye molecule in which the lowest structure is IV,: 


qT; 











N—(CH=CH),—CH=N—72. IV, 
Ty’ 


A number of intermediate structures exist of 
which the lowest in energy is probably IV,: 


T; 
\ + - 
N=CH—(CH=CH),—N—T7>. IV>p 
Ty’ 


The remaining (2n+1) intermediate structures 
are: 


Ti 
\t ~ 
N=CH—(CH=CH),—CH— 
“A 
Ty’ —(CH=CH)n_ (p41) -CH=N—7? 
Wavps 1): 


Since the set of (2n+3) structures IV,, IV, 
and W; interact among themselves in the same 
way as the set of (2n+3) structures for ion I, the 
neutral dye molecule may be considered as a 
limiting case of a highly unsymmetrical ion in 
which one of the ground structures, IV,, is of 
very much lower energy than the intermediate 
structures, but the other ground structure IV, is 
just a little lower in energy than the intermediate 
structures, Wop41). This case is handled, then, 
by treating R as small and R’ as large compared 
to unity. The results have already been discussed 
for the special case that RR’=1. Equation 70 of 
reference 5 gives the transition energy for Case 


IV as: 


1 13 e 
AEy= (x+—) — 2 cos —— 
a R’ N+1 





\~— _, fo 
‘wniedatiiadenige 
N+2 N+2 


(1—R)? 
—|(RR’-1). (70) 
R’-R 





— R2N+2 


As N increases indefinitely, AEy for the neutral 
molecule approaches a finite limit just as in the 
case of a highly unsymmetrical ion. 

In conclusion the author would like to express 
his appreciation for numerous discussions with 
Professor K. F. Herzfeld and Dr. L. G. S. 
Brooker, and also for the courtesies shown him 
by the Eastman Kodak Laboratories. 
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Transformations of the Dynamical Variables 
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Problems in quantum mechanics may be solved by canonical change of the representation 
in terms of which the dynamical variables are expressed. They may also be solved by contact 
transformations of the dynamical variables. The usual notation is modified slightly to make it 
more convenient for developing these two methods. Their parallel development shows their 


essential equivalence but formal difference. 





INTRODUCTION 


HE motion of a molecule is complicated and 
requires for its description a number of 
dynamical variables which may be chosen in 
various ways.! If we choose variables, such as the 
normal coordinates for small oscillations, which 
vary with the time nearly independently, a zero- 
order approximation to the motion is obtained 
by neglecting entirely the interaction between 
these variables.? If the electric moment and 
polarizability of the molecule are expressed in 
terms of them, the absorbtion, emission, and 
scattering of light may be calculated.’ The gross 
structure of the spectrum of the molecule can be 
explained in this approximation, but to explain 
the details, the interactions must be taken into 
account. This is conveniently done by perturba- 
tion theory or successive approximation, going to 
the first order, second order, or even to higher 
orders, in the small interactions.‘ This more 
accurate work is necessary also if we wish to 
obtain the forces between the atomic nuclei in 
the molecule from the observed spectrum with 
even moderate precision.5 


1H. Kronig, Band Spectra and Molecular Structure 
(Cambridge, 1930); H. Sponer, Molekulspectren II 
(Springer, 1936); G. Herzberg, Molecular Spectra and 
Molecular Structure I, (Prentice Hall, 1939); D. M. 
Dennison, Rev. Mod. Phys. 3, 280 (1931); 12, 175 (1940). 
2H. Kronig, reference 1, p. 31; H. Sponer, reference 1, p. 
42. 
3H. Kronig, reference 1, Chapter 3; G. Herzberg, refer- 
ence 1, Chapter 3. 
4 oi Heisenberg, and Jordan, Zeits. f. Physik 35, 557 
1925). 
5D. M. Dennison, reference 1, p. 177. 


In the customary method perturbation theory 
is applied to the representation of the motion in 
terms of the stationary states of the zero-order 
approximation.* When there are many variables 
this becomes cumbersome, and is liable to error 
because of the large number of separate terms to 
be considered. We have found in a number of 
problems that the work is considerably abbrevi- 
ated by applying transformations directly to the 
normal coordinates; and this liability to error is 
reduced.’ While this is one of the earliest methods 
of quantum mechanics and is explained in many 
textbooks, it was found less useful than the 
customary method in dealing with atomic spectra 
and it has been little used.* The transformations 
are parallel to the contact transformations used 
in dynamical astronomy to remove terms from 
the Hamiltonian function. 


1. NOTATION 


The usual notations are a little clumsy for 
dealing with contact transformations, and we 
modify them slightly as follows: 

A symbol in square brackets as 


Ls] 


6 A. Ruark and H.C. Urey, Atoms, Molecules and Quanta 
(McGraw-Hill, 1930), p. 597; P. A. M. Dirac, Quantum 
Mechanics, second edition, (Oxford, 1934), Chapter 8; E. C. 
Kemble, Fundamental Principles of Quantum Mechanics 
(McGraw-Hill, 1937), Chapter 11. 

7 Shaffer, Nielsen, and Thomas, Phys. Rev. 56, 895 
(1939); W. H. Shaffer, J. Chem. Phys. 9, 607 (1941); 
S. Silver and W. H. Shaffer, J. Chem. Phys. 9, 599 (194i). 

8 A. Ruark and H. C. Urey, reference 6, p. 595; P. A. M. 
— reference 6, p. 110; E. C. Kemble, reference 6, 
p. 373. 


(1.01) 


532 





by 


ory 
. in 
der 
sles 
ror 
3 to 
- of 
Vi- 
the 
ris 
ods 
any 
the 
‘tra 
ons 
ised 
“om 


for 
we 


ianta 
ntum 
E.C. 


anics 


895 
941); 
941). 
\. M. 
ce 6, 








will stand for a matrix with elements 
(m|s|u). (1.02) 


In this symbol m and uw are variables or sets of 
variables designating the row and column of the 
matrix; they may be given by several variables 
taking integral values, or continuous ranges of 
values, or both. In general there need be no cor- 
relation implied between the values of m and the 
values of u, but most usually u is given by cor- 
responding variables to those giving m taking 
the same values, and we then conveniently 
denote them by m and n or by n and n’. 

The physical systems in which we are inter- 
ested are described in terms of dynamical vari- 
ables, &, », ---, which combine according to a 
non-commutative algebra we suppose known, 
given by kinematical conditions which usually 
comprise commutation relations: e.g., a coor- 
dinate g and its canonically conjugate momentum 
p satisfying 


gp — pqg=th (1.10) 


(h, as usual will stand for Planck’s constant, 
6.62 X10-*? erg-sec.; i stands for h/27). 

The dynamical variables can be represented 
accordingly in various ways by Hermitian ma- 
trices, 


[a|é\a], [a|n\a], seas (1.11) 
with elements 
(Qm|€|@n), (dm|9|Gn), ***- (1.12) 


In these symbols m and n are the variables or 
sets of variables designating the row and column 
of the matrix; a designatés the representation; 
and é, n, ---, designate the dynamical variable 
which the matrix represents. Since [a|£|a] is 
Hermitian, (am|{€|@n) and (a,|£{am) are con- 
jugate complex numbers. 

E.g., if qy/(2uw/h) and px/(2/uwh) are repre- 
sented by the matrices 


0 i 0 0 oO 
~— «© i/2 0 0 
~ivf2 0 in/3 0 

. 6¢. =f © t/ 


0 0 0 


~in/4 0 
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and 
0 1 0 0 0 
1 0 72 0 O 
0/2 0 v3 0 
00 V3 0 V4 
0 O 


0 /4 0 


and if we call the representation E, then 
LE|qvV/ (2uw/h)|E] and [E|py/(2/uwh)|E] are 
these matrices, and if we number the rows from 
the top by m=0, 1, 2, ---, and the columns from 
the left by n=0, 1, 2, ---, 


(Em|q/(2uw/h)| En) =i/(14m)  n=m+1 
m,n=0,1,2,---, 
=—in/(16m) m=n+1 
=0 otherwise. 


m may be given by several variables taking in- 
tegral values, or continuous ranges of values, or 
both; but ™ must be given bv corresponding 
variables taking the same values; and the sum- 
mation or integration giving the product of two 
matrices in the representation, which corresponds 
to the product of the dynamical variables in the 
algebra, must be specified or implied. We write 


Lalén|aj=[aléla}la|nla} (1.21) 


and 


(am | En | Qo) = Yon (Am|E|Gn)(An|n|@o), (1.22) 


with the understanding that >, is to be inter- 
preted as required. We abbreviate the unit 
matrix [a|1|a] to [a\a]. When there is no 
material ambiguity as to what representation is 
implied, (dm|£\an) may be abbreviated to 
(m|£|n), and since m and n take the same values 
it is often convenient to replace them by m and 
n', especially when they are given explicitly by 
sets of variables; e.g., (j, m|L|j’, m’). When 
convenient we designate the representation by 
the symbol of a dynamical variable or by those 
of a set of commuting dynamical variables, a, 
B, --:, say, represented by diagonal matrices for 
which the elements vanish when n#m. We may 
then use the diagonal elements, a’, a’, ---; 
B’, BY’, «++; +++; of these matrices, the charac- 
teristic values of the dynamical variables, to 
designate the rows and columns of the repre- 
sentation, and write 


La, B, ---|E|a, B, ---] 
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for the matrix representing ~ in this representa- 
tion, and 

(a’, B’, wee lela’, Pa ee -) 
for an arbitrary element of that matrix. The use 
of quantum numbers combines these notations, 
dropping the primes when there is no ambiguity. 
E.g., if g and are represented as above, and if 


E= ¥4(P?/ut+nw’¢”’) 
E/hw is represented by the matrix 


eco eé 


Mg 
0 
0 
0 


eek so 
oreo 


| 
. n\ 


and if we put E=(n+ 4)hw, we may write 


(n—1|qv/(2uw/h)|n) =ix/(n/2), n=0, 1, 2, --- 
(n| qv/(2pw/h) |n—1) = —ir/(n/2), 


other elements vanishing (though the phases of 
the elements are not determined by 


16 (b?/u+uw’g’) = (n+ 4)hw). 


If we have two representations, a and a*, for 
which the variables designating the rows and 
columns take the same values, and if dynamical 
variables £, 7, are represented in repesentation a 
by the matrices [a|é|a], [a|n|a], ---, we may 
define (different) dynamical variables £*, n*, ---, 
represented by the same matrices in the repre- 
sentation a*, so that 


(am* | &*|dn*) =(Am| €|an), 


(m* | n*|an*) = (dm | n | An) (1.32) 
or 
[a*| &|a*]=[a|é|a], 
(1.31) 


[a*| n*|a*]=[a|n\a] 


The dynamical variables &*, n*, ---, then satisfy 
kinematical conditions of the same form as those 
satisfied by é, n, ---. Conversely, if we have two 
sets of dynamical variables ~, 7, ---, and &*, 
n*, «++, satisfying kinematical conditions of the 
same form, and if £, 7, ---, are represented in a 
representation a by the matrices [a|é|a], 
[a|n|a], ---, we can define a representation a* 
in which £*, n*, ---, are represented by the same 
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matrices, so that 


(Am* | §*|an*) = (dm| E| an) 


(Am* | n* | n*) =(dm||n), (1.32) 
or 
La*| &*|a*]=[a|é|a], 
(1.31) 


La*| n*|a*]=[a|n\a] 


Such a symbol as F(é, 7, ---) will designate 
a function of the dynamical variables &, 7, ---; 
F(é, n, --+) will be itself a dynamical variable 
and may be defined in general by a matrix repre- 
senting it in some representation of the dynam- 
ical variables £, 7, ---, which matrix must be 
invariant for any change of representation that 
leaves the matrices representing the dynamical 
variables £, y, ---:, invariant; the most usual 
functions can be regarded as convergent power 
series in the dynamical variables, the order of 
factors in the terms being, however, important. 

If the dynamical variables £*, n*, ---, satisfy 
kinematical conditions of the same form as those 
satisfied by &, n, ---; and if F(é, 7, ---) denotes 
a function of the dynamical variables &, 7, ---; 
then F(é*, n*, ---) will denote the same function 


of the dynamical variables é&, n*, ---; e.g., if 
F(§, 9, ---)=a&+bn+-::-, 
F(&, n*, «+ +) =ak*+bn*+---. 
If we require the function of &, n*, ---, that is 


equal to F(é,7,---), we shall denote it by 
F*(&*, n*, ---) so that we shall have 


F*(, n*, + )=F(E, ny eee) 


and these may represent the same physical quan- 
tity, such as energy, electric moment, etc. E.g., 


(1.4) 


if g=q* cos e+ p* sin € 


p= —¢q* sin e+ p* cos ‘S40 = @-M-> 


and g*p* — p*q* =1h follow from each other), and 
if F(q, p) =a(¢—p*) +b(gp+pq), then 

F(q, p) = F*(q*, p*) = (acos 2e—5 sin 2€)(q?*+p™) 

+(a sin 2e+5 cos 2€)(q*p*+p*q"*). 

These notations are believed to be in current 


use ;’ except that it is not usual to write [a] £|@] 
for the matrix representing the dynamical vari- 


9P. A. M. Dirac, reference 6, Chapter 3, and p. 109; 
E. C. Kemble, reference 6, p. 268. 
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able ~ in the representation a, distinguished on 
the one hand from its elements and on the other 
from matrices representing ~ in other representa- 
tions; and that the idea of a function of dynamical 
variables is often confined to a function of a set 
of commuting dynamical variables. 


2. CHANGE OF REPRESENTATION AND 
TRANSFORMATION OF DYNAMICAL 
VARIABLES" 

We change from one representation a with 
rows and columns denoted by m, n, to another 
representation b with rows and columns denoted 
by u, v, by using transformation matrices [6|a ], 
[a|b], such that for any dynamical variable £, 


[b|§|6]=(b|a)La|é|a]La|b], 
[a| &\a]=[a| 6 J[b| €|bJ[b| a). 


[b|a) and [a|6] must be such that (5,|a,) and 
(am|b.) are complex conjugates and that 


[b|\a}La|b]=[b/b), 
[a|6j[b|a}=[a\a]. 


In general there need be no correlation implied 
between the values of m or n and the values of 
u or v. If a one-to-one correspondence is given 
between these values so that the methods of 
summation or integration correspond, we may 
conveniently denote them by the same symbol 
m or n, and denote the representations by a and 
a*; and we shall call the change of representation 
a canonical change of representation. (This 
specialization is not usually made explicitly.) The 
transformation matrix is then unitary and can be 
expressed in the form 


[a|a* ]=eil*) (2.03) 


(2.01) 


(2.02) 


so that 
[a*|a]J=est 


where [s] is a Hermitian matrix with elements 
(m|s\n) such that (m|s|n) and (n|s|m) are 
complex conjugates. 

We shall regard 


[a*|é|a*]=e-"I[a|élajet} (2.04) 
or 

[a| ¢|a]=ertI[a*| g|ar eH 
as the standard form for a canonical change of 
representation. P 





*” P. A. M. Dirac, reference 6, Chapter 3; E. C. Kemble, 
reference 6, p. 394. 


We transform from one set of dynamical 
variables £, 7, ---, to another set a, B, ---, by 
defining each set as functions of the other. In 
general the new set will satisfy kinematical con- 
ditions of different forms from those satisfied by 
the old. If they satisfy kinematical conditions 
of the same form, we may conveniently denote 
them by &*, n*, ---; and if &, 7, -- - are represented 
by matrices in a representation a, then £*, n*, --- 
will be represented by the same matrices in 
another representation a*. The change from a to 
a* is canonical and 


[a| é|a}=[a|a*][a*| &*|a*]Ca*|a] 
“"s[ala*Ifalélelerlay. (2-5) 


Thus if U is a unitary function of & 7, ---, 
which has matrix elements in the representation 
a given by" 


(dm | U| an) =(Am|an*), 


so that 
La| U\a]=[a\a*], (2.06) 
we have 
*= UEU-, 
~ x ae (2.071) 


and indeed, for any function F(é, 7, ---), 
F(&, n*, ---)=UF(E, 9, ---)U-' (2.072) 
(contrast this with 
P(e", 9”, ---) = PCE, 9g, ->-) (1.4) 


by definition). (an|U-|am) and (a@m|U|a,) are 
complex conjugates and UU-'= U-'!U=1. Such 
a transformation from dynamical variables £, 
n, -**, to dynamical variables £*, n*, ---, satis- 
fying the same kinematical conditions is called 
a contact transformation of the dynamical 
variables. 

U being unitary, it can be expressed in the form 


U=eissh, (2.08) 


where S is a real function of the dynamical 
variables (i.e., the matrices representing it are 
Hermitian). If 

[a\a* ]=e'l*), (2.03) 
then 

(dm|S| an) =h(m|s|n) 

or 

[a|S\|a]=AL[s]. (2.09) 

11 P, A. M. Dirac, reference 6, p. 109. 








We shall write 
S= S(é, 7, ** +) = S*(é, n*, iit +) (2.10) 
and shall regard the form 
gt = eiSlhte-iSih — * = eiSlhye-iSIh| ... (2.11) 


or 
t= eitSlh Ex giSlh n= e~tSlhyteiSih 


J 


as the standard form for a contact transforma- 
tion. A function of the original dynamical vari- 
ables é, n, ---, is then transformed as follows: 


F*(é, n*, sie -)=F(é, °° -) 
= e-iSINF(E*, gt, «+ -Je'SM, (2.12) 


E.g., if go—pq=th, S=cp gives 


g* = ercplhge—teplh = qt+c, 
p* ~— etc Plh pe-icrlh = p 


and F*(p*, g*) = F(p, g) = F(p*, g*—0).” 

[a*| F\a*]=e-‘lI[a| Fla jes! (from (2.04)) 
and 
F*(e*, nt, +++) =e-SIMF(E*, og, + -)eiS/@ (2.12) 


are actually just different ways of writing the 
same result. The formula 


1 
e- (WK Fee ih) K = F+-(FK—KF) 


+:(-) (#k-KAK-KEK-KF)+ ‘an 
= FALK, FI+4CK, CK, FI}+--- 
$00K 00 0K, FD (2.13) 
where the Poisson bracket [K, F’] is defined by 
[K, F]= ~“(KF-FR), (2.14) 


is useful for evaluating these expressions. If 
F(é, n, ---) and S*(&€, n*, ---) are known as 





[a|a* |=e*!*) to (from (2.03)) 


[a|H\a]=[a|Ho\a]+ela| Ai \a)]+e[a|H2\a]+---, (from (3.2)) 


2P. A. M. Dirac, reference 6, p. 94. 
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13 P. A. M. Dirac, reference 6, Chapter 8; E.C. Kemble, reference 6, p. 395; A. Ruark and H. C. Urey, reference 6, p. 596. 















power series, this formula gives 


F°(g*, 9°, ---) = P(E, 9°, ---) 
+([S(é, n*, i -), F(&, vw. wie ‘)J 
+ 4LS(é€, n*, ro ), LS(é, n*, ia *), 
xX F(é, n*, oo) }+-*-. (2.15) 


3. THE MECHANICAL PROBLEM. PERTUR- 
BATION THEORY 

We suppose that the Hamiltonian function HH 
of a system, and the various physical properties 
in which we are interested, such as electric 
moment e and polarizability , are given in 
terms of dynamical variables the kinematical 
relations between which are known. Our object 
is to find matrices representing H, e, ®, etc., in 
a representation 6 in which H has only diagonal 
elements: to find values for u, v, and for E,, such 
that 


(b.|H|b,)=E, v=u (3.1) 
=0 otherwise, 
and to find (b,/e|b.), (b.|®|b,), --- etc. We 


divide this into a zero-order problem, and per- 
turbation theory." 

The zero-order problem is to transform to 
dynamical variables £, , - - -, such as components 
of angular momentum and normal coordinates 
and momenta, in terms of which H has nearly a 
form Hp» for which we know such a representa- 
tion, a. H is conveniently expressed as a power 
series in € 


H(é, n; ---)=Ho(é, Up -++)+eHH,(&, Ny ++) 
+ €H2(§, n, see)pees, (3.2) 


where ¢ is small (e.g., e=h/vJ where » is a vibra- 
tion frequency, J a moment of inertia), and 
where values for m, n, and Eo, are known, so 


(Am |Ho|\an)=Eom m=n (3.3) 
=0 otherwise. 





We must also be able to find the matrix com- 
ponents in the representation a of functions of 
&, », ---, in particular (@m|Hi|an), (@m|H2|a2), 
-++, (Gm|@|Gn), (dm|®|a,), ---. If we apply the 
canonical change of representation for which 
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[a* | H|a* ]=[a*| Ho*|a*]+La*| H,*|a*]+€[a*| H.*\a*]+--- 


we obtain, expanding in powers of e," 


[a*| Ho*|a* }=[a|Ho|a], 


[a*|H,*|a* ]=[a|Hi\a]+<(La|Mo\a js ]—[s ]La| Ho\a)), 


(3.41) 
(3.42) 


[a* | H2*|\a* |=(a| H2|a]+¢(La| Hi|a Ls ]—[s ]la| Hi\a])) 
+ 1?((La| Ho|a }Ls]—(s JLa| Ho|a])(s]—[s](Ca| Mo\a)[s]—[s]la| Ho|a))) 


=[a| H2|a]+ 13i(La|+H*\a)[s]—[s][a|i+H*\a)), 


If we apply the equivalent contact transformation from &, 7, 


to (from (2.08) 
H(é&, q, °° -)=Hd(E, ~** -)+eH7,(é, nN, ° 


transforming it to 


H*(&, n*, ae ) = Fo*(&*, n*, °© )+eF*(e, 7, 


we obtain, expanding in powers of e, 


Fo*(&, n*, a +) = H)(é, n*, ee -), 
F,*(&, w ii -) =A, (E, n*, mF ‘)+LS(é, n*, 


o -), Hy(&, n*, oe -)], 


(3.43) 


15 


reo to &*, , oie, given by U=eisslh 


»-)+€H2(E, 9, +++) +--+ (from (3.2)) 
+ )beF (Ee, nt, +s) pee ot 


(3.51) 
(3.52) 


Fy*(&, n*, met -)=H.(£, n*, we )+LS(E, n*, pias -), Ay (&, n*, 7 ‘)] 


+ 4 S(é, at a *), [S(é, n*, 
= H,(&, n*, “2 J+ ULS(, n*, ih *), 


These results we shall write 


H,* =H, 
H,y*=H,+([S, Hy), 


Ho*=Hot+ 16S, H,+H;*], 
H;*=H;+ (5S, Hy ]+ %LS, LS, 2H2+H* }], 


We remove terms from H by taking suitable 
forms for s or S. If the zero-order approximation 
is not degenerate, i.e., has no coincident energy 
values, we remove first terms of order e; then 
terms up to é, then to e’, and so on; only terms 
that are functions of Ho remain, i.e., that are 
diagonal in the representation a*, a**, ---, that 
we have reached. If the zero-order approximation 
is degenerate, and if H,, He, --- contain terms 
that remove the degeneracy, these terms also 
remain, and we must solve a zero-order per- 
turbation problem, changing the representation 


abo Heisenberg, and Jordon, Zeits. f. Physik 35, 565 
® Often applied: e.g., J. H. Van Vleck, Phys. Rev. 33, 
484 (1929); O. M. Jordahl, Phys. Rev. 43, 87 (1934). 
© Note that 
Hy*(é*, n*, sake *) = Hy(é*, n*, oes )+eLS(E*, n*, ai ), 
Hy(é*, n*, sah -)J+ sit cle # Fo*(é*, n*, iy -). 





+++), Ho(&*, n*, +) 7] 


H,(&, i sik -)+ Fi*(é*, n*, whe +). (3.53) 


(3.51) 
(3.52) 
(3.53) 
(3.54) 





or transforming the dynamical variables suitably, 
to approach the representation b. 

If this process converged we should approach 
the desired representation b, and could apply the 
same transformations to e and ®. In any case, if 
¢ is small, the first few transformations must give 
results with errors of order &, ¢, «8, ---, in the 
elements of e and @ as well as in E,,. 

Successive canonical changes of representation 
are especially useful when we want results for 
only a few values of m or n; contact trans- 
formations where the functions can be expanded 
in power series in the dynamical variables. We 
shall give tables of forms for S which give various 
values of [.S, Ho] for standard forms of Ho, and 
of the resulting higher order terms that are 
introduced into HZ. 
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The non-commutative algebra of polynomials in a coordinate and its conjugate momentum 
is reduced to common algebra by Weyl’s method, and tables are given for facilitating its use. 
It is shown how the problem of the anharmonic oscillator can be solved by contact trans- 
formations, and tables are given for removing terms up to the third degree in its coordinate and 
momentum and for finding the modified Hamiltonian up to the fourth degree, which is as far 
as is ordinarily required in molecular theory. To illustrate the power of the method, the energy 
is computed up to terms of the eighth degree for energy containing terms to that degree (in 
the coordinate only), obtaining Dunham's result together with the constant terms of that 


order not given by him. 





4. POWER SERIES IN A COORDINATE AND 
ITS CONJUGATE MOMENTUM 
[‘ 


gb— pq=th 
and if we deal only with terms up to the fourth 
degree in these variables, we may apply non- 
commutative algebra directly to expressions such 
as pgp’ and take explicit account of relations such 
as p'g= pqp*—2ihp*. It is, however, convenient 
even here, and almost essential for extension to 
terms of higher degree, to adopt a definite con- 
vention for a standard combination of such terms 
as p*g, p’qp, pap’, and qgp*, and to construct a 
general theorem for combining these combina- 


tions. 
We define {p’g*} by the expansion 


enritem 2. a pit’s = 


ris! 


(4.1) 





(4.2) 


where a and @ are ordinary numbers, so that, for 
instance, 


(P’?} =“ (p’e+papqt+pept+aqpap+ Pr’? + 9b’) 


and {p’q*} is r!s!/(r+s)! times the sum of the 
(r+s)!/r!s! terms distinct in non-commutative 
algebra which become ’q‘ in ordinary algebra.! 

1 This definition follows H. Weyl, Zeits. f. Physik 46, 1 
(1927). The definition given by M. Born and P. Jordan, 


Zeits. f. Physik 34, 858 (1925) and expanded by Born, 
Jordan, and Heisenberg, Zeits. f. Physik 35, 557 (1926), 


We shall write 
om a,sp'q’} => ars{p'g'} 


for a linear combination of terms of non-com- 
mutative algebra constructed in this way. With 
the aid of 4.1 any polynomial of non-commuta- 
tive algebra can be expressed in this form in one 
and only one way (Table I). E.g., 


p*q= |p'g— 3gthp*} 
= 6(p'q+P’ap+t par’ +4p*) — 3sthp’. 


TABLE I. h=1, [g, p]=1. 











1={1} ={q'} 
p= |g p+ 3oi¢"} 
g={q} qg= {|g pt eigq?} 
p=\P} in| \oig*} 
pg =\gp—3oig*} 
g={7} Pp =\¢p?+2igp—}o} 
gp =\qp+ li} gpqp = gp? +igp} 
pq =\qp— 41} gP'g=\9p?+1)o} 
P={p"} PP p=\|PP?+'4} 
pqpa=\¢p?—igp} 
g={g} Pe =| ¢p?—2igp— 4} 
Pp={Ppt+ig} gp* = | gp? + 391p} 
gpq=\¢p} - pap? =| qp>+ doip*} 
Pg =| ¢p—iq} P'qp = | gp® — eip*} 
gp? = \qp*+ip} P'q= | gp*— 3eip*} 
bgp =| qp*} pi={p'} 
P’'a=\qp?—ip} 
p={p*} 








Reintroduce & in terms of degree 1/gp of the leading term, #? in terms 
of degree 1/g*p?, and so on. 


leads to considerably more complicated results. It should 
be noted that we use differentiation only of the polynomials 
of common algebra. 
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TABLE IT. h=1, [g, p]=1, [{A}, (Bi J=1C. 














ke 

A NS q ? g qp p* @ gp qp* p§ 
q 0 1 0 q 2p 0 a 2gp 3p? 
p =i 0 —2q eg 0 —3¢ —2qp wae 0 
¢ 0 2q 0 2¢ 4qp 0 2¢° 4¢°p bgp? 
qp | —49 p —2¢ 0 2p? — 3q° —¢p qP* 3p 
P | —2p 0 —4gp —2p 0 —6¢*p —4qp* — 27° 0 
g 0 37° 0 3q° 6q°p 0 3q° 6q°p 9¢°p*?— 39 
gp | -¢ 2qp —2¢° gp 4qp* —3q' 0 3¢p*?+ 4 6gp* 
gp? | —2gp P —4¢p —gp 2p' —6q°p —3q¢p?— 15 0 3p" 
p* —3p* 0 —6¢P —_ — 3p 0 — 9p? + 39 — 6gp* —3p' 0 











Reintroduce h? in terms of degree 1/q?p? of the leading term, #4 in terms of degree 1/q‘p‘, and so on. 


We make use of the lemma? 
errtBa = eareBighhab = ebigare—thaB (4.31) 


or rather its corollary 


exPtBagypt+bq — e(aty) Pe (B+5) aehih(By—ad) (4.32) 


to obtain the result, fundamental for working 
with polynomials in p and g, (D,(A) means 0/dq 
acting on A but not on B, and so on), 

{A} { B} = {elit als )Dp(B)—Dp(A)Dq(B)) 4 B} 


0A 0B 0A OB 
={AB\+1 si oprah | 
0 


0°A 0°B 


dg” dp* 
A @B @2A 0°B| 


+ — 
dpdg 8pdq ap? aq?) 
0A 0°B 
+ 36(34i0)5| —— 
dq Op? 
0A 0°B 0A a°B 
3 +3 
0q’0p Agdp? = Agdp* dg?dp 


0A 0B 
————}4..-, (44) 
dp* dq’ 


In particular the Poisson bracket 





+ 14( 14th)? 


—2 

















0A 0B 0A OB 
[1A}, (B}}= | —— 
0g Op Op aq 
0°A 0°B A 0B 
+16(34ih)®  — -— 
0q*® Op*® dg*dp dgdp? 
7A 8B 0°A 0B 
+3 we =| (4.5) 
Ogdp* dq’dp dap* dq 





* H. Weyl, reference 1; the lemma follows from the proper- 
ties of the group. 


This last result connects quantum theory contact 
transformations with those of classical theory, 
enabling correction terms to be written down 


(Table IT). 


5. THE HARMONIC OSCILLATOR—EXCURSUS 
ON NOTATION 





1 mw? 
If H=—p?+ q’, (5.1) 
2m 2 
where 
gp — pq=th (4.1) 


the usual reduction of H to a diagonal matrix 
gives? 
(n|H/kw\|n)=n+1¢ n=0,1,2,--- (5.21) 


the remaining components vanishing, 


1 
(»-+1) —o/ vim) -+ivima n 





= /(2(n+1)) n=0,1,2,--- (5.22) 


the remaining components vanishing, 


n+1) 


= /(2(n+1)) n=0,1,2,--- (5.23) 


the remaining components vanishing. 


1 
(»|—cw/vime) —ty/ (mw)q) 


IV 








3 W. Heisenberg, Zeits. f. Physik 33, 879 (1925); M. Born 
and P. Jordan, Zeits. f. Physik 34, 858 (1925). 

The derivation is:—if H= 4p?+ 4q?, gpb—pq=1; write 
ptiq=a, p—iq=8, and 4H =a8+ 8a, aB—Ba=2; Ha—aH 
reduces to a; hence either (”|a|m)=0 or H,, —He =1;and 
(n|a|m)(m|B|n )= | (nx m) P= 2H. —1=2H,»+1. The 
series «++, Hn, Hn—-1, Hm—2, +++ must terminate since it 
must be positive, and, can do so only for the value 4%; the 
values of H are ---, 54, 39,°%; the representation can be 
chosen so that (n+1|a|m) and (n|8|n+1) are real and 
positive, and the result follows. 
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These variables and matrices have been used for 
the examples of Part I of this paper. 

Here we have used m for the mass of a particle 
moving under a restoring force proportional to 
the displacement such as to give a classical fre- 
quency of oscillation »=w/27. This is customary 
in theoretical work. The use of # and w removes 
most 27’s from the formulae.* 

In applications to diatomic molecules it is usual 
to use uw for the reduced mass in place of m and 
to write —kq for the restoring force so that we 
write 


1 1 


H=—p*+-—kq’. (5.11) 
"Sg 2 , 


It is usual to reserve v for the actual wave number 
of a line (in cm) rather than for the frequency 
(in rev./sec.), and to reserve w for a value (in 
cm~'!) corresponding to the classical frequency 
for small oscillations. If they occur in the same 
context we shall write voce and Wose for the 
former, frequency of oscillation in rev./sec. and 
radians/sec., and 7 and @ for the latter, observed 
wave number and corrected wave number, 
respectively. 





(n+r|(p+ig)"(p—1g)*|n+s) = 


the remaining components vanishing. From 


we obtain 
(n+r| {(p+ig)"(p—ig)*} |n+s) 
V((n+r)!(n+s)!) ( 
= Qr+s/2 1 





n! 


_V((n+r)(n+s) J 
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For polyatomic molecules it is convenient to 
use dimensionless coordinates, s= \/(Mwose/h)q 
with the corresponding momenta 


p= V (h/moose)P 
so that 


1 1 ah 
eT ee ae ees | 
c 


for the harmonic approximation for a single 
normal coordinate.® 

For the purposes of this paper it is convenient 
to choose the unit of mass so that m=1, that of 
time so that w=1, and that of length so that 
h=1, (this has already been done in Tables I 
and II). Thus we work with 


H= lop? + 1sq", 
gb — pq=t. 


To return to the usual units, replace p by 
b/V/(mhw), g by g\/(mw/h), H by H/he. 

We can obtain the matrix components of 
{p“q’} in the representation given by 5.21, 5.22, 
5.23, as follows: From 5.22, 5.23, 


(5.13) 
(5.14) 


Vint) nts) 








+s /2 

n! 
n, r, S, positive integers or zero, (5.3) 
et (ptiaq)+B(p—ia) = eaBea(ptia) eB (p—ia) (4.31) 

‘Ss 1 r(r—1)s(s—1) 1\? 
a ) 
1-(n+1)\2 1-2-(n+1)(n+2)\2 

rtei3F(—7, —3; 8+; 34) (5.4) 


n! 


n+r,n+s20, the remaining components vanishing (if is negative, the first —” terms of the series 


4P. A. M. Dirac, Quantum Mechanics, second edition (Oxford, 1934), pp. 90, 133; H. Weyl, The Theory of Groups and 


Quantum Mechanics (Methuen, 1931), p. 51. 


5G. Herzberg, Molecular Spectra and Molecular Structure | (Prentice Hall, 1939), Chapter 3. 


6S. Silver and W. H. Shaffer, J. Chem. Phys. 9, 599 (1941). 








it 
yf 
it 


ies 


ind 
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vanish); which, since 











ut+v (p+1igq)"t’—*(p —ig)* uly!(—1) 
wort | 
s=0 ante b b'\(r—s+b)!(s—b)!(v—))! 
oe ( +1 )sto-o/( — )s ! 
eg sini tM tin. | 2, nests <2) (5.5) 
s=0 2ete s\(r—s)! 





enables us to write down the matrix components of | p"q"}, e.g., returning to the usual units, 


hn 5/2 
(n-+4| [p%a4|In-+1) = V(mw)(—) }o(n+3)V/((n+4)(n+3)(n+2)), (n+2)20. 


6. THE ANHARMONIC OSCILLATOR 


If H=\6p?+4¢?+eHi(p, gq), (6.1) 
where 


gpb— pq=1 (5.14) 


we may remove a term c{p’q*} from 17, by a contact transformation to p*, g*, with a function ec{.S} 
given by 


[1op*?+ 1sq*?, |S} ]={p*q**} (6.21) 

i.e., (4.5) “ - 
f-—-f—s fy" (6.22) 

op* ag* 


so that the equation to find S is for the anharmonic oscillator of exactly the same form as in classical 
dynamics. (In returning to the usual units, replace S by S/h). We then have (3.53) 


H* =H = \op**+ log? +e(Ai(p*, gt) + [ce S}, op? + eq" )) +eLe{S}, Ai(p*, g*) 
é 
+4[c{S}, pe+ req JI+—> Let Sh, Cco{S!, Hi(p*, g*)+15Lel S}, sp**+ log" ]]]+---. (6.3) 


To save writing it is convenient to drop the asterisks, remembering that both {.S} and the trans- 
formed H are in terms of the transformed variables. 
Thus for s odd we may take 











prtig- — (s—1) (s—1)-+-2 
area ermemeane- fle Mie ili 2 —or |S (6.41) 
(r+1) (r+1)(r+3) (r+1)(r+3)---(r+s) 
and for r odd 
r—1,y8+1 —1| —| is -2 
if. Oe iptv (r—1) — (6.42) 
(st+1) (s+1)(s+3) (s+1)(s+3)-+-(s+r) 


If ry and s are both odd, these two expressions differ by a multiple of (p?+°)’**/?, and the sum of 
either and any multiple of this may be taken. (The transformation function reducing H to a function 
of 14p*+ 14q9** is itself in fact only determined up to an additive function of 14p**+ 14q*, though 
the reduced form of H is completely determined.) 

If y and s are both even, c{p"g*} cannot be removed, but can be reduced to C{(p?+q)"**/*}. If ~ 
means can be reduced to, , 
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TABLE III. 

{q?} w 46 {(p?+49")} =4(P’+¢@) 

(pj 1ot(p?+¢")} =lo(p?+¢") 

{gt} ~ 3¢ {(p?+q")?} ms} AY a 
tgp} wig {(p2+q*)*| = 16 ((p24+-g*)?+ 1) 

{pt} w3¢{(p?+q")*} =3¢((p?+ 97)? +1) 

[gin 6 l(PP+q)3} = 6((p? +9")? +5(p?+¢")) 
gip?} sw 46 | (p?+¢")5} =i i 6 ((p+q2)*+5(p?-+9°)) 

{q*p*} w 46 | (p?+4")?} = Yel (p+ q*+5(p 2+-@?)) 

(p8}m 56 l(p?+q")8} = 6((b? +97) +5(p?+¢")) 

{q8} wm 35408 { (p?+-9")*} = 35408 ((p? +97)! + 14(p?+¢7)° +9) 
{g° at Sd = Ko8((p?+¢q?)'+ 14(p?+¢7)? +9) 
{q*p* pre | Coad aye = 54 oe (pttat)'-+-14(pe-+-g2)2-49) 
{g°p®} mw 5408 {(p?+q")4} = 54 08((p?+9")°+ 14(p? +9")? +9) 

{p§} 3508 (ete) = 35498 ((p?+q")'+ 14(p?+9")?+9) 








Replace p by p/v (mhw), g by qv (mw/h). 


TABLE IV. h=1, [g, pJ=1, [44(p?+@), {S}J= ieunad 














Ai § Ai S 7 

ine . g—36(P +) 5gq°p-+ 36ap* 
gp — 550q'+ 3 167°?’ + 34 2p" 

: P er KiP+eyY — — ea'pt Hea 
p mg =. s(pesotye so2f a gal 3 op* 
q 2—16 3 (p?+4q?) ap P q 36 p— qp* ; 
3 14 (p?-+q?) —— ‘. Upt $39p> + sP* 

: ; ap? satpt +2) sp 
- Cpt ?6p op? 21 5¢°— Yeah 
gp — hq qp' 1sp 
ap* BP p —8isa°— $60°P*— 9p" 

—34¢°—g 











Replace p by p/ ¥ (mhw), q by gv (mw/h), H by H/hw, S by S/h. 














2n—1 (2n—1)(2n— 5), 
ign} {q?"-2p?} mw 7 g?"— tps }n- -om{p}el [(g?+p?)"} 
n 1 n(n—1) 1-3 1-3---(2n—1) ; 
i + $e tt]={@+p"} (6.5) 
1(2n—1) 1:2. (2n—1)(2n—3) 2:4-+-2n 


(Tables III and IV). 

For the anharmonic oscillator the difference from classical mechanics enters in the Poisson 
brackets, (e.g., in the usual units, [{q*p}, {gp*} ]= {15¢*p*— 6hq?p? — 3¢h?} of which the first term 
only would occur in classical mechanics). 

When A, can be expressed as a power series in p and q, a series for S is obtained reducing H toa 
power series in {(p?+q?)"} together with terms of higher order than e; a second transformation 
reduces H to a power series in {(p?+q?)"} together with terms of higher order than ¢é, and so on. 

(n|H|n) can then be found from 5.3, or reduction formulae giving {(p?+g")"} in terms of (p?+4q”) 
can be built up (Table III). If terms beyond é& can be neglected and if the terms of order ¢ in // 
are of order three or lower in p and q, (n|H|m) can be written down immediately from Table V. 


7. THE SLIGHTLY ANHARMONIC OSCILLATOR 


1 1 
If H= (<r*+ mere’) +ecsq?+ecagit+-:-, (7.1) 
2m 2 








50 
rm 


oa 
ion 


on. 


1 H 


7.1) 
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TABLE V. Hi=Hy'+Hy", rte [S'|]=Hi— Hh'*, [6(p+@), Hh'*]=0, 
MOU S"}, Mi" +HM'"*)= 

















































































































Sy” 
B {q) {p} {q?} {ap} {p?} {q"} {7p} {1p"} {9} 
HN - = 7 = os aaa io Es 
{7} {-—}9} 0 {-1eq} {—}4p} {—}4q} {—99?} {—ap} {—!¢p*} 0 
—l4 0 0 0 0 —34 (p?+q?) 0 —14 (p?+q?) 0 
{p} 0 {-%) {-¢p} {—}4q)} {-—} 9 0 {—!49"} {—gp} _{>3ar"} 
0 —12 0 0 0 —4P +e) — 0 4 (P+ ) 
te) | (Sal | tah tsa 3en% 0 {3402+ den?) (fae) Ii 49°?) {'4qp?} (49°) 
0 0 0 0 | 0 | 
a a, a eee ——— —-|-—--___|-- .___. 
{ap} | {—14p} {-14q} {—eqp} |{—4a?— M492) {-} (ar) {-34¢@ 2p) [-Ke- Lap?) | { -aq*p—!4p'} {—34qp?} 
0 0 0 —V4 (p2+q2) 0 0 } 0 | 0 
{ p?} {'4a} {-l4p} {'oq?—}4p?} 0 {!eq?—ep?} {493} {'4q*p} | a - Yan} { -34p4} 
0 0 0 0 0 0 ) 0 0 
{P} |{—3¢q?—p?} {ap} ce” 2q°} {¢q*p—p*} {—40°)  |{—34q'—3q2p?+14} | { —2qp*} | {s92p?+14—ps} (3p? 
~3;(p2-+q2) 0 0 0 —154 6 (p2+q2)2—7 0 ms 84 6 (pt-+q2)2-+54 6 
{q°p} 0 (—'40") | {—1492p) {—1493} {—4q*p} | 0 { -14q {-¢°p} {—34q¢p?+4} 
0 —14 (p?+q?) 0 0 0 | 0 -ovete 7 fo} 0 —346(p?+¢q2)?+Ke 
{ap} | {—1p"} 0 ty qp?} {—!4p)} | {—!ap*) | {—3¢q%?+14) | {—¢p"} | {—% 0 
—4(p?-+¢") 0 0 0 —{h¢ (p* Feet 0 ‘abd Set P HP) ts 346 0 
{*) | fap} fa? 3 ap*})  {—1ep3} | {-a34+2¢qp"} | {—1ep3} | { ‘a Cc = 9+ 3592p +14) { —2q%p} | —3a*p?-+ 9g — 9 apt} 
0 |=84(pt-+q2) 0 ™ » @ | ‘olm+aPPt Mel 0 —346(p?-+42)?—34 6 
| | 
Replace p we p/¥ (enh, q by ¢/vV (mw/h), H by H/hw, S by s/h 
where 
gp — pq=th (4.1) 


transformation with 
p*q** 2 gps 
S= ca —— +4+-— —— (Table IV) (7.2) 
mw? 3 mo! 
vives to order é 











1 1 3 C3" 3¢3" h*c3? ; 
H*= | Pt matgt + e((« or pegt+ )+ tee (Table V) (7.3) 
m 


2 mw? mg)4 2m*w! 


which a further transformation reduces to 


a ( 1 1 mw 
** — fy pr24— —_ me) 
2mhw 2h 


3 C4 15 C3” 1 or 
+e((- ie )((—pre:+ maar) +8) +— ~~) (Tables III, V) (7.4) 
8 m*w? °16 mw! mw 2m>, 


giving for the nth state, since 


























1 1 mw m 4 
p**?+— —¢ *2—y 14, 
Qmho 2h ; 
eh? C32 £30n?+30n+11 . 
=(n+}; 36¢4(n?-+-n-+ 14) —- ( ))+ vee (7.5) 
mw? mw? 8 . 
Further, if 
e=e)+eig+€*e2q°+ - ated (7.6) 


transformation with 7.2 gives 





” p*? 
=er(pt, at) ert eg" te( engt—csei(——+2 ))t (Table I1). (7.7) 


mw? mw! 
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To this order the further transformation makes no change, and the matrix components of e in a 
system in which // is diagonal can be written down at once.’ 


8. DUNHAM’S RESULTS 


To obtain the terms of the next two orders in E for the anharmonic oscillator, Dunham’s? results, 
Table II (formula 4.5) were extended to cover [{A}, {B}]={C} up to terms of order q* in C. If 


H= Wyp?+ 46g’ +aq? +bgq'+cq>+dq'+eq' +far+::-, (8.01) 
where 
qb— pq=1, (8.02) 
transformation with 
S=a\ pq?+ 25p*} (Table IV) (8.03) 
removes aq giving by 3.51, 3.52, --- and using Table II repeatedly, 


H= | 34(¢?+p*) +bg'+ 13a*(— 3q*— 6q*p? + 1) +q°+ab( — 4° — 8q*p? + 4q) 
+ ¥4ga8(129°+ 129%p? — 18q¢+ 24gp*) +dq°+ac( — 5q°— 109*p? + 10g") + }3a°b(209°+ 32g*p? 
— 56q’+489*p* — 16p*) + 1ga*( — 609° — 108q*p? + 1569? + 96q*p* + 120p? — 48p°) + eq’ 
+ad(—6q’ — 129°? + 20g*) + }a°c(30g’ + 609°p? — 140° + 809*p* — 80gp") 
+ 1¢a°b( — 120g’ — 2409°p? + 576° + 329%p* + 576gp? — 192gp*) + 130a°(3609' + 7209°p? 
— 1728q?+ 1440q*p* — 960gp*) +fq?+ae( — 7q°— 149°p?+35q') + 1ga7d(42qg'+ 96q*p? 
— 300g!+ 120q*p* — 240q*p?+ 20) + 1éa*c( — 210g* — 480g*p? + 16204 — 200g*p4 + 192097p? 
— 180 —480g2p> + 240p*) + 15404b(840g8+ 19209%p? — 6528q!+ 2560g'p! — 3840q°p? + 864 
— 2304g*p* — 2560p1+ 384p°) + 14440°( — 2520g* — 5760g*p? + 19584g+ 30528q"?p? — 1728 
—1920092p°— 5760p'+1920p*)+---}. (8.04) 
A second transformation with 
S= {b(96q%p+ 36ap*) + 140°(— 26a" —136qp") +e(a'P + $50°D*+ 515") 
+-ab(—4q'p—8q%p3— 195p5-+4p) + 14a%(12g'p+2092p?+%45p°—18p)} (Table IV) (8.05) 
reduces 
{bg*+ }4a*(—3gq'—6g*p?+ 1) +c9g>+ab( — 4q°— 89*p?+4q) + 140°(12g°+ 129*p* — 18q+ 24gp*) } 


to 


{ 36b(q?-+*)?+ 1 60*(—15(g?+p*)?-+8) |, 
giving, by 3.51, 3.52, ---, and using Table II repeatedly, 
H= {36(g-+p*) + 360g? +P")?+Hiea*( — 15(g?-+p*)* +8) +dg*-+ac( —5q°— 10g*p?+ 10g?) 
+ 14207b(4679°+713q*p? — 12749¢?+ 7772p — 58p? — 45 p®) + i 28a*( — 131196 —2781q*p* 
+3900g?+7719*p'+ 840p? — 543p*) + 142b*( —55q°— 69q*p? + 72¢° + 27q°p* — 36p? + 9p) +eq' 
+ad(— 6g’ — 129°p?+ 20g*) +a°c142(6819’ + 13419g°p? — 32459°+ 16649*p4 — 1856gp? 
+ 256gp*) +4%b!54(— 1497g’ —416195p? + 82479 — 3296g*p!+9210gp? — 2028¢p*) 


7M. Born and P. Jordan, Zeits. f. Physik 34, 858 (1925). 
8 J. L. Dunham, Phys. Rev. 41, 713, 721 (1932). 











)4) 


)5) 
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+a°!9(2970q7 +7020g*p? — 1831893+ 11430g3p! — 19890gp?2+ 3460gp*) 
+be'48(— 207g’ —555q°p? +507q' — 256g*p4 + 768gp*) +. ab? 4 (69g7 + 2439°p? — 307g" 
+1949*p* — 396gp") +fq> +ae(— 7q°— 149°p? + 35q') + 14bd(— 15g°— 279°p?+45q') 
+ 1¢c2(— 159° — 60g*p? + 609! — 4094p! + 240g2p2 — 12) + 1ga2d (1959 +519g%p? — 1425! 
+480q'p!— 960g2p?+80) + 4 sabc(465g8+ 12159%p? — 30519°+1570g'p* — 385892p2 
+ 288+ 384q2p* — 4484) + 1425b9(4259q°+ 620g%p? — 2310¢!+930g'p'+612g2p2+ 1089p! 
— 342p'+27p*) + '2sa%c( — 1725q% — 5595q%p? + 139599! — 3850g*p' + 16794q?p* — 1488 
— 1664g2p*+ 2688p! — 256p*) + 1456a2b?( — 2312598 — 601169? + 185450! — 65150g'p4 
+ 2080740q?2p? — 20480 — 55380g3p°+45242p — 405°) + 147204b(5866149%+ 19271289%p? 
— 5241036g'+ 1897540g*p4 — 6127800q2p? + 700416 + 246744g2p* — 1905004p! + 216438p') 
+1{924a°(—89903g* — 339260g*p?+ 841278q' — 259930g'p!+ 1210284922 — 131072 

— 59324q°p°+358574p!—42543p5)+---}. (8.06) 


A further transformation removes the terms of the order of eg’ and gq’ and reduces H to the form 
(Table ITT) 
H= 35(q2-+p*) + 16b(3(q2+P°)2+3) +H 60*(—15(g2+ f°)? 7) + Hod 5 (G+ 2)*-+25(g2+p2)) 

+ eac(—35(q?+ p?)* —95(q?+ p*)) + 1420°b(225 (g? + p?)? +459(g' + p?)) 

+M28a*(—705(g?+ p*)* — 1155(g?+ p*)) + 1526?(— 17(g?+ p*)* — 67 (9° + p?)) 

+ Mos f(35(g?+p2)!-+490(g?-+p?)? +315) + 2sae( —315(g?-+p")*—2730(g?-+p*)?— 1155) 

+ osbd(—165(g2+p?)*—1770(g?+p*)®— 945) +1456¢2( —315(g?+p2)4—2170(g2+p)? 

— 1107) +1456a7d(2715(g?+ p?)*+17070(qg?+ p?)? +6055) + 4 2sabc(2415(g?+ p?)4 

+14670(q?+ p*)? +5667) + 1256b*(6000(q?+ p*)*+ 13656(g?+ p?)?+ 1539) 

+ '456a*c( —9765(qg?+ p?)* —47730(q?+ p?)? — 14777) +14 1207b?( — 24945 (g?+ p?)4 

— 248052(g?+ p2)?— 40261) +! 02404b(116325(g2+ p2)'+479970(g2+ p2)2+ 131817) 

+ 4o48a°( — 115755(qg?+ p?)*—418110(g?+ p*)?— 101479) +terms of order a’. (8.07) 

This result agrees with that of Dunham by the W. K. B. method;* the terms Mosf 315 

— Mes ae 1155 — Mog bd 945 —--- — 14048 a8 101479, not given by Dunham, check with an inde- 
pendent calculation by Dunham’s method. The method used here may be regarded as differing from 
that of Dunham essentially only by expanding the potential energy in a power series first, and then 
expanding the transformation in powers of fh, while the W. K. B. method expands in powers of hf 
first and Dunham expands the potential energy later. Except for the terms not given by Dunham, 
his method obtains the form for H more briefly than this method; but these terms, the only ones of 
order h?, require in his method nearly as much time as all the rest; in this method they are obtained 
with the others. We see also that, while for perturbation theory we require only the infinitesimal 
classical contact transformations, if we tried to determine the finite contact transformation cor- 


responding to an arbitrary potential energy V(g) we should be led back to the phase-integral and for 
the quantum theory corrections to the W. K. B. method. 








AUGUST, 1942 


JOURNAL OF CHEMICAL PHYSICS 


VOLUME 10 


Statistical Distribution Laws for Rate Processes 


Bruce LONGTIN 
Illinois Institute of Technology, Chicago, Illinois 


(Received May 4, 1942) 


A re-examination of the fundamental hypotheses of statistical mechanics shows that the 
principle of the most probable distribution may be interpreted in such a way as to be applicable 
as well to systems which are not at equilibrium. The general method of deriving distribution 
laws for systems not at equilibrium is outlined. The use of this method is illustrated by deriving 
a general distribution law for systems in which the velocities are required to be non-Maxwellian. 
The distribution of velocities of molecules in any group having a specified average velocity is 
Maxwellian relative to the group average velocity. The average total energy of such a system 
is shown to be the sum of the kinetic energies of mass motion of the several groups of molecules 
plus the Maxwellian average energy corresponding to the given volume and temperature. 





ROSS physical properties of a system of 

molecules may be calculated from the mo- 
lecular constants either by kinetic theory meth- 
ods or by those of statistical mechanics. The 
kinetic theory method as applied by Enskog,! 
Chapman,”* and others consists in solving the 
Maxwell-Boltzmann integro-differential equation 
by successive approximations. The zeroth ap- 
proximation is Maxwell’s law‘ of velocity distri- 
bution. The second and third approximations 
have been obtained only with difficulty,’ and 
little work has been done on higher approxima- 
tions. The method at present is further limited 
to gases at low pressures. 

The method of statistical mechanics, whether 
based on the canonical ensemble of Gibbs® or the 
principle of maximum probability® is more direct, 
giving the form of the statistical distribution 
function free of approximations. However, this 
method has been developed only for systems at 
equilibrium. The limitation to systems at equi- 
librium is not necessary, and may be removed by 
re-examining the fundamental hypotheses. 

In developing a statistical mechanics which is 
applicable to rate processes, it is desirable to 
insure at the same time that it will apply as well 


1D. Enskog, Kinetische Theorie der Vorginge in mdssig 
verdiinnten Gasen (Dissertation, Upsala, 1917). 

2S. Chapman, Phil. Trans. Roy. Soc. A211, 433 (1912); 
216, 279 (1916); 217, 115 (1917). 

3S. Chapman and T. G. Cowling, The Mathematical 
Hoste of Non- Uniform Gases (Cambridge University Press, 

4 J. C. Maxwell, Phil. Trans. Roy. Soc. 157, 49 (1867). 

5 J. Willard Gibbs, The Collected Works (Longmans, Green 
and Company, New York, 1938), Vol. 2. 

6 J. H. Jeans, Phil. Mag. [6] 5, 597 (1903). 


to condensed systems. For this reason, the 
method of Jeans,® in which the element of the 
statistical ensemble is the whole system of mole- 
cules, will be followed. The language used is that 
of classical mechanics; however, only minor 
changes are required by the substitution of 
relativistic and wave mechanics. 


GENERALIZED STATISTICAL MECHANICS 
Definitions 


In establishing a law connecting the gross 
physical properties of a system, a series of con- 
nected measurements are made on the same sys- 
tem. The procedure is then repeated on the same 
or another practically identical system under 
comparable conditions. The number of repeti- 
tions is increased until the experimenter is con- 
fident that further repetitions will not materi- 
ally alter the average values of the measured 


properties. 
The requirement of comparable conditions im- 
plies that certain gross parameters g1, go, «-~ of 


the system are the same in each repetition of the 
series of measurements. This specification is in- 
adequate, as evidenced by the different values of 
the properties x, y, --- obtained in the several 
sets of measurements. A complete specification 
would describe the system as existing in the state 
i whose properties are precisely x;, yi, ---. The 
inadequate specification leaves the system to as- 
sume any one of the detailed states 1, 2, ---, 4, 
+++ by chance. 

The system may happen to assume the 7 state 
in a number nw, out of n repetitions. This defini- 
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tion of the relative frequencies w; is such that 
a wWi= 1, ; (1) 


while the average value of the property x; for this 
particular set of repetitions is 


t= } 2 XW; (2) 


If the property x is constant over the states 7, 
then x;=x;=--:=Z; the average value is the 
same as that for each state. In particular, >>; w; 
is the average value, 1, of unity. 

As the number » of repetitions is increased, 
w; may approach a limit denoted by W; if the 
limit exists. Then since x; depends only on the 
definition of the state i, the average value Z will 
approach a limit X given by Eq. 2 when W;, is 
substituted for w;. Let the operator D denote the 
value of the quantity operated on minus its 
limiting value. Thus Dw;=w;—W,, Di=z%—-X, 
etc. By applying this operator to Eqs. 1 and 2, 


one finds 
ae Dw;=D1i=0, (3) 


and 
Di= p x: Dw. (4) 


In the limit of large n, the operator D must be 
interpreted as denoting the change in the value 
operated on due to a trial variation in the values 
of the W;,’s. 


Fundamental Hypothesis 


In order to give any meaning to statistical 
methods, it is necessary to make some assump- 
tion which in effect requires that 

the set of limiting frequencies W; exists 
and is determined by the conditions im- 
posed on the set of repeated experiments. 
The relation between the values W; and the 
experimental conditions constitutes the statis- 
tical distribution law by means of which the 
limiting average values X, Y, of various 
physical properties may be calculated. 

Let a; denote the a priori probability that in a 
single experiment the system will assume the 
particular detailed state 7. In general this quan- 
tity is to be interpreted as the quantum weight 
of the state i, although other interpretations may 
be necessary in some cases. Further let 28,(w1, 





We, -- +) denote the probability that out of a series 
of m experiments the state 1 was assumed exactly 
nw, times, the state 2 in mw» times, etc. If the 
performance of n repeated experiments consti- 
tutes a series of unrelated events, then 


Hi ,,(wi, w2:--)=n! II (a"*/nw;!). (5) 


Using Stirling’s approximation, if the number n 
of experiments is very large, Eq. 5 may be 
written in the form 


log 78, = n(log a/w)w, (6) 
where (log a/w), represents 
> a w; log a;/wj. 


The hypothesis of a limiting distribution im- 
plies that in the limit of a very large number of 
experiments, the distribution w;=W,, we= Ws, 

- becomes very much (i.e., infinitely) more 
probable than any other distribution. For this to 
be true it is necessary that (log a/w)4 have a 
value for this distribution which is at least in- 
finitesimally greater than that for any other dis- 
tribution. If (log a/w), is a continuous function 
of the w,’s, then its value for the limiting distri- 
bution is a true maximum, and its variation 


} (log a;/W;—1)Dw;=0 (7) 


for all small values of the Dw,’s which are con- 
sistent with the experimental conditions. Ac- 
tually, (log a/w), is defined only for integral 
values of the numbers nw;, and may be considered 
a continuous function only in the limit of very 
large values of n. 

Equation 7 is analogous to the usual principle 
of maximum probability,*’7 but must be inter- 
preted somewhat differently. In the usual treat- 
ment the distribution law concerns the number 
of independent molecules out of a very large 
total number which are to be found in a particu- 
lar quantum state at equilibrium. The distribu- 
tion implied by Eq. 7 is not concerned with 
numbers of molecules, but rather with the num- 
ber of experiments out of a large number 2, in 
which the system as a whole is observed to as- 

7See, for example, E. H. Kennard, Kinetic Theory of 


Gases (McGraw-Hill Book Company, New York, 1938), 
Section 199. 
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sume a particular detailed state. The statistics of 
independent molecules may be considered as a 
special case in which the experimental system 
consists of only a single molecule. Furthermore, 
no assumption has been made which in any way 
limits the result to equilibrium states. 


Distribution Law 


To obtain the distribution law appropriate to 
any particular set of experimental conditions, 
Eq. 7 is to be solved subject to these conditions. 
The form of the resulting distribution law will 
depend on the nature of the restrictions which 
the system must satisfy. 

In obtaining the solution, it is only necessary 
to consider a set of independent restrictions. If 
one of the restrictions implies a second, only the 
first need be included; the second will necessarily 
be satisfied. By using this principle extensively, 
calculation of the distribution law is greatly 
simplified. 

The major features of the resulting distribution 
law are determined by equations which contain 
either the w,’s or the W,’s explicitly. Of these 
Eq. 7 is the only one which has been found to be 
generally true for systems at equilibrium. Any 
other such restriction involving the w,’s or W,’s 
explicitly can therefore be valid only under spe- 
cial conditions which are imposed by the ex- 
perimenter. 

The arbitrary restrictions imposed by the ex- 
perimenter may often be expressed in the form 


} B x,Dw; => 0, (8) 


where the property x; depends only on the state 
of the system and not explicitly on the frequen- 
cies w; or W;. For example, the rejection of data 
for which the average value ~ does not tend 
toward a definite limit leaves a set of data for 
which D&=0 in the limit; Eq. 8 is satisfied. This 
is simply the requirement that the property x; 
have a statistically reproducible value. 

As a working hypothesis it will be assumed 
that no arbitrary restriction containing the w,’s 
or W,’s explicitly need be considered which can- 
not be expressed in the form of Eq. 8. In all cases 
so far investigated, of the various possible ways 
of expressing a particular experimental condi- 
tion, the one most closely representing the phys- 
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ical facts has been in the form of Eq. 8. Since in 
practice such arbitrary restrictions can be exe- 
cuted only by rejecting data until the number of 
experiments falling in each detailed state has the 
specified value, it would also appear that any 
form other than Eq. 8 implies an unduly biased 
selection of data. The possible occurrence of such 
restrictions will therefore be left to be treated as 
special cases. 

The remaining restrictions do not involve the 
W,’s or w,’s explicitly. Hence the required solu- 
tion of Eq. 7 may be obtained by first eliminating 
the Dw,’s from an independent set of restrictions 
of the form of Eq. 8, together with Eq. 7. The 
resulting equation is then to be solved simul- 
taneously with the remaining restrictions. 

The result obtained in the first step by the 
method of Lagrange multipliers® * 9 is 


log W/ai=dAxxi tA it “ey +A1, (9) 


in which each coefficient \ must have the same 
value for all states of the system, but is otherwise 
as yet unrestricted. These coefficients may there- 
fore prove to be functions of the parameters g; 
whose values are independent of the detailed 
states of the system. The term ), enters in order 
to insure that Eq. 3 can be satisfied. Each of the 
other terms corresponds to one of the independ- 
ent equations, there being one term for each 
equation. 

Of the remaining restrictions, those which con- 
tain x,, y;:-- explicitly impose certain conditions 
on these variables. The others, which do not con- 
tain these variables, can enter only through the 
\’s. It is usually convenient to leave the distribu- 
tion law in the form of Eq. 9, remembering that 
additional restrictions remain to be satisfied. 


Systems at Equilibrium 


When the only experimental restriction is that 
the total energy E; of the system have a statis- 
tically reproducible value (i.e., DE=0 in the 
limit), then Eq. 9 reduces to 


log W;/ai=dzE;+ 1. (10) 


Equation 10 represents the Gibbs® canonical en- 


8 See, for example, R. B. Lindsay, Introduction to Physical 
Statistics (John Wiley and Sons, New York, 1941), p. 54. 
9J. E. Mayer and M. G. Mayer, Statistical Mechanics 
(John Wiley and Sons, New York, 1940), Chapters 3 to 5. 
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semble, in which the coefficient \z has the ex- 
perimental value? —1/kT, while \, has the value 
A/kT, where A is the Helmholtz free energy, k 
the Boltzmann constant, and JT the thermo- 
dynamic temperature. 

In applying Eq. 10 it is assumed that the num- 
ber of molecules Nz; of each component B has a 
definite value. If instead it is only assumed that 
the value is statistically reproducible, then 
DN x= 0 in the limit for each component. In this 
case, Eq. 10 also contains a term AgN x; for each 
component so restricted, and then represents the 
grand canonical ensemble® in which each \g is 
(1/kT) times the partial molal free energy (i.e., 
chemical potential) of the component. Other 
restrictions lead to other important forms of the 
ensemble.” 


NON-MAXWELLIAN VELOCITY DISTRIBUTIONS 


The most general problem of velocity distribu- 
tions is that in which it is required that each 
molecule / have a different statistically reproduci- 
ble average momentum p;; the restriction Dp, =0 
is satisfied by each molecule in the limit. Since 
these restrictions are independent of the require- 
ment DE=0, Eq. 10 must now contain an addi- 
tional term 2,-p; for each molecule. By the 
classical method it may be shown that the coeffi- 
cient Ag still has the value 1/k7. Expressing E; 
as the sum of a potential energy V; and the 
kinetic energy, the distribution law is found to be 


1 
log W;/ai= ——(Vit+> pi:?/2m,) 
kT l 


+> dr-puitrAr, (11) 
I 


where m, is the mass of molecule /. 
By writing w,; in place of mkT2:, Eq. 11 may 
be expressed in the form 


log W;/ai= — Vi/kRT—X(pui—wi)?/2mkT 
l 
+(1-> u.?/2mkT). (12) 
l 


The last term contains only parameters which 
are independent of the states 7, and hence is a 
new arbitrary constant 


b=i-D u?/2m kT. (13) 
l 


°R. H. Fowler and E. A. Guggenheim, Statistical Thermo- 
dynamics (University Press, Cambridge, 1939), Section 611. 


It is thus seen that log W;/a; is the same function 
of the (pi:—w.)’s and Ao for the non-Maxwellian 
distribution velocities that it would be of the 
Pi s and \, for a Maxwellian distribution at the 
same temperature. 

The summation over states of any function f; 
which contains the momenta and the constants 
wu, only in the differences (p;;—w,) will be in- 
variant to a change of the w,’s, provided only that 
the array of detailed states 7 is always trans- 
formed into itself by any shift of the origin of 
momenta in phase space. This requirement is 
satisfied a fortiori when the spectrum of per- 
missible velocities forms an unlimited continuum. 

In particular, from Eq. 1 it is seen that 


exp (—Yo) = ai 
Xexp [(-1/kT)(VitD (pii—wi)?/2m,)], (14) 
l 


which is such a summation. Since exp (— Ao) must 
therefore be invariant to a change in the w,’s, its 
value can be obtained by setting all of the w,’s 
equal to zero. It is then seen that Ao is equal to 
Ao/kT, where Ao is the value of the Helmholtz 
free energy for the system at the given tempera- 
ture and in the specified volume, for a Maxwellian 
distribution of velocities. 

Similarly, the value of (p;:— uy) must be in- 
variant to a change of the u,’s and hence equal to 
its Maxwellian value. This value is zero since 
(pii— ur\wis then the Maxwellian average momen- 
tum. In general (p;;— us) may be expanded as 
Pii— ui, since uw is constant over the states 17. 
Hence it has been shown that 


Pi =u; (15) 


the constants uw; must be equal to the specified 
average momentum. 

The distribution law may now be expressed in 
the complete form 


1 
log W:/a;=—[_Ao— Vi 
kT 


—> (pri—pr)?/2m,]. (16) 
1 


It is interesting to note that essentially Eq. 16 
was obtained by Maxwell,*:4 but for some reason 
the constant p, has always been interpreted as 








miVo, where Vo is the velocity of a system of refer- 
ence coordinates. For this reason the terms in p, 
have always been dropped. From the present 
derivation it is obvious that such a narrow inter- 
pretation was unwarranted. 

In any particular velocity distribution problem 
the requirements will probably take the form 
that the total momentum of molecules of a par- 
ticular group should have a statistically repro- 
ducible value. A detailed treatment of this case 
shows that the appropriate distribution law is 
obtained from Eq. 16 by setting each p; equal to 
the group average molecular momentum. Thus 
Eq. 16 represents the general distribution law for 
all problems of non-Maxwellian distribution of 
momenta, provided the array of detailed states 
i possesses the requisite transformation property. 

Several important features of this distribution 
law may be observed. First it is noted that the 
distribution of velocities of molecules about the 
specified group average velocity is Maxwellian. 
In particular, if all molecules in the system are to 
have the same group average velocity, the law 
represents a mass of molecules moving with this 
velocity whose distribution of velocities relative 
to the center of gravity of the mass is Max- 
wellian. Furthermore, in general the distribution 
law is invariant to any change in velocity of the 
reference coordinates. 

One application of this distribution law is 
illustrated by calculating the average energy of 
a system in which the velocity distribution is 
non-Maxwellian. Following the same procedure 
used to evaluate the w,’s, one may define a quan- 
tity [Vit do: (pii—pr)?/2mi Jw which is invariant 
to a change in the p;’s, and whose value is Eo, 
the Maxwellian average total energy of the sys- 
tem for the same temperature and volume. By 
expanding the separate terms of this average, the 
result, 


(Vitd (pii?/2m1) Inv 
2 (Pes Bo/mi)w tL (p.2/2mi)w=Eo (17) 


is obtained. The first term is the required aver- 
age energy E while the second has the value 
— >. p?/m,, and the third the value }>; p?/2mi. 
Making these substitutions, the average energy 
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E is found to have the value 
E=E +d pr2/2m. (18) 
l 


Equation 18 expresses the average total energy 
of the system as the sum of the Maxwellian total 
energy, which consists of potential and thermal 
energy, plus the kinetic energy of mass motion of 
groups of molecules. The correctness of this result 
may be seen from considering the special case of, 
for example, a planetary system, in which the 
total energy of the system consists of the Max- 
wellian total energies of the individual planets 
plus the kinetic energy of mass motion of each 
planet. 

This result is to be compared with the second 
approximation in the method of Enskog and 
Chapman,*”? which gives W£(1+})0. bi-pi:) for 
the distribution function. Here W;° is the Max- 
wellian distribution function, and b; is a small 
arbitrary constant. Although this form agrees 
with Eq. 16 to this degree of approximation, it 
requires E=E, which cannot be generally true. 
It is only with the inclusion of the third and 
higher approximations that the expression for the 
average energy begins to take the proper form. 


CONCLUSIONS 


The method of deriving statistical distribution 
laws by solving Eq. 7 subject to the expressed 
independent experimental restrictions is a simple 
extension of the method used for systems at 
equilibrium. It is, however, in no way limited to 
equilibrium conditions. While the application of 
this method has been illustrated only by the 
derivation of a distribution law for systems in 
which the velocities are required to be non- 
Maxwellian, it may also be used to solve other 
statistical problems encountered in the study of 
rate processes. 

It is of course recognized that the true test of 
this as well as of any theory is whether it agrees 
satisfactorily with experimental results, par- 
ticularly when extended into new fields. At 
present this test has been applied only qualita- 
tively, but with promising results. Consequently 
the author has no illusions that the distribution 
law for rate processes has been found in its 
ultimate form, but only hopes that a step has 
been taken in the right direction. 
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The methods developed in the previous paper are used to determine the most general form of 
statistical distribution law suitable for the treatment of rate processes. The result is particularly 
applicable to relatively dense systems composed of soft molecules. It is based on the assump- 
tions that each variable whose rate of change is to be measured must be statistically repro- 
ducible at all times, and that the molecular motions are not seriously discontinuous. Various 
special forms of the distribution law result from different ways of expressing the rate variables 
and experimental conditions. One form is derived for the case in which the rates of change of 
the average position and kinetic energies of the molecules are to be reproducible. Qualitatively, 
it agrees with the known frequent behavior of rate experiments, showing frictional dissipation 
of energy, and the tendency to return to equilibrium. Furthermore, it reduces to the classical 


distribution law under equilibrium conditions. 





ar the previous paper! a general method of 
deriving statistical distribution laws was de- 
veloped, and was applied to the case of a non- 
Maxwellian distribution of velocities. In this 
application, time was not considered explicitly 
as a variable. In the general study of rate 
processes, specific attention must be given to the 
way in which the experimental system changes 
with time. In so doing, the meaning of a statis- 
tical distribution and the form of the resulting 
distribution law must be somewhat modified. 


DISTRIBUTION OF RATE EXPERIMENTS 


The statistical method in general is concerned 
with the distribution of certain statistical ele- 
ments among several categories. In the study of 
a rate process, the statistical element is a rate 
experiment on a single system, in which the 
significant properties x, y, of the system 
changed in a particular way with time. The 
categories among which a series of rate experi- 
ments may be distributed are the different 
possible ways in which the set of properties 
x,y, +++ of the experimental system may have 
changed. Thus the statistical treatment of rate 
processes is concerned with the distribution of 
a series of repeated rate experiments among 
various detailed behaviors. 

When a particular experiment is said to belong 
to the detailed behavior 7 it is implied that the 





‘B. Longtin, J. Chem. Phys. 10, 546 (1942). 
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properties x, y, --- of the system have changed 
with time in the way ascribed to this detailed 
behavior. If the properties of the system have 
changed in some other way, then the particular 
rate experiment belongs to a different detailed 
behavior. 

In particular, the variation of the property x 
with time for the detailed behavior i may be 
represented by a Taylor’s series 


xi= QD (d™x,/dt™)o(t—to)”/m!, (1) 


at least up to the time at which the first discon- 
tinuity appears in any one of the derivatives. 
The quantities (d”x;/di")o are values of the 
various time derivatives of x; appropriate to the 
reference time tp. Any change in the numerical 
values of any of these quantities will cause x; to 
change in a different way with time. Thus a 
particular rate experiment belongs to the de- 
tailed behavior 7 if in that experiment all of the 
quantities (d"x/dt")o, (d™y/dt™)o, --- have the 
values appropriate to this detailed behavior. If 
any one of them is different from the value 
required of the detailed behavior 7, then the 
experiment belongs to another detailed behavior. 

In classical mechanics, if the coordinates and 
momenta of all elementary particles in a system 
are given specified values at the time fo, the 
behavior of the system is completely determined. 
These coordinates and momenta represent only 
a limited number of the (d™x/dt™)o’s; the rest 





are not independent variables and hence need not 
be specified separately. In quantum mechanics a 
somewhat different situation prevails. 

After all of the m repeated rate experiments 
have been completed and the data recorded, a 
statistical analysis is in order. At this time the 
number nw; of rate experiments belonging to 
each detailed behavior 7 is counted. It is at once 
evident that the value of mw; is not in any way 
determined by the length of time elapsed during 
any portion of an individual experiment, nor 
does it depend on time in any other way. Nor 
does the total number of experiments included 
in the statistical analysis depend on time. Thus 
each of the distribution frequencies w; must be 
independent of time. Likewise the limit W; 
approached by the frequency w; as the number 
is increased? must be independent of time. This 
constancy of the number of experiments belong- 
ing toa particular detailed behavior will be recog- 
nized as an analogue of Liouville’s theorem.*4 


Average Values 


The average value & of any property x; is 
defined as 

== x XxX {W;. (2) 

If the Taylor series expression given in Eq. 1 is 
valid, then 


B= ¥ (dxi/dé)ow,(t—to)™/m! 


i m=0 


=F (t—t)"/m! x (dmx;/de™)ow;:. (3) 


Since neither the w,’s nor the (d”x;/dt”)’s are 
functions of time, Eq. 3 expresses Z as a function 
of time by a Taylor series. Such an expression 
is valid only up to the first discontinuity. A more 
general form must be sought when the motion of 
the system is seriously discontinuous. 

One of the important requirements in a series 
of rate experiments is that each of the properties 
x, y, -*: whose rate of change is being studied 
must have a statistically reproducible value at 

2 Here one is comparing a set of statistical analyses each 
of which includes a different number n of rate experiments. 

3See R. H. Fowler, Statistical Mechanics (University 
Press, Cambridge, 1929), Section 1.5. 


4 J. W. Gibbs, The Collected Works (Longmans, Green 
andCompany, New York, 1928), Vol. 2, Part 1, Chapter 1. 
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all times during the progress of the experiment. 
This implies! that at all times 


Dz=0 (4) 


in the limit as m is increased, where D denotes 
the change in the quantity operated on, due to 
passage from the limiting distribution Wi, We, 
-++ to some other distribution wi, we, - 
If the Taylor series representation is used, 
DE=>Dd (t—to)™/m! > (d"x;/dt”™)oDw;, (5) 
since only w; is involved in a variation of the 
distribution. If Dz is to vanish at all times, the 
coefficients of each (t—f9)”"/m! must be zero; 


¥ (d"x;/dt") Dw; =0. (6) 


If Eq. 6 is satisfied for all integral values of m 
from zero up, then the quantity x will certainly 
be statistically reproducible at all times up to the 
first discontinuity. 


Distribution Law 


A statistical method is meaningless if the 
limiting frequencies W; do not exist. In order 
that they shall exist, it is necessary! that 


D(log a/w)15j =0 (7) 


in the limit, (log a/w), being the average value 
of log a;/w:, while a; is the a priori probability 
that the experimental system will exhibit the 
detailed behavior 7. The distribution law appro- 
priate to any particular case is to be found by 
solving Eq. 7 subject to the experimental re- 
strictions, and the requirement that the average 
value of 1 be unity. 

In the case of rate experiments these restric- 
tions include Eqs. 6. If none of the additional 
restrictions involves either the w,’s or W;’s ex- 
plicitly, the method of solution used in the 
previous paper! leads to the distribution law 


log W;/a:=A+ DY [Azo ™ (dx;/dt”) 6 
m=0 


+ryo™ (d"y,;/di™)o+---]. (8) 


The undetermined coefficients Azo, Azo, - °°; 
Ayo, Ayo, «++ have values appropriate to the 
choice of to as the time at which the derivatives 
d™x;/dt™ are to have reproducible values, these 
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coefficients being the same for all detailed be- 
haviors. 

Equation 8 expresses the frequency of occur- 
rence of rate experiments in which the various 
derivatives (d"x;/dt™) will happen to have any 
particular set of values (d"x;/dé™)» at the time fp. 
For example, if the quantities x; include molecu- 
lar velocities, Eq. 8 expresses the distribution of 
velocities, of accelerations, and of higher deriva- 
tives of velocity needed at this instant in order 
that measurement of the quantities x; shall give 
statistically reproducible results at all times. If 
the time ¢; had been chosen as the reference time 
rather than ¢, the distribution of velocities, etc., 
needed at the time /; would be given by exactly 
the same equation, , 


log W;/a; =A+ > [rr (d™x,;/di”) 1 


m=0 


+A yi (d"y;/dt") 1+ tis J; (9) 


except that the values of \,“ and (d”x;/dt”) are 
now those appropriate to the time 4. 

If the motion of the experimental system can 
be expressed by Eq. 1, then each of the deriva- 
tives (d"x,;/di")) may be expressed in terms of 
the (d"x;/dt”),’s by a Taylor series: 


(dmx, /dt") = ¥ (de+rx;/dt™) (to—ty)"/r!. (10) 
r=0 


By substituting these values into Eq. 8 and 
collecting terms containing the same derivative, 
it is transformed to 


log W;/a;x=A+>D, > (d™x;/dt™); 


z m=0 


s=m 


x DX Azo (to—b1)"—*/(m—s)!, (11) 
s=0 


where }°, implies summation over the different 
properties x, y, ---. 

Equation 11 expresses the frequency of occur- 
rence of any particular set of values (d”x;/dt”) 
at the time ¢, in a series of rate experiments 
which had been distributed at the time fo in the 
Way required by Eq. 8, and then allowed to 
change in their own natural way. If the distribu- 
tion law is to be more than a mathematical fig- 
ment, this result must be the same as that 
required by Eq. 9. This will be recognized as 
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analogous to the requirement in the Boltzmann- 
Enskog-Chapman' treatment of rate processes, 
that the number of molecules in any particular 
state of motion be stationary. However, it does 
not limit the possible applications of the method 
as severely. 

Remembering that the coefficients cannot be 
functions of the (d”x;/dt”)’s, it is evident that 
the coefficients of like derivatives in Eqs. 9 and 
11 must be identical. Hence 


Ari ™ = DY Azo’ (to— hs)” */(m—s)!. (12) 
s=0 
Equation 12 appears to express \z1°") as a func- 
tion of to as well as of t;. However, since neither 
log W; nor the (d”x;/dt”),’s depend on to, while 
the (d™x,/dt");’s are presumed to be linearly 
independent in order that Eqs. 6 may all be 
independent, it is evident that the A,;°’s can 
depend only on ¢; and not on tp. Hence it must 
be possible to combine the various terms in Eq. 
12 in such a way that fp vanishes. 
On differentiating Eq. 12 with respect to ¢; it 
is found that 


s=m—1 
An” /dti= — DY Azwo“(to—ts)™-**/(m—1—s)! 
s=0 
=n, (13) 


since \z9) is a function only of to and not of f;. 
On differentiating with respect to ¢) and collect- 
ing terms containing the same power of (fo—1/1), 
one finds 


Ade /Ato= ¥ (Azo +AAz0 /Alo) 


s=0 


X (to— ts)" —*(m—s)!. (14) 


In collecting the terms, two summations over s 
were combined in one of which the first term was 
missing ; this fact is expressed by defining 


Azo) =0. (15) 


In order that A,1° be independent of fo, it is 
necessary that 0A,1°"/dto be zero for all values 
of to. From Eq. 14 it is evident that this can be 
true only if 


dA20 /dto= —)Az9 8"). (16) 


5S. Chapman and T. G. Cowling, The Mathematical 
Theory of Non- Uniform Gases (The University Press, Cam- 
bridge, 1939). 
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Comparison of Eqs. 13 and 16, which are identi- 
cal in form shows that the requirements that 

(a) the A.0“”’s be independent of f1, 

(b) the A,1°’s be independent of fo, 

(c) log W; be independent of both to and f, 
and 

(d) that Eqs. 8 and 9 lead to the same result 
are self-consistent. 

From the foregoing it may be concluded that 
the general distribution law appropriate for any 
rate process is expressed by Eq. 8, in a form 
which is independent of the choice of the refer- 
ence temperature provided the undetermined 
coefficients satisfy Eqs. 14 and 15. This distribu- 
tion law is limited only in that it may not apply 
to processes in which the rate variables show 
serious discontinuities. For example, it will not 
apply to any molecular model assuming infinitely 
hard molecules, but will generally be applicable 
to any model in which, the molecules are more 
or less ‘‘soft.’’ Thus, while the Chapman-Enskog® 
method is most suited to the study of rarefied 
gases, the present method is more suited to the 
treatment of dense phases. 


SPECIAL FORMS 


Rate studies in general are concerned with the 
transport of material and energy, and with 
chemical reactions which may be expressed 
mathematically in terms of the change of certain 
reaction coordinates. The most detailed condi- 
tions of statistical reproducibility which might 
be needed to handle rate processes would thus 
require that the average position q., average 
momentum p,, and average kinetic energy 
(pi?)/2m, of each elementary particle / be re- 
producible at all times. Furthermore, if the 
experimental system is isolated, the total energy 
E;, will be a constant of the motion whose 
average value must be reproducible. 

So long as the momentum p); is expressible as 
mid qi;/dt, the requirement of a reproducible 
average momentum is implied by the require- 
ment of a reproducible average position, and 
hence is not an independent restriction. On the 
other hand, the average kinetic energy is not in 
general determined by the average momentum 

6 Cf. Glasstone, Laidler, and Eyring, The Theory of Rate 


oo (McGraw-Hill Book Company, New York, 
41). 
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(e.g., (pi?) is not equal to (pi)a?, so that the 
requirement of a reproducible average kinetic 
energy is independent. Thus the distribution law 
for this set of restrictions is 


log W;/a;=A+)k{E; 


+E EQ -d™qu/dem 
lL m=0 
+K)d"p,,?/di" }}, (17) 


the summation with respect to / extending over 
all elementary particles. 

The coefficients \zQ and \zK ,“” are the \’s 
corresponding to the m time derivatives of posi- 
tion and momentum. They must thus satisfy 
Eqs. 15 and 16 in order that the form of distribu- 
tion shall be invariant to a change in the refer- 
ence time. Equation 17 would contain terms in 
the derivatives d"£;/dt" but for the fact that 
these derivatives are all zero so long as E; is a 
constant of the motion. 

Less detailed conditions of reproducibility 
would require only that the average position, 
(© 1q1:/N)w, average total momentum, (>>: pii)a, 
and average total kinetic energy (>>: pi?/mi) 
of a group of N molecules be reproducible. A re- 
striction of this sort leads to the same result, 
except that Q,;, and likewise m,K,™, has the 
same value for all molecules of the group. 


Closed Form 


Equation 17 is not in a form convenient for 
practical calculations, since it contains an in- 
finite number of terms. It may be reduced to a 
closed form by assuming that for each of the 
detailed behaviors, the total momentum and 
total kinetic energy of each group are changing 
at a uniform rate. This assumption will not be 
too severe as a first approximation, provided the 
number of molecules in each group is very large. 

As a result of this assumption, the terms in 
dp ,2/dt” can be arranged in groups so that each 
group vanishes if m is 2 or greater. Similarly the 
terms in d”q,;/di” may be arranged in groups all 
of which vanish for any values of m equal to 3 or 
greater. Equation 17 then reduces to 


log W;/a;=A+Agk{E; 
+h CPr-qui—Pr- pri/mi+Qi: dis 


—K.ipi2/2mit+K pi: pii/mil}. (18) 
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In this equation it was assumed that the mo- 
mentum p;; can be expressed as m,dq./dt where 
m, is a generalized mass, while p:; is written for 
dpi:/dt. The coefficients Q,;, Q:, Q,@, K,, 
and K,® of Eq. 16 have been expressed as 
P,, —P., m.Q:, Ki/2mi, and K./2mi, respec- 
tively. 

Equations 15 and 16 imply that dz is inde- 
pendent of the reference time to. Since they also 
require that A4zQ: be —d0AcQ)/dto, and that 
hzQi® be —d0AcQ,/dto, it is further implied 
that 


P,=m dQ ./dbo, (19) 
and 
P,=dP,/dto, (20) 
while 
dP ,/dty=0. (21) 


In the above equations total derivatives are used 
since in Eqs. 15 and 16 the partial derivatives 
indicated simply that ¢; was being considered 
constant; it is not necessary to express this 
restriction when it is understood that the coeffi- 
cients do not depend on ¢; in any way. Finally 
it is implied that 


K.=dK 1/dto, (22) 
and 


dK j/dt)=0. (23) 


The particular form given for Eq. 18 was 
chosen to emphasize that P; must have the 
dimensions of momentum. Comparison with the 
results given in the previous paper show that at 
least in some cases P; is just the average par- 
ticle momentum for the group to which the 
particle / belongs. In general, due to the appear- 
ance of the terms in p;;- :;p, it may be just a little 
different from the average momentum; for con- 
venience call it the characteristic momentum of 
the group. 

Similarly Q; must have the dimensions of a 
coordinate (vector), while P,;/m, has the dimen- 
sions of an acceleration. Equations 19 and 20 
show that Q; and P,/m, are the coordinate and 
acceleration which correspond to the momentum 
P;. They may thus be referred to as the charac- 
teristic coordinate and acceleration of the group, 
respectively, and will sometimes but not always 
be equal to the average coordinate and ac- 
celeration. 
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In the classical distribution law dg is evalu- 
ated’ as —1/kT where k is the Boltzmann con- 
stant and T the absolute temperature of the 
system. In the case of the non-homogeneous sys- 
tems encountered in the study of rate processes, 
there may be some doubt as to the exact meaning 
of temperature. It is convenient, however, to 
write —1/kT° in place of \z, where T° is some 
temperature characterizing the conditions of the 
experiment, whose exact meaning remains to be 
discovered in each particular case. The coeffi- 
cients K, and K, have no previous analogue; 
the first has the dimensions of time while the 
second is dimensionless. 


Significance of the Terms 


In the distribution of rate experiments among 
detailed behaviors which is described by Eq. 18, 
those detailed behaviors will occur most fre- 
quently in which each term is large positive or 
small negative. Since \g is —1/kT°, and hence 
generally negative, each of the bracketed terms 
will tend toward a minimum value in the most 
frequent behavior. 

At equilibrium only the terms \ and AzgE; 
would appear in Eq. 18. Thus the added terms 
represent a displacement from the equilibrium 
conditions. Each term P;-q:; describes a displace- 
ment of the molecules from the equilibrium 
position, which will tend most frequently to 
make P;-q:; negative. Thus the most frequent 
displacement will be in a direction more or less 
opposite to the characteristic acceleration; the 
characteristic acceleration is in such a direction 
as to restore the molecules to their equilibrium 
position. The terms Q,-p,; describe the same 
fact in a different form, namely that the most 
frequent acceleration is in a direction opposite 
to the characteristic displacement. 

The terms —P,-p;;/m; were introduced and 
discussed in the previous paper. There it was 
shown that they represent a non-Maxwellian 
distribution of velocity in which each group of 
molecules has a Maxwellian distribution about 
their characteristic velocity. The termsin p;?/2m, 
may be combined with the kinetic energy part of 
E;, with the result that the kinetic energy of 
each group is multiplied by (1—K,1)/kT® instead 


7Cf. J. E. Mayer and M.G. Mayer, Statistical Mechanics 
(John Wiley and Sons, New York, 1940). 
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of 1/kT. In effect, each group of molecules now 
has a velocity distribution corresponding to a 
temperature 7°/(1—K,) instead of the tempera- 
ture 7°. Thus the terms —K.p:2/2m, allow for 
thermal inhomogeneities. 

The terms K ,pi;-pi;/m. describe the way in 
which the rate of conversion of potential into 
kinetic energy (i.e., the rate of production 
dp.2/midt of kinetic energy) departs from the 
normal. Thus they describe frictional and colli- 
sion phenomena where the previous terms de- 
scribed heterogeneity and driving forces. From 
this standpoint it is perhaps the term most im- 
portant in determining such properties as vis- 
cosity, and diffusional and electric resistances. 

Since the term tends as a whole to be negative, 
the most frequent rate of increase of kinetic 
energy must be opposite in sign to K;. When the 
most frequent value of the kinetic energy is less 
than the steady state value, it must tend to in- 
crease, while if greater it must tend to decrease. 
Thus the sign of K; depends on whether the most 
frequent value of the kinetic energy of the group 
is greater or less than the steady state value. 

Qualitatively, the distribution law described 
by Eq. 18 appears to agree in major detail with 
the known most frequent behavior of rate experi- 
ments. Quantitatively, it can be expected to be 
correct only insofar as the assumption of constant 
rate of change of total momentum and kinetic 
energy is justified. This assumption is certainly 
justified at least in the mean for systems in which 
the rate processes have reached a steady state. 
A more complete quantitative treatment is best 
handled in connection with specific problems. 
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Other Forms 


While Eqs. 17 and 18 are adequate to handle 
the majority of rate processes, they need not 
always be the most convenient. For example, 
the simple reaction 


H+H-—H; 


may be treated by requiring that the H—H 
distance have an average value which was 
initially large and changes reproducibly. Stated 
in this way, the statistical problem may be 
handled by means of Eq. 17 or 18. On the other 
hand it may be more convenient to require that 
the average number N»2 of molecules of H: be 
initially negligible, and increase with time in a 
reproducible way. In this case one must return 
to Eq. 8, including N2 as one of the variables 
Xi, Vi, +++ to obtain a distribution law which is 
an extension of Gibbs’ grand canonical ensemble.* 
The second of above two methods is essentially 
that used by Eyring® and co-workers in their 
theory of absolute reaction rates. However, their 
distribution law does not include terms beyond 
m=(0. 

At present it would be preferable in any par- 
ticular problem to start with Eq. 8. A considera- 
tion of various possible ways of expressing the 
experimental restrictions will lead to different 
forms of the distribution law among which 
certainly one or two should be conveniently 
usable. 


8See R. H. Fowler and E. A. Guggenheim, Statistical 
Thermodynamics (University Press, Cambridge, 1939), 
Chapter 6. 
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The Structure of Substituted Ethylenes and 
Their Isomerization Polymerization and 
“Peroxide Addition” Reactions 


ROBERT A. HARMAN AND HENRY EyYRING 
Princeton University, Princeton, New Jersey 
July 14, 1942 


HE double bond in substituted ethylenes can dis- 
sociate into a single bond in two ways.! 

(1) By a transition from its singlet ground state (in 
which the two electrons have opposite spins) to a triplet 
state (in which the spins are the same) about 25 kcal. 
above the ground state—the triplet mechanism. 

(2) By a rotational adiabatic transition to an upper 
singlet state about 40 kcal. above the ground state—the 
singlet mechanism. 

In both cases the excited products have the great reac- 
tivity of free radicals. In case (1) the activated product 
will be of an atomic type because the electron spins are the 
same: consequently the two electrons can never occupy the 
same orbital on one of the carbons to give an ionic type 
radical but will occupy two low orbitals having little ionic 
character. In case (2) the electron spins are opposite; the 
two electrons therefore occupy the same orbital which, 
depending upon the other substituents, may have appre- 
ciable ionic character. In spite of the much lower activation 
energy required for mechanism (1), actual systems may 
react sometimes by one mechanism and sometimes by the 
other because the probability of oan transition required for 
mechanism (1) is low.! 

On the basis of these considerations, it is possible to 
correlate a considerable number of observations involv- 
ing cis-trans-isomerization, polymerization, and “peroxide 
addition” reactions of ethylenic compounds. All of these 
reactions reveal the two types of mechanism in their 
behavior with catalysts and in the rate expressions for the 
reactions. 

It would be expected? that paramagnetic substances, 
free radicals, and atoms would catalyze mechanism (1) by 
providing a non-homogeneous magnetic field which will act 
differently on each of the two magnetic dipoles arising from 
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the spin of the two electrons of the double bond. An 
activated complex is thus formed with the paramagnetic 
substances involving the state with the spins the same, as 
well as the state with the spins opposed. The probability 
of the transition from singlet to triplet type of bond is thus 
increased without appreciable effect on the activation 
energy. 

That this expected effect of paramagnetic catalysts is 
actually found is demonstrated by the rate expression for 
initiating the polymerization of styrene (a) in an atmos- 
phere of nitrogen, (b) in the presence of benzoyl peroxide 
(which can react or decompose to give paramagnetic 
products), and (c) in the presence of free radicals.* 


(a) k’ 1.5104 e-28300/R7, 
(b) k’ 1010.85 7 23500/RT 
(c) k’ 1012-2 e@23500/RT | 


The effect of paramagnetic substances including atoms 
in catalyzing cis-trans-isomerization*® and “abnormal” 
addition‘ is likewise to be explained as a catalysis of the 
transition to the triplet state. The reversal of addition from 
the Markonikow rule when paramagnetic catalysts are used 
can be understood in view of the fact that the triplet 
activated complex is non-ionic as discussed above; con- 
sequently small differences in charge on the carbon atoms 
have no effect on the sense of addition of the incoming 
uncharged atoms or radicals, but that activated complex 
is favored which has the least steric hindrance and which 
has the greatest entropy. These effects will in general be 
secured when the heaviest entering group attacks the 
carbon with the fewest heavy groups on it (reverse of the 
Markonikow rule). 

The second mechanism, by contrast to the first, proceeds 


, by addition of a positively charged atom or group, and this 


fact accounts for catalysis of the second mechanism by the 
Friedel-Crafts’ type catalysts in isomerization,® poly- 
merizations,’? and addition‘ reactions. The Friedel-Crafts’ 
catalyst lowers the activation energy for forming ions by 
combining with the negative partner of the positive 
addend. The effects of resonance and induction which are 
found in substituted benzenes® operate in ethylene addition 
reactions by mechanism (2) to give the Markonikow addi- 
tion products for ortho-para-directing groups and the reverse 
addition for meta-directing groups. 

It is thus possible to outline a general theoretical rela- 
tionship between structures of substituted ethylenes and 
their reactivities. 

1J. L. Magee, W. Shand, and H. Eyring, J. Am. Chem. Soc. 63, 677 
OOF. Wigner, Zeits. f. physik. Chemie B23, 28 (1933). 


3 (a) G. V. — and E. Husemann, Zeits. f. physik. Chemie B43, 
385 (1939); (b) G. V. Schulz, H. Dinglinger, and E. Husemann, ibid. 


B39, 246 (1938); (c) G. V. Schulz, Naturwiss. 27, 659 (1939). 
4F, R. Mayo ‘and C. Wolling, Chem. Rev. 27, 351 (1940). 
5B. Tamamushi and H. Akiyama, Zeits. f. Elektrochemie 43, 156 
(1937). 
6C. C. Price and M. Meister, J. Am. Chem. Soc. 61, 1595 (1939). 
7H. Mark and R. Raff, High Polymeric Reactions (1941). 
8 T. Ri and H. Eyring, J. Chem. Phys. 8, 433 (1940). 





